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FOREWORD 


This volume (Volume II) of the proceedings of the GEOS-2 
Program Review Meeting held at NASA Goddard Space Flight Center 
on 22-24 June 1970 presents the results to date of the Goddard 
Space Flight Center tracking intercomparison tests conducted 
with GEOS-2. The volume is composed of a series of technical papers 
prepared by various investigators on the tracking intercomparisons 
conducted using the optical and electronic subsystems on the GEOS-2 
spacecraft . 


John Berbert 
GSFC 

Principal Investigator 



COMPARISON OF C-BAND, SECOR, AND TRANE T 
WITH A COLLOCATED LASER ON 35 TRACKS OF GEOS-2 


Prepared for the GEOS-2 
REVIEW CONFERENCE 
June 22 - 24, 1970 


By 


H. Parker 
D. Carney 
J . Berbert 



COMPARISON OF C-BAND, SECOR, AND TRANET 
WITH A COLLOCATED LASER ON 35 TRACKS OF GEOS-2 

H. Parker* D. Carney* J. Berbert 


As part of the GEOS Observation Systems Intercomparison Investigation at Goddard, 
several of the geodetic satellite tracking systems used with GEOS-2, including the GSFC 
experimental Laser (GOD LAS), an Army SECOR, and a Navy TRANET, were moved to the 
NASA Wallops Island station and located near the FPQ-6 and FPS-16 C-band radars there. 
Simultaneous tracking of GEOS-2 by all these systems was accomplished during April, May, 
and June 1968, to enable comparisons of the tracking data freed from the effects of un- 
certainties in survey, in the gravity field, and in systems time synchronization . Reference 
orbits were determined from the Laser data. Comparison of tracking data from the radio 
tracking systems with the 35 Laser reference orbits yielded residuals from which zero-set 
and timing biases were derived for each system . Preliminary results for the first 10 passes 
of simultaneous tracking of GEOS-2 were documented in X -514-68-458, "Comparison of 
C-Band, SECOR, and TRANET with a Collocated Laser on 10 tracks of GEOS-2", dated 
November 1968. These results were presented to the GEOS-2 C-Band Project Technical 
Conference at WTR in November 1968 and to the National Fall Meeting of the AGU in 
December 1968. 

Attempts to account for the various system anomalies reported in X-514-68-458 led 
to the following changes to the data sets. 

® The CDC-3200 Computer C-Band Preprocessing program was rewritten 
for the IBM 360 computer and expanded to include the computation of 
the calibration constants . 

© All the SECOR data were corrected by the Army Map Service for a scale 
factor error. 

• All the TRANET data were corrected by the Naval Weapons Laboratory 
for a Base Frequency Computation error . 

The effects of these changes to the data were analyzed for the same 10 GEOS-2 
passes and reported to the National Spring Meeting of the AGU in April 1969. Subsequently, all 
35 coincident Laser and radio tracking system passes were analyzed with the data modifi- 
cations included, and the results are the topic of this report. 
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The 35 Laser reference orbits were computed in the NONAME program using a priori 
measurement sigmas for relative data weighting of 2 meters for range and 200 arc seconds 
for both azimuth and elevation. The NONA ME program was developed for Goddard orbital 
studies and operates on the IBM 360 computer. Laser tracking data are summarized in 
Table 1, which gives the date, the coincident data, the time span and number of Laser range 
measurements, and the values of range and elevation at the start, end, and at the point of 
closest approach (PCA) for each of the 35 passes. 

The NONAME zero-set and timing bias recovery runs used tightly constrained 
orbital elements as determined by the Laser and assumed the a priori estimates of measure- 
ment and bias sigmas shown below. The a priori estimates of measurement sigmas were 
chosen to approximate the a posteriori e stimates resulting from the previous 10 pass studies. 
The a priori bias estimates were intentionally left unconstrained, except for the TRANET 
time bias, which was constrained due to the high correlation between the TRANET time and 
zero-set biases for these relatively short passes . 


MEASUREMENTS 


FPQ-6 RANGE 
FPQ-6 AZIMUTH 
FPQ-6 ELEVATION 
FPS-16 RANGE 
FPS-16 AZIMUTH 
FPS-16 ELEVATION 
SECOR RANGE 
TRANET RANGE RATE 


1.0 METER 

30 ARC SECONDS 
12 ARC SECONDS 
1.5 METERS 
50 ARC SECONDS 
20 ARC SECONDS 

2.0 METERS 
4 CM/SEC 


BIAS 

RANGE 
ANGLE 
RANGE RATE 
TRANET TIME 
TIME 


1000 METERS 
1000 ARC SECONDS 
1000 CM/SEC 
0.2 MILLISECOND 
1 SECOND 
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The results of the system intercomparisons , with anomalous data removed, are 
summarized in Table 2, which gives the sample size (number of coincident passes), and 
the average RMS noise, zero-set and timing biases as well as the total RMS of the residuals 
about the Laser orbits for each of the radio tracking systems . 

The time history of the derived zero-set biases over the duration of the experiment 
is given in Figure 1 for the range biases of each of the ranging systems and in Figure 2 for 
the TRANET range rate biases. 

The AN/FPQ-6 radar is a centroid tracker, and any difference in the pulsewidth 
used for calibration and the pulsewidth experienced in tracking will result in a range bias 
error. For GEOS-2 combined beacon/skin tracks this resulted in a 30 meter error in the 
beacon data. For these tracks, the skin track data near the center of the pass was removed 
and the beacon track data was corrected by adding 30 meters . The consistent SECOR range bias 
is larger than expected and still not explained . However, if the bias drift indicated in 
Figure 1 is extrapolated linearly to earlier times, the resulting bias would be near 0 meters 
just prior to launch . The TRANET data is submitted to the GSDS for both a low-frequency 
(162:324 MHz) and a high-frequency (324:972 MHz) pair. In this study the two frequency 
pairs were treated as two distinct stations . After the first several weeks of operation the 
variations in TRANET were markedly reduced. 

In summary, the results with all 35 of the WICE passes have the same general 
characteristics for the C-Band and SECOR data as reported earlier for the first 10 passes. 

The preprocessing changes noted in this report caused the C-Band bias values for a few 
passes to shift by several meters and caused the SECOR bias to be more negative by about 
2.5 meters on all passes. The preprocessing changes in the TRANET data reduced the 
previous biases from 15.6 and 17.6 cm/sec to -3.4 and 1.5 cm/sec. 
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Table 1 

this is a summary for station to52 
maximum elevation and minimum range measurements 
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Figure 1 
SER ORBITS 

-6, FPS-16 & SECOR Range Bias 











Figure 2 
LASER ORBITS 
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CARNARVON 

LASER COLLOCATION EXPERIMENT 
(CALACO) 

J. H. Berbert 
J. M. Hlavin 


INTRODUCTION 

In December, 1968, the Goddard mobile laser was shipped to Carnarvon, 
Australia, for the purpose of evaluating the accuracy of the Goddard Range and Range 
Rate (GRARR) System and the FPQ-6 C-band radar. An earlier laser/GRARR inter- 
comparison at the Rosman, North Carolina, station using GEOS-1 resulted in several 
recommendations for improving the GEOS-2 GRARR data. ^ This second laser/ 
GRARR test was needed to evaluate the improved GRARR data and to establish that the 
laser could support both day and night tracking at a remote site. At Carnarvon , the 
station code for the mobile laser was CRMLAS, and for the GRARR site it was 
CARVON. 

The laser arrived at the Carnarvon site in January, 1969. Setup operations 
were started immediately after arrival , and the first track was made on February 6, 
1969. 

The laser tracking system placed at Carnarvon uses an intense, highly colli- 
mated, short-duration light beam for illuminating the spacecraft being tracked. At the 
spacecraft, the beam is reflected back towards the ground station by an array of cubi- 
cal corner reflectors. The returning light is photoelectrically detected, and its travel 
time is measured to yield the range data. The data rate for the laser was 1 measure- 
ment per second in range , azimuth , and elevation. 

The GRARR system is a precision spacecraft-tracking system that determines 
range using sidetone ranging techniques. The range rate is determined using coherent 
Doppler principles. Each GRARR station uses an S-band system in conjunction with a 
1-, 2-, or 3-channel transponder on the spacecraft being tracked. The data rate for 
the GRARR was 1 sample per second in range and range rate. 

Data from the C-band radar system were not available for this report. 
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TRACKING DATA 


During the period from February 6 to May 16, 1969, 94 laser tracking passes 
were made. However, the data from only 84 of these passes were available for this 
report. The additional data from the other passes will be included in a later report. 

Of the 84 laser tracking passes investigated, 54 occurred during the day (sun 
elevation angle above -10°) and 30 occurred at night. These acquired data demon- 
strated that the laser could be regularly used for both day and night tracking when 
STADAN Minitrack data are used for initial spacecraft acquisition. 

There were 57 simultaneous GRARR/laser tracking passes during the test 
period. During the preprocessing and validation of the GRARR data, it was found that 
the range data taken April 29, 1969 , contained a systematic error due to an intermit- 
tent hardware failure. The range time interval decade counter which counted hundreds 
of nanoseconds failed in that it indicated a zero whenever a seven should have been in- 
dicated. Therefore , approximately 10% of the range measurements contained a 700- 
nanosecond error, equivalent to a 105-meter bias. The remaining 90% of the range 
data and all of the range rate data did not appear to be affected by the failure. How- 
ever, since there was a known hardware problem , data from the 13 GRARR passes 
after April 29, 1969, were discarded, and a set of 44 GRARR tracking passes were 
considered in this investigation. The analysis of the additional GRARR passes will be 
included in a later report. 

The GEOS-2 spacecraft has a 2-channel S-band transponder that is used in con- 
junction with the ground equipment. During the course of this experiment, both trans- 
ponder channels were used and evaluated; transponder channel A (1.4 MHz) was used 
for 12 GEOS-2 passes , and transponder channel C (3.2 MHz) was used for 32 GEOS-2 
passes. Channels A and C of the GRARR operate on up-link frequencies of 
2271. 9328 MHz and 2270. 1328 MHz, respectively. 

DATA ANALYSIS 

The raw laser data were preprocessed by the Optical Systems Branch at GSFC 
and deposited in the Goddard data center. The following systematic corrections were 

applied to the laser data. 

a. Range and elevation data were corrected for refraction. 

b. Range zero-set system delay correction was added. 

c. Outliers were removed at the five sigma level. 
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The GRARR station clock was synchronized to the cesium time standard , and 
the laser station time was maintained by the rubidium secondary time standard . The 
difference between these two standards was less than 0. 1 millisecond for the passes 
considered in this report. 

For data comparison purposes , it was assumed that there were no systematic 

errors in the laser range measurements. However , there may be a bias of up to 2 me- 

( 2 ) 

ters in the laser data ' ; , and this should be taken into account in evaluating the test 
results. 

The raw GBARR data were preprocessed with the following systematic correc- 
tions applied 

a. Range and range rate data were corrected for refraction. 

b. Range zero-set transponder delay correction was added. 

c. Range and range rate data were corrected for antenna offset. 

d. Range outliers were removed at the 10 sigma level. 

e. Station range zero-set bias was removed electronically by the station 
operator prior to each pass by matching survey and measured range. 

A careful check was maintained on the GEOS-2 transponder and GRARR ground 
station master oscillators. A small long-term (6 month) drift was noted in the on- 
board transponder oscillator. However, a preliminary analysis of the transponder os- 
cillator drift indicates that it would not be detectable in the measurements taken during 
this experiment. 

The procedure used to evaluate the GRARR data was to fit a short-arc orbit to 
the laser range , azimuth , and elevation measurements using a weighted least squares 
criterion. GRARR and laser residuals from the short-arc orbit were then computed 
and their statistics were tabulated. The GSFC NONAME ^ orbit determination pro- 
gram was used in this analysis. The a priori relative weight assigned to the orbital 
elements and the range , azimuth, and elevation measurements in forming the short- 
arc orbit are given in Table 1. The weights shown leave the orbital elements essen- 
tially unconstrained, allowing them to adjust to fit the laser data. 

The GRARR range and range rate, and laser range , azimuth, and elevation re- 
siduals about the laser reference orbit were computed for each short-arc. Because 
the GRARR data were edited such that they were within the laser time span, there was 
no extrapolation of the short-arc orbit used in determining the GRARR residuals. The 
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direct result of this procedure is to minimize any error contributed by the short-arc 

orbit. 

' Table 1 

Orbital Element and Laser Measurements Weight Factors 


Parameter 

Weight Factor 

cn 

i 

£ g 

Position (X, Y, Z) 

±10** meters 

w a 

Velocity (X, Y, Z) 

3 

±10 meters/second 

CO 

4-5 

Laser range 

±2 meters 

S3 

0) 

a 

Laser azimuth 

±200 arc seconds 

CD 

U 

0 

Laser elevation 

±200 arc seconds 

m 

Si 

CD 

GRARR range 

Not used in orbit formation 

S 

GRARR range and range rate 

Not used in orbit formation 


The mean and rms values of the residuals were calculated on a per pass basis. 
The rms values derived in the above runs are about the laser orbit and are indicative of 
the total error in the GRARR System (systematic errors plus noise) . The detailed per 
pass statistics of the range and range rate measurements are shown in Tables 2 and 3. 
The first column in Tables 2 and 3, EPOCH, gives the year, month, day, and hour of 
the pass under consideration. The second column, ID, indicates a day or night pass 
by D or N, respectively, and indicates the GRARR transponder channel , A or C, used 
when simultaneous data were taken. 

Average values and their variations for the rms and mean values of the range , 
range rate, azimuth, and elevation residuals with respect to the laser orbits for the 
44 GRARR passes and the 84 laser passes were calculated and are listed in Table 4. 

An examination of the laser data given in Table 4 indicates that there was no 
essential difference in the laser residual statistics between day and night tracking. 

An examination of the GRARR statistics in Table 4 indicates an average range 
total rms value of 7. 3 meters and an average range rate total rms value of 1. 7 centi- 
meters per second was established for the 44 passes. The variation in the range rate 
rms value between the two transponder channels is not considered significant. How- 
ever, a significant difference in range rms values between the two transponder chan- 
nels is indicated. It should again be noted that the average GRARR range and range 
rate rms values are from residuals about the laser orbit, contain systematic errors, 
and are representative of the total unmodeled uncertainty in the GRARR measurements. 
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Table 2 

NONAME Range Measurement Residual 
• Mean and rms Values in Meters . 

NONAME RANGE MEASUREMENT RESIDUAL MEAN AND RMS VALUES 

(VALUES IN METERS) 




• C R NLA S 

CARVON 

NCARVN 

• CRP4L AS 

CAR V3N 

FPOCH 

10^ 

ME AH 

MEAN 

MEAN 

RMS 

RMS 

69C20623 

0 

0@0 



I ©0 5 3 


69020709 

D 

O©002 



I©301 


69020723 

D 

# 0 



1 © 1 99 


69C2081 0 

D 

0® O 



1 © 343 


69021 51 0 

D 

0© 004 



1*773 


69021 61 1 

D 

d © o 



1 ® 235 


69C21 622 

0 

0 © GO 1 



I ©243 


69021 71 1 

D 

— G « 0 



0 ©893 


69021 723 

D 

0© GO 1 



1 ©589 


69021 B23 

DC 

0 ® 0 0 4 

-2*324/ 


1 @445 

5*595 

69C21 'n 0 

nc 

*“*0 © 0 

-9 © 1 1 8 


1 ©200 

9 @ 85 3 

69021 023 

DC 

■ *®0 ® 00 l 

“6 ©769 


1 ©41 1 

7*566 

69022010 

DC 

-0® 004 

1 ©743 


1 ©508 

3 ® 40 0 

690221 00 

D 

0© 0 



1 ©043 


69C221 l 0 

DC 

0*001 

—4 ® 455 


1 © 5 4 5 

5*295 

6 9C22323 

DC 

0 © 00 2 

-1 <®68 l 


I ©745 

3 *264 

6902251 C 

DC 

0 @ 0 

-1 ©239 


I ©070 

3*691 

69022523 

DC 

0® 00 2 

—5 ©0 75 


I @387 

7 @ 75 6 

6902261 0 

DC 

0 © 00 1 

0 ©698 


1® 354 

5® 25.2 

69022700 

DC 

0©0O! 

-10 ®I34 


1 ©0 81 

I 1 ©260 

6902271 0 

DC 

~*O©00! 

-7 #753 


1 ©40 1 

9 © 44 5 

6902280© 

DC 


-I6®©72 • 


1 © 2 0 3 

16©51 5 

69C228I 1 

DC 

0 @ ©0 1 • 

<3,, S ® 283 


1 ©1 16 

9® 736 

69030223 

DC 


—9 ©9 33 


. 0 @9 00 

10® 57 8 

6903031 0 

DC 

-0©O 

-I ®25© 


I ® 222 

3*524 

6903051 0 

D 

1 



1 ©627 


69030600 

DA 

0® © 

—3 ®657 


1 ©130 

5 © 005 

6903061 1 

D 

0© 005 



1 ©273 


6903091 0 

DC 

0©0 

-5 ©635 


@©953 

6 © 9 1 t 

69C31 01 0 

D 

-0*0 



I ©220 


69031 1 00 

DC 

-OeOOl 

© ®0 19 



5*600 

69031111 

DC 

-0®@©3 

-It ®387 


I ©4M 

12*512 

69031211 

NC 


-1S©839 


1 ©232 

1 6 © 349 

69031311 

NC , 

0® 003 

■ <*4 ©77 1 


1 ©217 

5 © 6© I 

69031 323 

DC 

© ® 0 

— 7 • 10 l 


1 ©0 4 1 

7® 649 

69C314I I 

NC “ 

®”0 ® 0 

-I2®944 


t ®S#4 

13® 453 

69C31 700 

DA 

— '0 ® 00 1 

-1 ®5t7 


t @274 

3©457 

69031 71 1 

NA 


—4 @961 


I © 540 

£©01 7 

69C31 80C 

DA 

0® 00 S 

-4 @671 


1® 130 

5® 729 

69031 81 I 

NA 

—0© 006 

© © S66 


1 © 5 40 

3 ©290 

69031 901 

DC 

0® 0 

0 ®7©3 


I ©107 

3*015 

69031911 

N 

0©G01 



1 ©249 


69031 923 

D 

0© 0 



1 © 2 1 2 


69032023 

DC 

™0 © © 

-ID ©637 


1 ©230 

11©I84 

69 €322 00 

DC 

-0® 001 

—8 ©710 


1 ® 60 4 

9*201 

6903241 1 

NC 

0© 0 

—4 ©9 © 6 


I © 46 1 

s* see 

€9032501 

DC 

— 0 © 0 

-7 ©762 


1 ©402 

8*272 

69C3251 1 

NC 

**■ 0 ® 0 

2 ©91 7 


1 #240 

4 © 6 1 3 

690331 1 1 

N 

0® 00 6 



1 ©416 


69C40I 12 

N 

-0® 002 



I ©8 74 


69040501 

D 

0 • ©0 2 



'1*376 


6904051 1 

N 

0 



& • 555 


69040601 

D 

O©0 



I #078 



NCAR VN 
RMS 
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Table 2 

NONAME Range Measurement Residual 
Mean and rms Values in Meters (Cow t ) 



NONA ME 

RANGE MEASUREMENT 

RESICU AL MEAN 

AND RMS 

VAL ues 




C VALUES 

IN METE RSI 





crmlas 

CARVON 

i 

NCARVN 

CRML AS 

CAR V3N 

EPOCH 

ID* 

MEAN 

MEAN 

MEAN 

RMS 

RMS 

6904061 1 

NC 

-0© 00 3 

-0 ©248 


1 ® 485 

2 © 63 C 

69040812 

NC 

-0© 0 

-10 ©325 


1 © 4£4 

i l © 449 

6904091 1 

NC 

0®'0 

— 12 ©870 


1*314 

1 3® 174 

6904! 000 

D 

0®0 



0 © 987 


69041011. 

~NC 

0®0©I - 

-8*074 


1*341- 

9® 173 

69C4I1 01 

0 

—0® 0 



1 ©0 57 


090421 00 

DA 

-0®0 

8 ©7 1 0 


1 « I 49 

9® 175 

690421 1 1 

NA 

0®0 

-5 ©8 25 


I ©086 

6 ©62 8 

090421 13 

NA 

0®0 

—3 ©776 


0 « 977 

4 « 7C 7 

69042211 

NA 

-0*00! 

-4 ©860 


! © 40 7 

5 © 67 7 

0904231 1 

NA 

—0 ® © 1 

—2 © 6 1 7 


1 ©072 

4© 162 

69042512 

NA 

-»0® ©04 

0 ©637 


1*366 

3 © 04 5 

69042001 

DA 

0®0 

1 ©169 


1 ©161 

3 © 29 9 

69042912 

N 

•0® 0 



I ©048 


6 % OA 30 01 

0 

~0®0 



1 ©019 


69043012 

N 

~0®001 


. 

0 ©98 1 


690501 02 

D 

—0® 0 



0 ®986 


090501 12 

N 

-0®0 



1®240 


09050212 

N 

*0® 00 1 - 



I ©300 


69050411 

N 

0®0 



1 ©073 


6905050! 

D 

-o®o 



I ©116 


69050512 

N 



4 

1 ©0 20 


69050601 

0 

-0®00t 

V 


I ©1 39 


69050612 

N 

0© 00 2 



1 ©0 27 


09051001 

0 

0®0 



1 ©168 


69051 1 12 

N 

0®0 



l © 383 


69 €51 2 01 

D 

-O©001 



1 « 1 90 


69051212 

N 

—0 © 0 



l ©285 


09051402 

D 

0®0 



1 ©162 


09,051 51 3 

N 

0®O 



1 a 557 


69051001 

D 

-0®001 



1 © 1 25 




CRMLAS^.. 

CARVON 

NCARVN 

CRML AS 

CAR VD N 



mu AN 

MEAN 

■ mm 

RMS 

RMS 

A VFRAGE 


0® ©0© 

-4 ©992 


1 © 258 

7*283 

ST® OS V® 


0® 002 

S © 194 


0 ©2 10 

3 © 652 

NO® OF PASSES 

84 

44 





* id ; time of pass; d-cay, n-night 

G PARR T RAMS PONDER CHANNEL - A 3H C 


NCAH VN 
RMS 


NCARVN 

RMS 
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Table 3 

NONAME Range Rate Measurement Residual 
Mean and rms Values in Centimeters per Second 


NCNAME R RATH MEASUREMENT RES IDUAL MS AN AN 0 RMS VALUES 
(VALUES IN CMS/S EC I 




CARVON 

C4RV0N 

EPOCH 

ID* 

MEAN 

RM S 

690 218 23 

DC 

l .60 5 

2.700 

6902 19 10 

OC 

0 .336 

1.149 

69021923 

DC 

0 .639 

1.403 

69022C 10 

DC 

C .836 

1.222 

69022110 

DC 

2.626 

2.905 

69022323 

cc 

0.367 

1.044 

69C22510 

cc 

1 .757 

2.097 

69022523 

DC 

0 .888 

1 .902 

69022610 

DC 

1.312 

2.458 

69022700 

CC 

0 .222 

1.020 

69022710 

DC 

C .466 

1 • 635 

69022800 

DC 

1 .287 

1.616 

6902 28 1 1 

DC 

0 .821 

1.786 

69 C 30 2 23 

CC 

0 .821 

1.810 

690 3 0 3 10 

CC 

1 .179 

1.973 

69C3C60C 

DA 

1 .40 2 

1.795 

690309 10 

DC 

0.338 

0 .953 

69031 100 

CC 

0 .30 4 

1.227 

69C 31111 

CC 

0 .731 

1.343 

69031211 

NC 

0 .8 1 4 

1 .46a 

69031311 

NC 

C .303 

1.225 

69C31 323 

CC 

0 .4 33 

1.187 

690 3 14 1 1 

NC 

0 .643 

1.062 

69031700 

DA 

0 .30 1 

1.173 

6903171 1 

NA 

-0 .229 

1.725 

69031800 

OA 

0 .287 

0.810 

69031811 

NA 

-0.118 

1 .245 

6903190 I 

DC 

0 • 107 

0 .866 

69032023 

OC 

0 .50 5 

C .997 

69032200 

CC 

-0 .639 

2.683 

690324 1 1 

NC 

1.443 

2.127 

6903250 1 

DC 

1 .659 

2.158 

690 3251 1 

NC 

-3.240 

4.731 

690406 1 1 

NC 

0 .865 

1.294 

690408 1 2 

NC 

0.146 

3.179 

6 9 C 4 0 9 1 1 

NC 

C .46 7 

1 • lo9 

69C4 1 0 1 1 

NC 

C .952 

1.420 

69C42 100 

DA 

-0 .205 

0.714 

69042111 

NA 

-C .CO 8 

1 .373 

69042113 

NA 

1 .486 

2.080 

69042211 

NA 

0.352 

1.740 

69042311 

N A 

0 .687 

1.514 

69042512 

NA 

0 .760 

2.235 

69042801 

OA 

0 .695 

1.410 



CARVON 
• MEAN 

CARVON 

RMS 

AVERAGE 

0 .600 

1.672 

ST. DEV. 

0.854 

0 . 749 

NO. OF PASSES 

44 



* ID 


TIME OF PASS: D-DAY * N-NIGHT 

GRARR TRANSPONDER CHANNEL - A OR C 



Table 4 

Summary of Laser and GRARR Statistics About Laser Orbit 


| Measurement 

84 Track Average 

54 Daylight 
Track Average 

30 Nighttime 
Track Average 


Range rms 

1.3± 0.2 m 

1.2 ± 0.2 m 

1. 3 ± 0. 2 m 


Range mean 

0. 0 ± 0.1 m 

0. 0 ± 0.1 m 

0. 0 ± 0.1m 

K 

m 

Azimuth rms 

31. ±17. arc sec 

32. ±16. arc sec 

29.0 ±18. arc sec 

m 

< 

Azimuth mean 

0. ±1. arc sec 

0. ± 1. arc sec 

0. ±1. arc sec 

1*3 

Elevation rms 

22. ±11. arc sec 

23. ±10. arc sec 

20. ±11. arc sec 


Elevation mean 

-9. ±14. arc sec 

-8. ±16. arc sec 

-10. ±11. arc sec 



Measurement 

44 Track Average 

Channel A (12 Track) 
Average 

Channel C (32 Track) 
Average 


Range rms (total) 

7.3± 3.7 m 

5.0± 1.8 m 

8. 0 ± 3. 8 m 

m 

Range mean 

-5.0± 5. 2 m 

-1. 7 ± 4.0 m 

-6. 1 ± 5. 0 m 

<3 

K 

Range rate rms (total) 

1.7± 0.8 cm/sec 

1.5± 0.5 cm/sec 

1. 7± 0. 8 cm/sec 

O 

Range ratemean 

0. 6± 0. 9 cm/sec 

0.5± 0.6 cm/sec 

0.7± 0. 9 cm/sec 


An inspection of the GRARR residuals plots indicates that small systematic 
errors existed within the GRARR data. An error model consisting of range and range 
rate zero-set biases and a common station timing bias was selected to explain the ob- 
served systematic error. The GRARR data were evaluated for a second time against 
the laser-determined short-arc orbits using the weighted least squares technique and 
the above error model to determine the true rms noise and bias values of the range and 
range rate measurements. The relative weights used in these regression runs are 
listed in Table 5. 


Table 5 

Weight of Orbit Elements , GRARR Measurements , and GRARR 
Biases Used in Regression Analysis 


Parameter 

Weight 

Orbital 

Ele- 

ments 

Position (X, Y, Z) 
Velocity (X, Y, Z) 

-3 

±1 x 10 meters 

±1 x 10 meters/second 

Meas- 

ure- 

ments 

GRARR range 
GRARR range rate 

±5 meters 

±0.05 meters/second 

Bias 

GRARR range 
GRARR range rate 
GRARR time 

±1000 meters 

±10 meters/second ! 

±1 second 
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The range rms noise, range bias, range ra;e rms noise, range rate bias, and 
station timing bias were computed for each pass. 

Average values and their valuations for the GllAItR noise and bias parameters 
determined in the regression runs for the 44 passes were calculated and are listed in 
Table 6. 


Table 6 

Regression Analysis Results 


GRARR 

Measurement 

44 Track Average 

Channel A (12 Track) 
Average 

Channel C (32 Tracks) 
Average 

Range rms 
Range bias 
Range rate rms 
Range rate bias 
Time bias 

2. 9 ± 0. 2 m 
-4. 3 ± 6. 0 m 
1. 3 ± 0. 6 cm/sec 
1. 2 ± 2. 3 cm/sec 
-0.4 ± 1. 1 msec 

3. 0 ± 0. 2 m 
-2. 0 i-4.0 m 
1.3 ±0.4 cm/sec 
0.0 ±0.7 cm/sec 
0.2 ±0.4 msec 

2.9 ±0.2 m 
-5. 2 ± 6. 5 m 
1.3 ± 0.6 cm/sec 
1. 7 ±2.5 cm/scc 
-0. 6 ± 1. 2 msec 


The results of the error model regression runs listed in Table 6 could account 
for the difference in the GRARR range rms values between the two channels. Here it 
is noted that the range noise rms values for both channels were the same (approxi- 
mately 3 meters). However , trie channel A transponder had a -2. 0-meter bias and 
the channel C transponder had i • 5.2-meter bias. The difference between the two 
transponder biases is probably due to drift in the transponder electronics since its 
prelaunch calibration. 

The regression values of average range rate rms noise given in Table 6 display 
no significant difference between the two transponder channels. 

Both the time bias average values in Table 6 and the variations about these 
values are unbelievably large. As noted earlier, a very tight control was maintained 
between the station clocks and on the GRARR system master oscillators. The results 
of this careful control should limit the pass-to-pass time biases to less than 0. 2 milli- 
second . It was noted that for these short-arc passes , the time biases were deter- 
mined primarily by the range data and thus may be the result of some other syste- 
matic error in the range data incorrectly modeled as a time bias. 

The range rate biases as determined in the regression show a substantial 
difference between the two transponder channels. However, it should be pointed out 
that there is a high correlation between time bias and range rate bias for the short-arc 
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geometry used in this intercomparison and that the extracted range rate bias for 
channel C may be the result of the unlikely channel C time bias. The negative average 
time bias shown is in the correct direction and has about the right magnitude to ac- 
count for most of the positive average range rate bias that is given. 

CONCLUSIONS 

It was concluded that the laser could support scheduled tracking, both day and 
night, at a remote site, and the laser measurement statistics show no significant dif- 
ference between day and night tracking. The basic laser measurement statistics were: 

® Laser range rms value . 1. 3 meters 

• Laser azimuth rms value 31. arc seconds 

© Laser elevation rms values 22. arc seconds 

In addition, it was concluded that the Carnarvon GRARR data had the following 

uncertainties: 

® Total range uncertainty of ±7. 3 meters consisting of: 

— Zero-set bias of -2.0 meters for channel A 
— Zero-set bias of -5. 2 meters for channel C 
— Measurement noise level of 2. 9 meters 

— Plus other systematic errors (such as time bias, refraction, etc.). 

• Total range rate uncertainty of ±1. 7 cm/ sec consisting of: 

— Measurement noise level of 1. 3 cm/sec 

— Systematic errors (such as zero-set bias, time bias, refraction, etc.). 

The above results are in substantial agreement with the earlier Rosman test, 
where the Rosman GRARR showed a range bias of -5. 3 ± 12.4 meters , a systematic 
S-shaped variation in the range rate residuals between ±6 cm/sec, and a 2. 1-ms time 

bias. 

The range measurement shows a small negative bias in both tests , with varia- 
tions in the Carnarvon basis being substantially smaller. The systematic component 
of the range rate residuals and the time bias errors determined by the regression runs 
are also substantially smaller for the Carnarvon test than those determined at Rosman. 
The S-shaped variation in the Carnarvon range rate residuals was absent in most of the 
passes and, when present, the effect was smaller than that documented in the Rosman 
test. The above results were probably due to the better control exerted over the sta- 
tion clock and over the transponder down-link oscillator as recommended in the Rosman 
test result's. 
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The above results augmented with the additional laser, GRARR, and the 
Carnarvon C-band data will be published in more detail in a future report. 
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LASER/LASER COMPARISONS FROM THE 
GORF, ARLACO, AND GOLACO TESTS 

J. Berbert, M. Hlavin, R. Reich 


I NTRODUCTION 

The Goddard Experimental Laser (GEL) has been used since 1966 for participa- 
tion in geodetic data observation campaigns and for intercomparison with other geodetic 
tracking systems (References 1,2, and 3). In late 1968, the Goddard Mobile Laser 
(GML) became available. This afforded an opportunity to collocate and intercompare 
the GEL and the GML to help determine whether systematic errors were apparent in 
the data from either one. Since the GML was developed, it has participated in three 
such laser/laser intercomparison tests: twice with the GEL at Goddard and once with 
the SAO laser at Mt. Hopkins, Arizona (see Figure 1). 

This paper summarizes preliminary results from all three tests. The GML 
data were used for forming the reference short-arc orbits, since the GML data were 
common to all three tests. The resulting reference orbits were tightly constrained 
and used to determine SAO or GEL measurement residuals from which the relative 
range biases and bias uncertainties were determined (see Figure 2). The a priori 
weights of the orbital elements and the GML range, azimuth, and elevation measure- 
ments, used in forming the short-arc orbits, are given in Figure 3. The a priori GEL 
measurement and bias weights are given in Figure 4. 

In these tests, relative time biases were not derived because the laser system 
clocks were known to be synchronized to within 0.05 millisecond. 

RESULTS 

A chronological plot of range bias of the GEL or SAO laser, relative to the GML, 
is plotted for each pass of the three intercomparison tests: GORF 1 (see Figure 5) , 
ARLACO - Phases 1 and 2 (see Figures 6 and 7), and GOLACO (see Figure 8). A sum- 
mary of the average range noise and average relative bias over all passes within each 
intercomparison test is given in Figure 9. 

GORF-1 

The large GEL rms noise value of 1.86 meters in the first laser intercomparison 
test (GORF-1) was attributed to using an older refurbished ruby rod at least an inch 
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shorter than the optimum 6 Inches. The shorter rod results in a longer pulse , which 
Increases the ms noise. All the GEL 6 -inch ruby rods had been provided to the GML. 


The large and relatively stable 4. 09-meter bias between the GEL and the GML 
exhibited in the five GORF-1 passes was an order of magnitude larger than had been 
expected. A review of the calibration techniques revealed that the GML bore sight 
tower calibrations had been done near the threshold levels of the return pulse energy; 
whereas , the GEL calibrations were done near the expected return-pulse signal level 
from the satellite. This would cause the GML system to trigger earlier on the higher- 
energy return pulses from the satellite than on the calibration pulses, thus leading to 



For this test, preparations were made to analyze the results of the first few 
passes as soon as possible in order to determine and correct any deficiencies. These 
precautions proved worthwhile , since quick-look analysis of the first pass data uncov- 
ered a 100-millisecond time bias and a 5. 5-meter range bias between the two lasers 
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(not shown in Figure 6). The time bias was determined to be due to an intentional 
offset in the time tag of the SAO data which was overlooked in the preprocessing of the 
data. A large part of the range bias, 4.8 meters, was due to a change in the SAO 
system internal delay since the previous calibration. These errors were detected, 
verified, and corrected by the end of the second pass. 

At the conclusion of the ARLACO test, SAO preprocessed their laser data and 
submitted them to the National Space Science Data Center (NSSDC) at GSFC These 
are the data denoted by the crosses in Figures 6 and 7 and averaged in Figure 9. No 
satisfactory explanation is presently available for the apparent 3. 3 -meter shift in bias 
between phases 1 and 2, or for the large outliers near the ends of both phases. 

The data denoted by crosses within circles in Figures 6 and 7, and also aver- 
aged in Figure 9, are the SAO data that were corrected for refraction. These data 
had to be corrected, since discussions with SAO indicated that the data in the NSSDC 
were not corrected for refraction (contrary to the refraction code indicator in the data 
format) . Refraction corrected data subsequently were submitted by SAO to the NSSDC 
in June 1970, after these analyses were done. 

The data denoted by triangles in these figures are the refraction-corrected SAO 
data compared with GML data using a constant 102. 2-nanosecond internal-delay calibra- 

f 

tion rather than the delay value measured for each pass. The 102.2-nanosecond value 
is the average measured GML internal delay for all the GOLACO passes, applied after 
the fact to the GML data from the ARLACO passes. As can be seen, this constant value 
increases the fluctuations in the bias for the first few passes of phase 1 and decreases 
the fluctuations for the last few passes of phase 2. The motivation for using a constant 
value of internal delay for the GML data arises from the GOLACO test results which 
are explained below. 

The decrease in the GML (MOBLAS) rms noise from 1.23 meters in the GORF-1 
test to 0. 99 to 1. 06 meters in the ARLACO and GOLACO tests is thought to be due 
to correction of an instability in the receiver circuit shortly after the GORF-1 test. 

GOLACO 

The results of the third laser inter comparison test (GOLACO) arc given in 
Figures 8 and 9. 
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The decrease in the GEL (GODLAS) ms noise from 1. 86 meters in the GORF-1 
test to 1.00 meter in the GOLACO test is attributed to replacing the shorter refurbished 
ruby rod with an optimum length 6-inch rod and to various other circuit improvements. 

It is interesting to note what appears to be a correlation between bias results 
and the training of the operator responsible for the GML bore sight tower calibration, 
indicating the somewhat subjective nature of the calibration. For the first three 
GOLACO passes, the GML calibrations were done by a member of the crew who had 
not previously had that responsibility. On about April 5th, the crew was retrained in 
the calibration procedures that were established at the conclusion of the GORF-1 test. 
This coincided with a noticeable decrease in the bias for the next several weeks. From 
about April 25th through the end of the test, the newly trained calibration operator was 
replaced by another operator. This coincided with a noticeable increase in the bias and 
bias fluctuations. 

A related but more specific effect is illustrated in Figure 10. Here , the derived 
range bias between the lasers is plotted against the GML internal delay measured prior 
to each pass during the bore sight tower calibration. The systematic nature of the bias , 
when plotted in this manner, indicates the high correlation between the range bias and 
the calibration measurement of internal system delay. In fact, a straight line with a 

t 

slope of 1. 5 meters per 10 nanoseconds (1. 0 meter per 1. 0 meter) and an intercept of 
106.2 nanoseconds provides a good least squares fit to the data in Figure 10. The rms 
of the bias variations about this line is only 0. 71 meter compared to an rms of 1. 25 
meters about the best fitting horizontal line. Thus , using a constant internal delay of 
106.2 nanoseconds, rather than the measured internal delay, reduces the bias 
variations. 

A reasonable physical interpretation of this result is that the GML internal 
system delay is more stable than the calibration measurements of it. For example , 
an increase of say 10 nanoseconds (1. 5 meters) in the GML internal-delay measure- 
ment (not accompanied by a real change in GML internal delay) would produce GML 
range measurements to the satellite shorter by 1. 5 meters. Consequently, the bias of 
the GEL relative to the GML would be larger by 1. 5 meters. Difficulties in matching 
the calibration return-pulse signal level to the satellite return-pulse signal level could 
cause such an effect. 

If the 106. 2-nanosecond intercept of the fitted line is chosen as the constant 
value for the GML internal system delay, then the average bias for the GOLACO data 
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is zero, with variations in bias of 0.71 meter, as shown by the variations about the 
fitted line. This 106.2-nanosecond value of internal system delay probably corresponds 
to the delay for a calibration return-pulse signal level equal to the average satellite 
return-pulse signal level for the GOLACO data. However, the GML system was not 
able to precisely measure or correct for return-pulse signal level during the GOLACO 
test. The return-pulse level, used in each prepass calibration, could be adjusted only 
in 10-db steps and at best was only an estimate of the expected average satellite return 
for that pass. Furthermore, considerable signal strength fluctuations are to be ex- 
pected from pulse to pulse within a pass. Because of this, the return signal should be 
measured or compensated for on every pulse. 

In Figure 11, the same range bias values are plotted against the GEL internal- 
delay measurements. In this case, the points appear randomly scattered, so that 
I correlation between these parameters is not apparent. This is probably due to the use 
of an experimental GEL device to measure the return-pulse level of each pulse during 
GOLACO passes. Using this device, a prepass boresight tower calibration of delay 
versus pulse level was applied to the measured satellite return pulses to correct for 
systematic errors in delay due to pulse level. 

| CONCLUSIONS 

i 

1. The relative bias between the GEL and the GML system was reduced from 4. 09 
to -1.19 meters from the GORF-1 to the GOLACO test. The rms noise of the 
GEL was reduced from 1. 86 meters to 1. 00 meter, and the rms noise of the 
GML was reduced from 1. 23 to 1.06 meters. These improvements are attrib- 
uted to several factors: 

• The better GML calibration procedure established after the GORF-1 test, 
which approximately matches the calibration return-pulse signal level to 
the average satellite return-pulse signal level. 

® The GEL pulse-level measurement and correction techniques. 

• The Improved system hardware in both the GEL and the GML (receiver 
upgrading, longer ruby rod, etc.) 

2. In ARLACO phases 1 and 2, the SAO system exhibited average pulse -to -pulse 
rms range measurement noises of 1.4 and 1.1 meters and average biases re- 
lative to the GML of -1.4 meters and +1.8 meters. 
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All 3 independent laser systems , involved in the ARLACO and GOLACO tests , 
were in agreement to better than 2 meters. This supports the hypothesis that 
these systems , at the time of the tests , had an absolute accuracy of better than 
2 meters. Furthermore, all three laser systems are undergoing signal proc- 
essing developments which are anticipated to improve their accuracy by almost 
an order of magnitude. 



K 

W _ 
ft ~ 
X w 

S 3 


tS » 

w ^ 
03 « 


00 

so 

\ 

rH 

O 

5 

w 

HH 

Eh 

< 

O 

O 

h-H 

CT> 

SO 

N 

t>* 

O 

Ph 

0 

c— 

N 

rH 

3 

w 

H 

O 

So 

M 

0 

t> 

\ 

\ 

rH 

t— 1 

ffl 

W 

*4 

\ 

rH 


W 

w 

rH 

rH 

O 

4 

rH 

rH 

O 

Q 

N 

O 

S 

0 

O 

O 

S 

So 

10 

Eh 

00 

Q 

W 

< 

Eh 

o> 

Eh 

cn 

p 

Ph 

Q 

Ph 

O 

Eh 

so 

\ 

< 

fcS 

O 

so 

\ 

so 

\ 

<3 

< 

0 

t- 

00 

Q 

w 

N 

so 

(M 

Q 

Q 

\ 

0 

Q 

0 

rH 

Q 

Q 

co 

\ 

0 

O 

« 

\ 

0 

\ 

C<J 

O 

O 

<N 

\ 

rH 

O 

<: 

rH 

rH 

O 

a 

00 


7 


33 


Figure 1. Laser/Laser Tests 
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Figure 2. Range Measurement Intercomparison 
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Figure 3. A priori Weighting of GML Data Used in Formini 
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Figure 4. A Priori Weighting of SAO and GEL Data Used in Bias Recovery 
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Figure 8. GOLACO, Range Bias vs Date, March— May 1969 






















Figure 11. Plot of Range Bias Between GEL and GML 
vs GEL Internal Delay 
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REFRACTION STUDIES USING THE WICE DATA 


J. Mallinckrodt, H, Parker, J. Berber! 

During the first few months of GEOS -2, Goddard sponsored the Wallops 
Island Collocation Experiment (WICE) as part of the GEOS Observation 
Systems Infcereomparison Investigation. The availability in the WICE 
experiment, of © large number of collocated redundant metric sensors 
experiencing predictably different degrees of refraction affords a unique 
opportunity to compare various methods of refraction correction. This 
paper will discuss the progress of refraction studies based on the WICE 
data. 


Figure I indicates the principal objectives of the study. The first 
objective involves formula comparisons. Since most practical refraction 
correction formulas involve varying degrees of approximation in themselves, 
irrespective of any errors in the external data inputs, a considerable 
effort has been made first to compare various refraction correction 
formulations under idealized common input data assumptions in order to 
become quite aware of their differences and to spot any gross deficiencies 
in some of the formulations. Within the GEOS program the refraction 
correction formulas in use by various agencies for correcting various 
data sources are not standardised. For the troposphere alone 4 such 
distinct formulations have been identified for elevation angle, 9 for 
range, and 7 for range rate. In addition to these, the comparisons 
have included the CRPL Standard Atmosphere calculations and a highly 
refined ray tracing integral program. The discrepancies have been found 
to be highly significant at low elevation angles and significant in some 
cases at all elevation angles. Figure 2 is presented to give am idea of 
the variability between various formulations . The report on this study 
has been written and is available (Ref. 1) . 

The second main source of variation in the refraction corrections 
is in the data source. For example, in the case of the ionosphere, 
options include ionospheric soundings , multi-frequency measurements such 
as with SECOR and TRANET, or a priori predictions . Figure 3 indicates 
the WICE data sources permitting refraction comparisons for the tropospheric 
and ionospheric errors . 

Recent effort has concentrated on the ionospheric comparisons and 
includes the following parts , 


Ref. 1: Herbert, J. BE. and Parker , H. C. , "GEOS Satellite Tracking 

Corrections for Refraction in 
the Troposphere" , GSFC X-514-70-55, 
February, 1970. 
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1* Composite Ionospheric Profiles 


- 2 - 


Figure 4 illustrates the problem o£ synthesizing a composite profile 
from available bottomside and topside (ALOUEIIE) sounder data* Coincidences 
in time and position between ALOUEIIE and GEOS were* of course 3 dependent 
on their different orbital parameters* Nevertheless , some 46 near 
coincidences during April-June, 1968 were identified* These were ranked 
quantitatively for degree of coincidence by evaluating a semi-empirical 
expression for ionospheric correlation as a function of spatial and 
temporal separation and from this and consideration of the availability 
of other data such as TRANE! , SECOR and Laser J^or Intercomparisons , the 
list was reduced to 17 candidate passes for extensive analysis* The 
ALOUEIIE data have been obtained for most of these cases , scaled and 
reduced to electron density profiles vs true height, as have the coincident 
Wallops Island bottomside sounder records* From these, and using appropriate 
interpolation between records , a composite best-estimate electron density 
profile has been formulated applicable to the point of closest approach 
on the GEOS trajectory for 16 of the 17 selected passes* 

2® Comparisons Between Profiles and 2-Frequency SECOR Corrections 

Ray tracing calculations have been carried out on the 16 composite 
profiles using the REEK program and implicitly assuming a spherically 
stratified ionosphere with no variation in electron density with latitude 
or longitude* These results have been compared with the two-frequency 
SECOR derived ionospheric corrections '(see Figures 5 and 6) * 

The SECOR correction at the maximum elevation point generally agrees 
within about 107* with the ray traced composite profile correction* However, 
initial results of these comparisons have shown that at times there is a 
considerably greater degree of North-South or rise-set assymetry than had 

been .anticipated* 

In Figure 5 the SECOR determinations show as much as 100% difference 
in range corrections at 35° elevation on the rise and set* Two possible 
explanations suggest themselves for this phenomenon* 

• Actual horizontal gradients of total electron content in the GEOS-2 
orbital plane* Sunset line effects seem the most likely explanation 
although preliminary studies indicate that even these would not 
normally be expected to induce more than a 30-40% change in electron 
content in a 700 mile span* 

* Magneto- Ionic effects® However, preliminary estimates do not 
explain the observed magnitude of assymetry® 
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3 • Comparisons Between SECOR and TRANET Multi -Fr eq uency Corrections 

Comparisons of differentiated SECOR with TRANET corrections and 
integrated TRANET with SECOR corrections have been made on several passes 
with excellent agreement* 

4 * Geographical Variation Predictions 

In pursuit of this question several approaches to predicting spatial 
variations or actual electron content horizontal gradients are being 
studied,, All of these utilize the published "Ionospheric Predictions 91 
or actually the harmonic expansion coefficient representation thereof 
as a basis for spatial scaling according to- the computed latitude and 
longitude of each height point along the line-of -sight ray to each of 
several representative points along the GEOS trajectory* 

5 * ALOUETTE Determination of Spatial Variations 

Additionally, the ALOUETTE topside reductions are being carried 
out at several points along the trajectory, representing roughly the rise 
to set trajectory above 30° elevation* This will permit a rough check 
of the gradients predicted as described above* If the discrepancies 
warrant , it may be advisable to attempt estimation of the electron content 
gradients on the basis of the ALOUETTE observations, however, this is 
still only being considered on. a contingency basis* 

In summary, the preliminary results show that, as expected, the MBS 
predicted ionosphere is unrealistic for determining tracking systems 
refraction corrections in the VHF/UHF band. Using a composite electron 
density profile, based on bottoms ide and topside ionograms, for determining 
the SECOR ionospheric refraction corrections, provides agreement with 
the 2-f requency SECOR refraction corrections to within XG7o at high 
elevation angles near the sounding station, but fails at low elevation 
angles, far removed from the area where the soundings were obtained,, 
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evaluation and as a reference for analysis of less precise tracking systems. 

Both functions were accomplished by forming precision orbits, the arc lengths 
being 2-4 days. The scheduling of the GRARR was oriented to this end throughout 
the GEOS-n lifetime. 


Having formed the orbit, data from the several GRARR stations were then evaluated 
and any inconsistencies noted. In those cases where operational problems were 


Bg jTjj i* Kira 


the orbit determination process was simply used as a tool to ensure that the 4 



Section 3.1*1. The additional corrections for GEOS-B are as follows: 
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1. - X-Y axis offset corrections to measured range and range rate 

2. Tropospheric and Ionospheric Refraction corrections to range and 
range rate (ref. Private communication; Parker to Herbert, Fall 1969), 
based on orbit determined elevation 

3. Range correction based on measured range rate and prelaunch 
calibration of spacecraft transponder 

4. Time synchronization corrections (as received) . 

3.2 GRARR DATA REDUCTION 

The orbit determination program used edits the tracking data based on the value of 
the measurement residuals from the orbit. The editing criteria employed are 90 
meters for range, 25 cra/see for range rate and 1 degree for X-Y angle. The final 
orbits contain from 8-20 spacecraft passes. 

3.3 QUALITY ASSURANCE RESULTS 

Several GRARR data errors were discovered and corrected as a result of the GEOS-II 
quality assurance program . Some of these are listed below . 

On approximately one half of the supported passes, the ALASKA GRARR exhibited 
a one second time error in the range rate information . The cause was found to be 
in the equipment turn on procedure. 

The Carnarvon station exhibited an abnormally large noise content in the range 
data when compared to other stations . The cause was found to be a faulty circuit 
board in one of the mixers; the symptoms could only be detected on targets possessing 
satellite dynamics. 

A procedural error in equipment turnoff was discovered at Rosman. 

Several cases of X-Y angle encoder slippage were detected at Carnarvon. A severe 
RF axis misalignment problem was also discovered at this station. 

3.4 GRARR VALIDATION STATUS 

Figure 2 contains the current status of the validation effort. The operational 
dates of each station are shown in Figure 3 . 
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MOTS OPTICAL DATA VALIDATION STATUS 
TECHNIQUE USED GSFC X-5 14-69-83 


GEOS-I MOTS OPTICAL VALIDATION REPORT 
DATA COLLECTION PERIOD FEB 1968 - DEC 1969 

NUMBER OF PLATES TAKEN 7971 

NUMBER OF PREOP ERATIONAL PLATES 41 

PLATES BROKEN OR NON-REDUCIBLE 2468 

PLATES REDUCIBLE 5642 

PLATES REQUESTED FOR REDUCTION 3571 

PLATES VALIDATED 2643 

PLATES REJECTED 35 

PLATES IN PROCESSING 222 

PLATES TO BE REDUCED 671 

PERCENT OF TASK COMPLETED 74. 9% 


SPEOPTS NETWORK CLOSED 


JAN 1969 



GRARR DATA VALIDATION STATUS 


NUMBER OF PASSES TAKEN 

NUMBER OF PASSES VALIDATED 

PERCENT DATA PASSING VALIDATION CRITERIA 


2024 

1140 

94% 


TIME PERIOD VALIDATED 


FEB 19 6 8-APRIL 19 69 



STATION OPERATIONAL DATES 


TOTAL TIME PERIOD 

• MADGAR 

• ROSMAN 

9 CARVON 

• ALASKA 


FEB 1968 - SEPT 1969 

FEB 1968 - MAY 1968 
SEPT 1968 - JAN 1969 

FEB 1968 - SEPT 1968 

FEB 1968 - JUNE 1969 

FEB 1968 - SEPT 1969 


Figure 3 


65/66 



IONOSPHERIC CORRECTIONS 


BASED ON VHP RANGE 
AND RANGE RATE MEASUREMENTS 

P . E. Schmid 

Goddard Space Flight Center 

S. Rangaswamy 
Nat* 1 Acad of Science 

C. W. Murray , Jr . 
Goddard Space Flight Center 


Prepared for 

GEOS-2 Review Conference 


June 22-24, 1970 


ABSTRACT 


The interpretation of radio tracking data will lead to bias 
in range and range rate computed from these data if one assumes 
that the radio signals travel with freespace velocity - that is, 
the speed of light in a vacuum. This propagation bias due to the 
ionospheric plasma is directly proportional to the number of 
electrons along the ray path linking the vehicle and tracking 
station. Consideration of appropriate equations for GEOS-2 
ionospheric corrections led to the conclusion that the integrated 
electron content can be extracted directly from VHF tracking 
data., This paper outlines the theoretical aspect of this unique 
ionospheric correction scheme and presents experimental verifica- 
tion based on analysis of Goddard Range and Range Rate VHF 
tracking data. 
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INTRODUCTION 


A careful review of the biasing effect of the Earth's 
ionosphere on G EOS-2 tracking data led to the conclusion that the 
integrated electron content required for ionospheric corrections 
is inherent, at least at VHF (nominal 150 MHz), in range, range 
rate, and angle tracking data. By refining the analysis pro- 
cedure this technique may prove useful during spacecraft tracking 
at frequencies as high as 2GHz even though at 2GHz the ionospheric 
effect is reduced by a factor of approximately 20 0. At this point 
it is worth recalling the reasons why VHF tracking is still widely 
used by NASA and others, even though based solely on radiowave 
propagation considerations, the higher frequencies 1GHz and above 
are much more desirable (reference 1) . The two distinct advantages 
of the lower frequencies such as used by the VHF Goddard Range 
and Range Rate tracking system (uplink 148 . 260 MHz, downlink 
136.140 MHz) are ease of signal acquisition and simplicity in 
spacecraft antenna design . Ease in signal acquisition is especially 
important during near-Earth passes of spin stabilized unmanned space- 
craft in highly elliptical orbits where lunar distance as well as 
near-Earth tracking is desired without antenna electronic or 
mechanical despinning and/or switchover. Such spacecraft are included 
in Goddard ' s Explorer Series and this paper will b e concerned 
with one of these spacecraft, namely. Explorer 34 (also termed 
IMP- 4, Interplanetary Monitoring Platform) . 

Spacecraft antenna simplicity (i.e. , no mechanical or electronic 
beam steering) and ease of acquisition by the ground tracking station 
at VHF both result from the following basic antenna principles: 

1. Antenna conical beamwidth is inversely proportional to 
the square-root of antenna power gain. 

2. For a given effective antenna aperture, as frequency increases 
beamwidth decreases . 


69 










3. Effective aperture is closely linked to physical aperture 
or cross-section. 

4. Energy transfer between an omnidirectional antenna and 

a fixed aperture antenna is independent of frequency providing the 
omnidirectional pattern falls within the beamwidth of the latter 
(reference 2) . 

Thus, the motivation for the analysis presented in this 
paper is two fold, namely: 

1. The experimental verification of GEOS-2 ionospheric 
correction equations as applied to range tone delay and carrier 
Doppler phase, and 

2. Investigation of a possible procedure for real-time 
ionospheric corrections to VHF Goddard Range and Range Rate 
tracking data such as associated with the forthcoming RAE-B 
(Radio Astronomy Explorer) scheduled for launch in late calendar 
1972. 

The basic tracking data handling for such electron content 
extraction is indicated in figure 1. Metric conversion is the 
conversion of time delays, Doppler count times and X and Y 
angles to the desired quantitities range, range rate, aximuth 
angle and elevation angle. The polynomial smoothing and differ- 
ententiation leads to those quantities used in the calculation of 
integrated electron content (elevation angle, elevation angle rate, 
range rate, and differentiated range) . It also provides an 
estimate of the noise on each of these parameters which in turn 
is used to calculate the one-sigma uncertainty associated with 
the calculated integrated electron content at the mid-point of 
the selected data arc . 

The midpoint of the data arc is selected because: 

1. Calculation of the one-sigma uncertainty is easier to 
compute at the center and. 


71 



2 . For uncorrelated noise and odd-ordered polynomials fit 
to the data, the uncertainty is a minimum at the center. For 
uncorrelated noise and even-ordered polynomials, the uncertainty 
is not a minimum at the exact data span center. However for any 
reasonable data stretch the difference between the ratio of the 
integrated electron content one-sigma uncertainty to the value of 
the integrated electron content as calculated at the center and 
as calculated as the exact point of minimum uncertainty is negligible. 
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TRACKING SYSTEM AND GEOMETRY 


Since the experimental verification o£ the GEOS- 2 ionospheric 
correction equations is based on Goddard VHF tracking data from 
the Explorer 34 spacecraft, a brief description of Explorer 34 
orbit parameters and Goddard Range and Range Rate (GRARR) tracking 
system characteristics is in order. 

ORBIT PARAMETERS: 

Explorer 34, a spin stabilized 160 pound spacecraft, was 
launched on May 24, 1967 into the highly elliptical four -day 
period orbit indicated in figure 2 . The vehicle reentered on 
May 3, 1969 . This spacecraft was also known as IMP -4 or IMP-F. 
Currently another vehicle, IMP-G, is collecting scientific data in 
an orbit quite similar to that employed by IMP-F . This affords 
the opportunity of collecting tracking data during passes which 
accentuate the biasing effect of the ionosphere while minimizing 
data noise. As will be discussed more fully later, the best 
results in terms of low variance and repeatibility in extracting 
electron content from VHF GRARR data have been during near-Earth 
(i . e . , perigee) passes. It is expected that use of the 100 kHz 
tone (instead of the usual 20 kHz) and the use of 10 second 
Doppler integration times (instead of 1 second) will lead 
to satisfactory electron content calculations at lunar distances. 

The effect of spacecraft spin on range rate derived from carrier 
Doppler measurements is a most important consideration which is 
discussed in detail in reference 3, 4 , and 5. This will not be 
elaborated upon at this time other than to state that the bias due 
to spacecraft spin has been completely modeled out in this analysis . 
In terms of magnitude , the spin bias for IMP -4 is on the order 
of 3 cm/sec whereas the bias due to the ionosphere during perigee 
passes is on the order of 30 cm/ sec corresponding to night time 
values of integrated content measured as 10-^ electrons/m2 . The 
absolute magnitude for the spin induced bias in a phase modulated 
tracking system is (from reference 4) : 
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CARRIER FREQUENCY = 136.140 MHz 
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where 

k = spacecraft transponder turnaround ratio 

A = up link transmission wavelength 

f = spacecraft spin rate. 

The + sign for a given spin sense depends upon relative 
spacecraft transponder receive-transmit polarizations. Note that 
the result is independent of the aspect angle. 

The absolute sign of the spin bias is positive if the transmit 
polarization is in the same sense as the spin (IMP-4), negative 
if opposite sense. For the VHF GRARR system k = 12/13, A = 2 m. 

For IMP-4 f s = 0.4 revolutions per second and the minus spin in 
equation (1) is appropriate. 

The VHF GRARR tracking sites are listed in figure 3 . Each 
site also has S-Band GRARR capability such as used to track GEOS-2 . 

At present simultaneous VHF and S-band tracking is not possible 
because of shared subsystems . 

VHF TRACKING SYSTEM PERFORMANCE : 

The GRARR tracking system, like many similar schemes, employs 
range tones narrow band phase modulated (modulation index less than 
unity) onto a carrier . This signal is transmitted from ground station 
to spacecraft where it is translated by an onboard oscillator and 
re- transmit ted to the ground site. The downlink carrier is derived 
from the onboard translation oscillator . Sidebands proportional 
to the Doppler shifted uplink frequency are phase modulated onto 
this carrier. After ground receiver combining and detection of 
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carrier and sidebands a Doppler output is obtained which corres- 
ponds to the result which except for spacecraft spin bias, would 
also have been obtained had a coherent (i.e., phase locked) 
transponder of turn-around ratio of k = 1 been employed . 

The key to the ionospheric data derived from raw VHF 
tracking data is the different manner in which the electron plasma 
affects the range tone (group delay) and carrier Doppler range 
rate (phase delay) . This will be discussed more fully in the next 
section . However, the important point, in terms of required 
tracking system performance is low random noise for angle, angle 
rate, range rate, and range measurements . The systematic range 
bias , as will be shown, is of no consequence since it is differ- 
entiated out. Also there is no reason to expect systematic range 
rate offsets (other than due to spacecraft spin and ionosphere) since 
the Doppler comparisons are made in real time with two-way 
radiowave propagation time short relative to local and reference 
oscillator short term stabilities (reference 6 ) . 

Figure 4 shows the expected range noise related to slant 
range from the tracking site . The calculation is based upon an 
assumed omnidirectional (OdB) spacecraft antenna, 1 watt radiated 
power and no radio frequency system losses. IMP -4 radiated 4 
watts (reference 7) via a turnstile antenna (nominal OdB gain rela- 
tive to isotropic) and the measured range noise levels during peri- 
gee passes range typically from 3 to 10 meters . This is considered 
reasonable agreement since the radio frequency losses plus variations 
in antenna gain versus aspect angle can easily add 5 to 10 dB 
fluctuations . Three important points are to be noted regarding 
GRARR VHF ranging - first, the range tone j itter becomes thermal 
noise limited at approximately 20,000 km for the assumed parameters - 
second, the entire noise level can be reduced by a factor of 5 by 
switching in the 100kHz range tone and finally the thermal limited 
noise can be reduced by a factor of 3 by switching in the 0 . 1 Hz 
(instead of 1 . 0 Hz) tracking filter. The range measurement noise 
is linked to range tone phase measurement noise by: 
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iR = (i) (r) 

where 

A$ = phase jitter (radians) 

= *- one wavelengths (reference 7) 


Using analogous link parameters one obtains the range rate 
noise characteristics indicated in figure 5. Usual near-Earth 
operation (i.e, , unless otherwise requested) corresponds to 1 
per second data at a 1 kHz loop bandwidth. Measured noise 
in this case , however, exceeds the 6 cm/ sec indicated in figure 5, 
due to the cyclic component to the Doppler measurement introduced 
by the non-uniform phase pattern of the rotating spacecraft antenna . 
If this harmonic component (typically 30 cm/ sec amplitude sinusoid) 
is modeled out range rate noise values corresponding to figure 
5 are achieved . However, at present only the spin induced Doppler 
bias is modeled out while the sinusoidal variation is taken out 
by means of the polynomial smoothing (figure 1) . 

Finally, the integrated electron content measurement while 
not a critical function of elevation angle (figure 6) , is quite 
sensitive to angle rates . Thus , boresight pointing biases are 
again differentiated outj however, it appears that monopulse angle 
data taken below 20° elevation begins to become corrupted in 
harmonic fashion by multipath propagation as a result of the 
relatively broad ground site antenna beamwidth (nominally 15°). 

In the data used for this paper, one sigma noise over a one to 
three minute data time span for angle data ranges from 0.01° 
to 0.07° 
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EXTRACTION OF INTEGRATED ELECTRON CONTENT FROM TRACKING DATA 


The following will show how the comparison of differentiated 
range based on sidetone ranging with range rate based on carrier 
Doppler leads to a value of the integrated electron content, I , 
in a vertical column above the tracking site. 


RANGE IONOSPHERIC BIAS: 


Figure 6 indicates the need for range corrections at VHF 
frequencies. Dyatime range biases easily exceed a kilometer. 

Range rate bias, which results from the changing phase path length 
as the ionosphere is scanned, can introduce meters per second of 
systematic error compared to system resolutions on the order of 
0.5 cm/sec (refer to figure 5). Because of the extreme variability 
of the daytime ionosphere (e.g. , factor of 3 day-to-day variations 
of total electron content) modeling of such biases has not been 
very successful. If VHF ranging is employed, real time ionospheric 
bias corrections are highly desirable. The magnitude of the range is 
bias can be estimated by: 


AR 



40.3 



40.3 


kf sin E 


(3) 


(references 1 and 8) 


where N = refractivity of ionosphere along path between ground 
station and spacecraft (index of refraction is given 
by n » 1 + n(10 n< 1.0 

f = frequency of transmission, for VHF GRARR use 
geometric mean i.e. 1/2 
(Reference 9) 


1 + 1 


1/2 


if ^ f 2 

' u £ d ) 


in Hz 
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3 

N g = electron density along ray path (electrons/m ) 

I = vertical integrated content above tracking station 

v 2 
(electrons/m ) 

E = local elevation angle at the ray path at an altitude of 
400 km = elevation angle at ground station 


Equation (3) expresses the measured excess in range over 
free space range when the speed of light is used in rotating 
range tone group delay to true range assuming tropospheric bias 
is modeled out. Tropospheric component is frequency independent 
in the usual range of interest (20 MHz to 20,000 MHz) and has 
a maximum of approximately 25 meters at 5° elevation (reference 1) . 
The troposphere range bias also has a 1/sin E dependence. At 
140 MHz the ionospheric bias is seen (figure 6) to be in general 
much greater than that attributable to the troposphere. One 
can differentiate equation (3) with respect to time to obtain: 


AR = 


AR. - AR^ 
1 2 

AT 


/40.3 cos E 
V sin^ E I 



(4) 


where: 

AR = ionospheric range rate bias based on range tone measure- 
ments . 

E = elevation angle rate (radians/sec) 
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and all other parameters are as previously defined. 

A more sophisticated model for the ionosphere (see equation 3) 
might be inserted for the integrated content along the ray path 
to obtain: 

= E jk /tysias (5) 


However, equation (4) is what was used in the preliminary 
experimental verifications presented in this paper. 

RANGE RATE IONOSPHERIC BIAS: 


The ionospheric bias on the range rate measurement based carrier 
Doppler, as will be shown, is approximated by: 


Ar = 


r 40.3 cos E 
2 

^ sin E 



* 

E 


( 6 ) 


where Ar is the ionospheric range rate bias based on carrier 
Doppler measurements 

It will be noted that the only difference between equations 
(5) and (6) is in sign. This is the basis of the integrated 
electron content calculation since (5) and (6) can be combined to 
yield: 
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sm E 


r 


dR 

4 - 


(7) 



dR = differentiated range measurement 
dt 


THEORETICAL BASIS FOR INTEGRATED ELECTRON CONTENT CALCULATION: 

The basic principle involved in the integrated electron 
content extraction can be derived by considering the somewhat 
idealized case of unity frequency turnaround at the spacecraft. 

The resultant ionospheric range tone and Doppler phase bias thus 
determined is identical to that obtained with th e actual GRARR 
tracking system except the frequency to be considered in the 
ionospheric bias calculation is a function of both up and down 
link frequencies (refer to equation 3) rather than a single 
frequency (reference 9) . 

Consider the following narrowband phase modulated signal: 

E = cos ( o) Q t + 3 sin oj t ) (8) 

where: 

= 2Trf Q angular uplink carrier frequency 
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3 = index of modulation ^ 1.0 


w m =2Trf m angular range tone frequency 

As is well known, the carrier and principal sidebands of 
such a signal are of the form: 

cos to t 
o 

cos (w + oo ) t (9) 

o m 

-COS ( 03 “ 03 ) t 

o m 


which is the type spectrum transmitted up to spacecraft and 
back during GRARR tracking. The highest frequency tone, f^, 
normally used with VHF GRARR is 20 kHz. Since this is a phase 
modulated scheme the absolute signal amplitudes during normal 
operation do not affect the range or Doppler measurement. 

Thus, the effect of two-way transit through the ionosphere on 
these three discrete frequencies is to be examined in terms of 
the Doppler and range phase demodulation at the ground site. 

Doppler Phase Measurement 

The effect of two-way transit through the ionosphere on 
Doppler phase measurement is as follows: 
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cos w It - - 


but n = 



(see equation 3) 

Therefore the range rate phase is given by: 
_ -2^ f 40.3I V 




R - 


sin E 



( 10 ) 


where : 

R = "true" slant range based on 

time delay 
C = speed of light 
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Range Tone Phase Measurement 


The effect of the same two-way transit on the range phase 
delay measurement is as follows: 

(“o + “ m ) (* ■ t/ n l ds )] f[( OOS “o~“J ■ i/ n 2 ds ). 

i SIDEBANDS 


cos 

CARRIER 


COS to 


o (t - i 



The foregoing results from the frequency dispersion introducted by 
the ionosphere which by way of equation 3 is linked to the respective 
index of refraction calculations in the following manner: 


c / nds 

C J 1 

-fn 0 ds 
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Differentiation of equation 10 results in the usual expression 
for range rate and in the GRARR system is approximated by averaging 
the phase change over the time required to count a fixed number 
of cycles of the difference frequency plus an inserted offset 
bias frequency which is required to permit range rate sign deter- 
mination. The difference between "time" range rate at the center 
of the count interval and the measured "average range rate" is 
a function of averaging time and the magnitudes of the odd time 
derivatives (i.e., 3, 5, 7...) of range (reference 10). For the 
IMP-4 data arcs described in the next section, the difference 
between "true" and "average range rate" was negligible. Equation 
11 presents the measured range tone phase. Differentiation of 
this result and comparison with carrier Doppler range rate results 
in equation (7) presented previously for calculation of integrated 
electron content I v . If the spacecraft is spin stabilized, then 
as indicated by equation (1) , the carrier Doppler determined 
range rate must be adjusted accordingly. 


EXPERIMENTAL RESULTS 

In order to facilitate data analysis the plotting routine 
indicated in figure 7 was written. Each data stretch is examined 
to provide the necessary inputs required for the calculations of 
integrated electron content, namely, elevation angle, elevation 
angle rate, measured range rate and differentiated range. A 

a 

typical film printout of the comparison of r and dR is given 

by figure 8. This film sequence which presents tft£ calculated 

electron content also includes such pertinent parameters as the 

one-sigma limits on range noise, range rate noise and integrated 

electron content. Station identification, epoch, and tracked 

spacecraft nomenclature are also part of the standard output. 

17 ) 

Figure 8 resulted in a calculated value of I v of 1.08 (10 
electons/m^ in a vertical column over Santiago, Chile at 10 PM 
local time on December 20, 1968. Subsequent values of I v cal- 
culated at Santiago over a period of 90 minutes (figure 9) show 
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good agreement with an observed RMS deviation from the mean of 
10%. This is also the calculated uncertainty based on angle 
range and range rate measurement noise. Figure 9 also includes 
two calculations from Carnarvon, Australia on August 29, 1968 
at approximately 11 PM local time. The possibilities for further 
reducing measurement noise are discussed in the conclusions 
presented as the final section of this paper. 

Figure 10, which is based on independent Faraday rotation 
data (reference 11) is presented as tentative verification of 
the electron content values extracted at, for example, Carnarvon 
on August 29, 1968, and Santiago on December 20, 1968. Since 

”| r l O o 

the total content varies between 10 to 10 electrons/rn 

depending on time of day, season, and state of solar activity, 

tlie calculated night time values of 10 electron/m 2 appear quite 

reasonable for 1968, which was near the peak of the 11 year 

sunspot cycle. Figure 11 shows the night time (10 PM local) 

1 7 

integrated content over Hawaii to be on the order 0.4 (10 ) 

2 

electron/m which reflects the lower sunspot activity of 1964 . 

The last solar sunspot minimum occurred in 1964 and the maximum 
in 1968. 

CONCLUSIONS 

At VHF range measurements uncorrected for the ionosphere can 
be in error by kilometers. Uncorrected range rate can be biased 
by the ionosphere on the order of meters per second. Since the 
VHF GRARR tracking system resolution is on the order of a few meters 
in range and with a 10 second integration less than one cm/sec 
in range rate these biases must be modeled out. Since the range 
and range rate ionsopheric induced bias is directly proportional 
to I v the integrated electron content, it is advantageous 
to measure I v in real time. 

The often used integrated electron content values based on 
average monthly predictions can easily be in error by 50%, 
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especially during the daylight hours. This paper has presented 
the theory and preliminary experimental verification whereby the 
desired quantity, I v , is extracted directly from tracking measure- 
ments of angle, range and range rate at a nominal 140 MHz . 

The accuracy of this electron content determination scheme 
can be improved and the concept perhaps extended to reliable 
lunar distance real time corrections . Such improvements are 
anticipated as a result of considerations to: 

1. Include earth curvature 

2 . Model out harmonic terms in data due to spacecraft spin 
and multipath 

3 . Use longer doppler averaging time ( 6 per minute setting) 

4 . Use higher frequency range tone (100 kHz setting) 

5 . Use elevation angle rates based on orbit computation 

6. Record VHF tracking data during more daytime passes 

(IMP-6) 

7. Intercompare with best estimate ray traces based on 
top and bottom side sounding 

8. Determine optimum data pass length 

9 . Investigate the calculation of electron content at 
points other than the midpoint of the data arc . 

This work is part of a continuing effort within the Mission 
and Trajectory Analysis Division (GSFC ) to improve the efficiency 
and accuracy of NASA tracking data utilization in orbit and trajectory 

computation . 
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A TWO-LASER COLLOCATION EXPERIMENT 


M. R. Pearlman, C. G. Lehr, G. M. Mendes, and M. R. Wolf 

INTRODUCTION 

A two-laser collocation experiment has been performed on the Geos 2 
satellite with the Smithsonian Astrophysical Observatory (SAO) prototype 
ranging system and the NASA mobile laser unit. The object was to deter- 
mine the relative accuracy of these two laser systems, both of which are 
being used in the routine collection of satellite geodetic data. The effort 
was conducted at the SAO Observatory at Mt. Hopkins, Arizona, during the 
period October 1969 to January 1970. The two laser systems demonstrated 
a relative accuracy of 1 to 2 m, determined from short-arc single-pass 
calculations. From a point-by-point analysis of near-simultaneous satellite 
returns, we were able to extract the displacement of the collocated lasers 
to an accuracy of a few meters over a baseline of ~ 130 m. 

The laser data have also been included in 2-week long-arc calculations to 
which Baker-Nunn camera data were added. Laser range residuals of 10 m 
or better and angular-data residuals of 3 to 6 arcsec were achieved. From 
this long-arc analysis, relative station positions were determined to better 
than 10 m, 


THE LASER SYSTEMS 

The SAO ranging system is the prototype of the units being deployed at 
the Observatory^ network of stations. It uses a ruby laser with an oscillator- 
amplifier configuration. The laser has an output of 7 . 5 J and a nominal pulse 

This work was supported in part by grant NGR 09 -015-002 from the National 
Aeronautics and Space Administration. 
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width of 1 5 nsec. The system has a static pointing mount that permits the 
laser to be pointed to predicted positions and fired by the clock at the appro- 
priate epoch. This mode of operation enables us to range to satellites when 
they are in the earth's shadow and during daylight hours. Range measure- 
ments are made with a time- interval counter with a 1 -nsec resolution. 

During the experiment, the laser was operated at a firing rate of 1 to 2 per 
min. Further details on the system are given by Lehr (1969) and Lehr and 
Pearlman (1969)* 

The NASA mobile unit (MOBLAS) used a ruby laser with an output of 
1 J and a pulse width of 30 nsec. The laser has a firing-rate capability of 
1 per sec. For the experiment the system had an on-site computer facility 
for tracking and searching purposes and used a 10 -nsec time -interval counter 
for range measurements. The mobile laser is discussed in detail by Johnson, 
Plotkin, and Spadin (1967). 

The NASA station was located approximately 100 m north of the SAO laser 
(see Figure 1). The experiment was divided into two phases, corresponding 
to two NASA locations (HOMLAS and HOMLA2) approximately 10 m apart on 
an east-west baseline. Phase I lasted through November 15, while Phase II 
comprised the remainder of the experiment. For the experiment, the SAO 
Mt. Hopkins laser was referred to as HOPLAS. Our intention was to retrieve 
the IGG-rn interstation separation and then to recover the 1 0 -m shift in the 
MOBLAS station location. 


SHORT-ARC ANALYSIS 

We used short-arc single-pass analyses to make a simple and direct 
intercomparison between the two laser systems (see Tables la and b). Owing 
to the relatively large data rate of the NASA laser, we used their data to com- 
pute a short reference arc with which to compare the SAO satellite returns. 
The relative bias between the two systems was taken as the difference between 
the system biases to the reference arc. The standard deviation was taken 
as the square root of the sample variance. 
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During both phases of the experiment, the two systems demonstrated a 
relative ranging accuracy of 1 to 2 m. The relatively small number of SAO 
points listed per pass is attributed mainly to the relative data rates, but it 
also reflects the method of analysis. Only the SAO satellite returns that 
were encompassed by the NASA data or that fell very close to the extremities 
of the NASA data have been included. 

The sign of the bias between the systems changed several times during 
the 4-month experiment, but it was predominently negative during Phase I 
and positive during Phase II. This bias component was probably introduced 
into one or both of the systems during the calibration procedure, which 
involved a determination of the system delay by ranging on a target at a known 
distance from each laser. 

Measurements of this nature are extremely sensitive to changes In return- 
pulse height caused through variations in output-power level, receiver gain 
(including photomultiplier settings), counter trigger threshold, or attenuation 
in the optical train that may have occurred during the calibration procedure. 
Even a change of a factor of 2 in the ratio of pulse peak amplitude to trigger 
level can introduce a change in the system calibration equal to 25% or more 
of the pulse halfwidth. In our experiment, this error could be 0,5 m for the 
SAO system and 1.0 m for MOBLAS. These errors would be considerably 
larger for larger changes in the ratio of peak-signal strength to trigger 
threshold. 

The standard deviations in the data for each system indicate that in the 
present configuration, each has a noise level of the order of 1 m. Consider- 
ing the size of the output-pulse widths and the pulse-height statistics, this 
value appears to be quite reasonable. These systems would be capable of an 
improved signal-to-noise value if shorter pulses were used or if calibration 
and satellite range measurements took some account of pulse size and shape 
(Lehr, Pearlman, and Scott, 1970). 
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Several system errors were detected by the early collocation results. 
These included errors arising from system misalignment, calibration taken 
infrequently or not taken immediately following an instrumentation modifica- 
tion, and even overshooting the target during calibration. Each of these 
involved an error amounting to no more than a few meters and it is indeed 
gratifying that the lasers were capable of detecting these offsets while ranging 
to satellites. 

DETERMINATION OF RELATIVE STATION POSITIONS 

We calculated relative station coordinates by examining the pertinent 
geometry for pairs of nearly simultaneous data points taken by the two laser 
systems (see Table 2). Both the interstation baseline and the 1 0 - rn NASA 
station shift are clearly discernible in the data. The values for station 
separation are solutions of a least -squares fit of the appropriate data, and 
the residuals are the standard deviations. The SAO firing epoch was adjusted 
to be as close as possible to that of the NASA laser. Epoch differences were 
compensated for by linear interpolation between successive NASA data points. 
With this interpolation procedure, separations in firing time of up to 10 msec 
appear to be consistent with a 1 -m accuracy. The 10-msec criterion was 
used in selecting data. Although the interpolations involved approximations 
and there were few data points, the results bear out the 1- to 2-m accuracy 
capability of the present laser systems. 

LONG -ARC ANALYSIS 

Long-arc calculations permit us to add significant amounts of Baker- 
Nunn data to the analysis and thereby test the relative accuracies of two 
independent techniques. In addition, long-term orbits are essential for the 
dynamical geodesy that we hope to do with the satellite data. Accordingly, 
it is of interest to see how well we can fit long arcs to the data. 

We used the 1969 Smithsonian Standard Earth (II) to calculate two-week 
orbits from laser and photoreduced Baker -Nunn data (see Table 3). To maintain 
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a reasonable balance in the data, we selected every tenth NASA return 
instead of using the total data. In the orbital solutions, the laser and camera 
observations have been given weights of 5 m and 4 arcsec, respectively. For 
these solutions, the standard errors of unit weight are in the vicinity of 1.0, 
indicating that the choice of weights is reasonable. It must be pointed out 
that we have a preponderance of laser data that are not well distributed; 
consequently, our orbits probably do not represent the best that could be 
obtained. 

The value chosen for the Baker-Nunn camera data is an estimate of its 
accuracy. The residuals obtained in the solution were in the neighborhood 
of 3 to 6 arcsec. 

The accuracy of the laser data appears from the short-arc analysis to 
be well below the 5 -m weight ascribed to it. However, the uncertainties in 
the gravity-field model, which are larger than the errors in satellite range 
measurements, dominate the solutions for the laser data. 

When the standard error of unit weight is very nearly equal to 1.0, the 
laser residuals are close to the 5-m weighting factor. For values below 
1.0, the residuals are small, indicating that the lasers may be weighted too 
heavily for this particular solution. Alternatively, for data in which the 
standard error of unit weight is above 1;G, the lasers are probably not being 
weighted enough, and indeed the accompanying range residuals are quite 
larger 

Relative station locations have been extracted from the long-arc solutions 
(see Table 4). Averaging over the limited amount of data does not permit 
us to overcome the existing modeling errors. We are able to separate the 
NASA and SAO station locations hut do not have enough accuracy to measure 
the MOBLAS station shift with any reasonable certainty. 
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CONCLUSIONS 


We conclude that lasers are capable of ranging to an accuracy of 1 to 2 m 
with respect to one another. We also find that the laser data, when compared 
to the Baker -Nunn camera data, appear to be at least as good as the gravity 
model used in the calculation. 

Considerably more data would be required to obtain better results with 
long-arc analyses. This conclusion comes as no surprise, as we know the 
enormous amounts of data that will be required to pursue any reasonable 
program in dynamical geodesy. 
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LONG-ARC RESIDUALS 
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INTERSTATION SEPARATION 
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MINITRACK SELF CALIBRATION 


Goddard’s Minitrack interferometers provide the bulk of satellite 
tracking for Goddard’s Space Tracking and Data Acquisition Network (STADAN). 
They are still calibrated by comparing Minitrack and optical observations 
of a number of aircraft passes. There are several possible sources of error 
in the Minitrack system but the principle one is a phase bias added by the 
interferometer to the phase difference that is being measured. Fortunately 
this phase bias , or zero set error* 1© fairly stable and may remain constant 
for many years. However, when it changes it can be a sudden change and all 
observations taken with that particular interferometer baseline will carry 
that bias until the next aircraft calibration -- a time period that may be 
as long as a year. 

It would clearly be desirable to devise a calibration scheme for 
Minitrack that would update calibration values more often than can be done 
with the aircraft. Several schemes have been tested but the one that offers 
the best hope is a bootstrap method that regresses on zero set values 
simultaneously with the determination of long arc orbits. Minitrack 
interferometers are located at eight STADAN sites. There are four primary 
baselines at each site, each having its own zero set value. A satellite 
in an orbit of sufficiently high inclination can, in time, be observed by 
all 32 baselines. Consequently in addition to the orbital elements there 
are 32 bias parameters to solve for. These additional unknowns require 
that the observational period be of sufficient length to allow for adequate 
statistical control. An accuracy in determining the orbit of about eight 
seconds of arc is sufficient for the determination of Minitrack zero set. 
Even less accuracy is tolerable if the errors in the computation have a low 
correlation between one pass at a given station and the next. 

The NONAME orbit generator that has been used extensively at Goddard 
for geodetic studies, taken with recent earth model refinements based on 
GEOS II observations , appeared to be a suitable instrument for testing 
the feasibility of the self calibration concept. NONAME has demonstrated 
an accuracy in fitting optical observations for periods up to a week that 
is more than sufficient for the purpose of Minitrack calibration. Further- 
more, one version of NONAME had the convenient capability of regressing on 
Minitrack zero sets simultaneously with the determination of orbital 
parameters . 

Sufficient data to perform this test was available from the GEOS II 
data bank. Minitrack and MOTS camera observations were chosen covering the 
period in 1968 from February 21 to March 27 . 

The basic strategy of the experiment was to take a week’s worth of 
Minitrack observations of GEOS II and, using the NONAME program, solve 
for Minitrack zero set. These zero sets were then to be applied to the 
preprocessing of the Minitrack observations of GEOS II made the following 
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week. Using NONAME again* an orbit would be fitted to the second week's 
batch of observations, this time solving for orbital parameters alone. 

This latter orbit would be compared with the MOTS camera observations taken 
during that time. Then another orbit would be computed for the second 
week, only in this case the Minitrack observations would have been pre- 
processed with zero sets taken from aircraft calibrations of the inter- 
ferometers, It, too, would be compared with the MOTS observations, (The 
MOTS observations on GEOS II, accurate to two seconds of arc, or better, 
are an excellent standard of comparison,) Th© orbit that best agreed with 
the MOTS sightings would tell which calibration method is best. (In 
actual practice this confidence would be enhanced by solving for zero 
sets using the observations of a number of different satellites in different 
orbits. This would nullify any interaction between the determination of 
zero sets and the gravity model errors that might be peculiar to a given 
orbit (resonant terms) and it would ensure good coverage for each set of 
baselines in the tracking network.) 

Actually, not just two but five successive weeks of data were studied. 

Not only was the basic experiment carried out but a number of variations 
and side paths were explored. 

In the determination of Minitrack zero set from orbital data no a 
priori information, i.e., previously determined aircraft calibration values , 
was employed. Initial estimates of phase bias were arbitrarily set at zero 
when, in fact, they might be up to +180° for a given baseline. This device 
of starting in a degenerate condition was employed to allay any doubt that 
a self calibration system might slowly diverge from a previously established, 
satisfactory set of calibration values. If one could converge to a proper 
set of estimates from a starting point of total ignorance, then the doubt 
would be unwarranted. 

To date, the results of this experiment have been negative, i„e., 
other than expected. The solutions for Minitrack zero sets converged to a 
set of incorrect values, and the simultaneously determined orbital parameters 
exhibited corresponding errors. These errors were large and systematic. 

The errors in zero set determination ranged as high as 4 to 5 times the errors 
that typically exist when the aircraft calibration values are used. The 
error in determining the orbit was principally along track. The error was 
primarily manifested as a declination error when the computed orbit was 
compared to MOTS camera observations. The bias in these declination residuals 
ranged up to a minute of arc. 

Table 1 shows the MOTS camera residuals when these observations are 
compared to GEOS II orbits whose parameters have been determined simultan- 
eously with the determination of Minitrack zero set. The optical data was 
given a weight of zero in this computation. Chiefly to be noted i® the 
consistent negative sign and the relatively large value of the average 
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errors in declination. The values in the table are an average of all the 
points from all the passes taken by the given camera during the week ? s 
period for which the orbit was computed. Each of the five trajectories 
was based on a different set of Minitrack observations and there was no 
overlap in either the Minitrack data or in the optical data from one week 
to the next. The bias in the declination residuals for weeks three and four 
are less than, those for the other weeks but the bias is in the same direction. 
It should be pointed out that the relative positions of sun, earth* and 
satellite orbit precluded any photographing of GEOS II by MOTS during this 
p&r iod when the ttaeolllta wins moving north to aioueh. 

Table 2 shows the difference between the Minitrack calibration values 
computed by the simultaneous regression on orbital elements and zero sets 
and those determined by aircraft calibration. The numbers given are averages 
of the values calculated separately for each of the five weeks. The two 
methods of calibration are in much better agreement for the east-west base™ 
lines than they are for the north-south baselines. This is contrary to what 
had been expected. The expectation was that* for a satellite in polar orbit 
like GEDS II, any adverse interaction between the computation of along track 
satellite position and baseline zero set would be self correcting. The 
north-south baseline intercepts north going passes on one side of the orbit 
and, twelve hours later, south going passes on the other. A balance of passes 
in the two directions was presumed to negate the possibility of compensating 
biases developing in the regression on zero set and orbital parameters. 
Apparently the interaction between these factors is more complex than the 
intuitive model allowed for. 

Figure 1 shows three curves that represent the along track, cross track 
and radial components of the differences between the reference orbit determined 
by MOTS camera observations and the orbit resulting from simultaneous 
regression on zero set values and orbital elements with Minitrack observations 
as sole input. The graph shows only about one orbital cycle but this same 
pattern is repeated throughout the seven day period covered by the orbital 
fit . 


The results of this experiment are still under investigation. Many 
computer runs have been made to check various hypotheses about the basic 
results. The orbits were recomputed with various combinations of Minitrack 
stations totally omitted from the computation. This produced no major change 
in the results. Numerous checks, some involving simulated data, have been 
made to ensure that no sources of large errors exist in the data processing 
or orbit computation systems. No gross errors have been found but the 
possibility of some form of mishandling of the Minitrack observations has 
not been dismissed. Just beginning is a study of the more subtle statistical 
elements in the distribution and character of the Minitrack observations. 
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Whether Minitrack self calibration is feasible is still a question. 
Certainly the attempt described here produced an unsatisfactory result. 

Why this was so is not understood. Gaining this understanding is the next 
step. Such an understanding should help in determining the proper technique 
for Minitrack self calibration if, indeed, this is a practical possibility. 
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Table 1 


OPTICAL RESIDUALS (ARC SEC) FOR 
MINITRACK ORBIT & BIASES SOLUTIONS 


STATION 


Hum 

Week #2 


Week #4 

Week #5 

wmrm 

Rt . Aac. 

-10,3 

25.5 


1 9 4 

0*0 

Dec 1 in 

-56.9 

-46.5 


7.8 

-49.1 

SNTAGO 

Rt „ As c * 
Dec 1 in 

11.9 

-74.8 



-7.4 

-6.3 

8.3 

-33.3 

MOJAVE 

Rt 6 As c . 

-1.6 

21.2 

14.0 

2.4 

5.3 

Dec 1 in 

-45.8 

-35.1 

-18.1 

3.2 

-57.6 

JOBURG 

Rt * As c o 
Be cl in 

15.8 

-66.3 



-1.8 

18.9 

-56.7 


Rt. Asc. 

-9.8 



11.5 

-1.5 





Table 1 (Continued) 


STATION 


Week #1 

Week #2 

Week #3 

Week #4 

Week : 

JAMAC 

jRtl O Ail 6 Q 

-7.3 

12.5 


-4.6 

-1.3 

Dec 1 in 

-26.9 

-20.8 


4.8 

-38.7 

GSFCN 

Rt. Aac. 

-5.0 

17.2 




Dec 1 in 

-55.1 

-32.2 




ROSMA 

Rt. Asc. 

-8.6 

22.3 

6.0 

-7.3 

0.3 

Declin 

-48.0 

-36.0 

-22.5 

2.0 

-40.2 

ULASKA 

Rt. Asc. 


-1.4 

-12.9 

-5.2 

-0.1 

Declin 


-15.9 

-5.2 

— 5 » 4 

-35.8 

TANANA 

Rt. Asc. 
Declin 



23.2 

-80.2 
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Table 2 


FIVE WEEK AVERAGE OF DIFFERENCE BETWEEN 
MINITRACK SELF CALIBRATION DETERMINED ZERO SET 
AND BEST ESTIMATE ZERO SET FROM AIRCRAFT 
CALIBRATION AND QUALITY ASSURANCE 


Station 

Antenna System 

FTMY1S 

Polar 

Equatorial 

QUITOE 

Polar 

Equatorial 

LIMAPU 

Polar 

Equatorial 

S NT AGO 

Polar 

Equatorial 

NEWFLD 

Polar 

Equatorial 

WNKFLD 

Polar 

Equatorial 

JOBURG 

Polar 


Equatorial 

ULASKA 

Polar 

Equatorial 

MADGAR 

Polar 

Equatorial 

ORORAL 

Polar 

Equatorial 


Difference 

Baaellne Wavelengths x 1Q~ 


EW 

2.0 

NS 

-19.6 

EW 

4.0 

NS 

-12.1 

EW 

5.9 

NS 

-4.7 

EW 

9.5 

NS 

4.0 

EW 

11.2 

NS 

-13.5 

EW 

8.3 

NS 

-3.0 

EW 

2.4 

NS 

-0.3 

EW 

-3.0 

NS 

-4.4 

EW 

-0.4 

NS 

-20.5 

EW 

-6.5 

NS 

-13.6 

EW 

-0.6 

NS 

-8.2 

EW 

-2.6 

NS 

-4.8 

EW 

-1.8 

NS 

-8.7 

EW 

0.9 

NS 

-5.9 

EW 

-2.0 

NS 

-6.9 

EW 

-1.6 

NS 

-7.2 

EW 

-3.4 

NS 

-5.0 

EW 

-7.3 

NS 

-5.5 

EW 

-3.0 

NS 

-15.3 

EW 

-1.0 

NS 

-3.6 
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AN INTERCOMPARISON OF NAVY TRANET DOPPLER DATA 
AND OPTICAL DATA FROM THE GEOS-I SATELLITE 

Russell W. Agreen 
James G. Marsh 

ABSTRACT 

Orbital solutions for the GEOS-I satellite obtained from U. S. Navy 
TRANET Doppler data and those from optical flash data recorded by the 
NASA Space Tracking and Data Acquisition Network and Smithsonian 
Astrophysical Observatory systems were intercompared. The orbital 
arcs used for this study were two days in length and were in the period 
July 9 through August 7, 1966. RMS of fits for the orbital solutions were 
on the order of 1.9 seconds of arc for the optical solutions and 2.7 cm/ 
sec for the Doppler solutions. Comparisons of the corresponding opti- 
cal and Doppler orbital ephemerides showed total RMS position differ- 
ences ranging from 20 to 40 meters. Biases in the base frequency value 
of the Doppler data were solved for; consistent biases in the base fre- 
quency around 9 cm/sec ±3 cm/sec were found. 


in 
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AN INTERCOMPARISON OF NAVY TRANET DOPPLER DATA 
AND OPTICAL DATA FROM THE GEOS-I SATELLITE 

by 

Russell W. Agreen and James G. Marsh 
Goddard Space Flight Center 


I. INTRODUCTION 

This report presents the results of an intercomparison between orbits determined from Navy 
TRANET Doppler data and those determined from optical flash data from the GEOS-I satellite. 
This investigation was conducted for the purpose of establishing our capability to handle the 
Doppler data for geodetic purposes. In the past, our geodetic studies have used optical data as a 
standard. 

The NONAME Orbit and Geodetic Parameter Estimation System (Reference 1) was used to de- 
termine all orbital solutions and to generate the intercomparisons. All of the data were obtained from 
the National Space Sciences Data Center at Goddard Space Flight Center. It was taken during the 
periods July 9— July 26 and July 31— August 7 of 1966. 

In addition to evaluating orbits generated from Doppler data, the data were examined for biases 
in the base frequency value and in the time tags associated with each measured frequency (or 
observation). 


II. OBSERVATIONAL DATA AVAILABLE 

Table 1 presents the orbital characteristics 
were from this satellite. 

GEOS-I was launched by the National Aer- 
onautics and Space Admi ni stration under their 
National Geodetic Satellite Program. The orbit 
was chosen to minimize the perturbative effects 
of air drag and solar radiation pressure. On 
board the satellite were flashing lamps, Doppler 
transponders, laser corner reflectors and 


of the GEOS-I satellite; all data used in this study 
Table 1 

Orbital Characteristics of the GEOS-I Satellite. 


Apogee Height 

2273.0 Kilometers j 

Perigee Height 

1116.0 Kilometers 

Eccentricity 

.07 

Inclination 

59.4 Degrees 

Anomalistic Period 

120.3 Minutes 


1 
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other electronic tracking instruments which enabled the numerous participating agencies to re- 
cord large amounts of valuable tracking data. 


Optical Data 

The GEOS-I optical flash data used in this analysis were obtained from the Geodetic Satellite Data 
Services of the National Space Science Data Center located at Goddard Space Flight Center. It is 
composed of observations from the NASA STADAN Minitrack Optical Tracking System (MOTS), 
the Goddard Special Optical Tracking System (SPEOPT), and the Smithsonian Astrophysical Obser- 
vatory (SAO) Baker- Nunn camera stations. All optical station positions are referenced to the SAO 
C-7 system (a e = 6378142 meters, Reference 2). 

The orbits generated from this optical data are used as the standard to which the TKANET 
Doppler orbits are compared because it is generally felt that the optical flash (or active) data from 
GEOS-I forms a very high-precision data set. The accuracy of the optical data is on the order of 
2 seconds of arc which amounts to a positional error of approximately 15 meters for GEOS-I. 

Among the reasons for this confidence are the use of a stable on-board clock to set off the optical 
beacon flash intervals, exacting the time of observation to millisecond accuracy, and the short 
duration (about 1.3 ms.) of the flashes, enabling the cameras to record the observations as point 
images against a background of reference stars rather than less exacting streak images (Refer- 
ence 3). 

Substantial evidence of this quality is seen in the orbital RMS T s of fit; for each of the thirteen 
2-day arcs generated from optical data, the RMS of fit is on the order of 2 seconds of arc. This is 
a very good fit considering that the NASA (MOTS and SPEOPT) and SAO (Baker-Nunn) data were 
processed through independent systems, yet weighted equally in the orbital solutions. The MOTS 
and SPEOPT data were obtained by MOTS 40 and PTH-100 cameras and all plates were reduced at 
the New Mexico State University to yield right ascension and declination pairs and U.T. C. time tags. 
The SAO data were taken by Baker-Nunn cameras on film and reduced by SAO to yield right ascen- 
sion and declination pairs and time tags in A.S. atomic time (Reference 4). 

Some preprocessing of the data received from the Space Science Data Center is done using the 
NONAME orbit determination system. MOTS and, SPEOPT data time tags are corrected for flash 
buildup time; no other preprocessing is necessary. SAO data has atomic time (A.S.) tags and 
NONAME is used to apply the conversion to Universal Time Coordinated (U.T.C.). Also, a transit 
time correction is applied to the SAO data to refer the observation from station time to satellite 
time. Finally, SAO observations are referred to the mean equator and equinox of 1950.0, and it is 
necessary to apply precession and nutation from that epoch to the true equator and equinox of epoch 
of the observations (Reference 1). 

Station network, name, number, location, and camera type are presented in Table 2 for all of 

the optical stations used in this study. 
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Table 2 

Optical and Doppler Stations Used in the Analysis. 


Network 

Station 

Number 

Camera Type 
(optical) 

Location 

ST A DA N 

1BPOIN 

1021 

MOTS 40" 

Blossom Point, Md. 


1FTMYR 

1022 

MOTS 40" 

Fort Myers, Fla. 


lOOMER 

1024 

MOTS 40" 

Woomera, Australia 


1MOJAV 

1030 

MOTS 40" 

Mojave, Calif. 


1JOBUR 

1031 

MOTS 40" 

Johannesburg, Union of S. Africa 


1GFORK 

1034 

MOTS 40" 

East Grand Forks, Minn. 


1ROSMA 

1042 

MOTS 40" 

Rosman, N. C. 


1TANAN 

1043 

MOTS 40" 

Tananarive, Madagascar 

SPEOPT 

1EDINB 

7036 

MOTS 40" 

Edinburg, Texas 


1COLBA 

7037 

MOTS 40" 

Columbia, Mo. 


1BERMD 

7039 

MOTS 40" 

Bermuda 


1PURIO 

7040 

MOTS 40" 

San Juan, Puerto Rico 


1GSFCP 

7043 

PTH-100 

GSFC, Greenbelt, Md. 


1DENVR 

7045 

MOTS 40" 

Denver, Calif. 


1SUDBR 

7075 

MOTS 40" 

Sudbury, Ontario 


1JAMAC 

7076 

MOTS 40" 

Jamaica, B.W. I. 

SAO 

lORGAN 

9001 

Baker-Nunn 

Organ Pass, N. M. 


10LFAN 

9002 

Baker-Nunn 

Olifantsfontein, Union of S. Africa 


1 SPAIN 

9004 

Baker-Nunn 

San Fernando, Spain 


1NATOL 

9006 

Baker-Nunn 

Naini Tal, India 


1QUIPA 

9007 

Baker-Nunn 

Arequipa, Peru 


1SHRAZ 

9008 

Baker-Nunn 

Shiraz, Iran 


1CURAC 

9009 

Baker-Nunn 

Curacas, Lesser Antilles 


1JUPTR 

9010 

Baker-Nunn 

Jupiter, Fla. 


1VILDO 

9011 

Baker-Nunn 

Villa Dolores, Argentina 


1MAUIO 

9012 

Baker-Nunn 

Maui, Hawaii 


AUSBAK 

9023 

Baker-Nunn 

Woomera, Australia 

NAVY 

ANCHOR 

2014 


Anchorage, Alaska 

TRANET 

DOPPLER 

WAHIWA 

2100 


South Point, Hawaii 


LACRES 

2103 


Las Cruces, N. M. 


LASHM2 

2106 


Lasham, England 


APLMND 

2111 


APL Howard County, Md. 
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Doppler Data 


The GEOS- 1 TRANE T Doppler data used in this study were also obtained from the National 
Space Science Data Center. Data from only the 5 stations listed in Table 2 are used for the period 
under consideration. Data are available from the Doppler station in American Samoa (TAFUNA, 
#2017) but it is not used because of uncertainties in the station position. The data available from 
McMurdo Sound, Antarctica (MCMRDO, #2019) are not used because of the very low maximum ele- 
vation angles on all of the passes. The 5 Doppler stations used are referenced to the SAO-C7 
system. 

All of the Doppler data used in this study are converted from frequency measurements to 
range rate using the following equation for one-way Doppler data: 


R 


c(F b -F m ) 


where 

f b = base frequency 
f m = measured frequency 
c = speed of light (2.997925 x 10 8 m/sec) 


Among the preprocessing done on the data before it was submitted to the Data Center was a 
first-order ionospheric refraction correction applied at the tracking stations (Reference 5). The 
NONAME system is used to further correct on the range rate values for tropospheric refraction 

as follows (Reference 1): 


AR 


2.77N cos E 

s 

328. 5( . 026 + sin E) 2 


meters/second 


where 

N s = (surface index of refraction - 1.) x 10 6 , = 328.5 in the absence of a better value for 
the surface index of refraction 

E - elevation angle computed from the initial estimate of the trajectory 
E = computed rate of change of elevation 

A transit time correction is also applied to the observations to put the time tags at the satellite. 

In addition^ if is felt necessary to include in the NONAME preprocessing an adjustment to the 
base frequency (of the spacecraft oscillator) for each pass over a station. Even though a nominal 
value of the satellite oscillator frequency exists for GEOS-I, it was modified by NWL for each pass 
of data (Reference 6). NWL computed a reference orbit with their ASTRO Computer Program and 
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derived an expected satellite frequency for each observation time. Then (o-c) values were calcu- 
lated and used to produce a corrected nominal satellite oscillator frequency for each pass of data. 

This "base frequency” was included in the Doppler data submitted to the Data Center. Due to the dif- 
ferences (gravity model, station positions, and other constants) between the ASTRO and NONAME 
Orbital Computation Systems, better orbits are obtained in NONAME when the base frequency for 
each pass of Doppler data is adjusted along with the six orbital elements (1 & v). This is done in 
NONAME by adjusting on the range rate measurement bias for each pass of data since a linear re- 
lationship exists between R and f b in the formula for converting Doppler data to range rate. It is 
recognized that residual refraction effects, unmodelled orbital errors, and other small unknown 
biases may be absorbed into this base frequency adjustment. 

III. INTERCOMPARISON OF DOPPLER AND OPTICAL ORBITS 

NQNAME and Perturbations Applied 

Designed to provide accuracy in geodetic studies, the NONAME system at Goddard consists of 
a definitive orbit and geodetic parameter estimation program with a number of auxiliary programs. 
The main program of the system, the NONAME ODP, can operate in either the data reduction or 
orbit generation modes (Reference 2). 

In the data reduction mode, the NONAME ODP can estimate the following parameters from 
satellite tracking data: 

1. the six orbital parameters x, y, z, x, y, z for some specified epoch 

2. certain physical constants relating to atmospheric drag or solar radiation pressure 

3. tracking station co-ordinates relative to the center of mass of the earth 

4. tracking Instrument errors— zero set bias or timing bias 

5. geopotential coefficients 

All observation time tags are transformed to U.T.C. time at the satellite and numerous preproc- 
essing options exist for various data types (i.e. right ascension, declination, range, range rate, di- 
rection cosines, X and Y angles, azimuth, elevation). 

The orbit is numerically integrated (CowelPs method) in fixed steps and interpolated to get 
computed observations for residuals, (o-c) values. A Bayesian least squares estimation scheme 
and a Newton-Raphson iteration formula are used in correcting on the six orbital elements and any 
other specified parameters. There are convergence criteria for the iterations, rejection criteria 
for observations, and observational data weighting schemes optional to the user. 

In the orbit generation mode, an initial epoch and position and velocity vectors are input and 
the equations of motion are numerically integrated (Cowell) to give an ephemeris of position, velocity, 
and time. 

The potential of the earth is represented by a normal potential of an ellipsoid of revolution 
(SAO C-5, C-6, 07, etc.) and small variations, expressed by a set of spherical harmonics (SAG 
M-l, APL 3.5, etc.). In addition, the following perturbations may be represented as disturbing 
functions as optioned by the user: 

1. solar gravitation 

2. lunar gravitation 
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3, solar radiation pressure 

4* atmospheric drag force (NONAME uses the Jacchia-Nicolet model for the atmosphere) 

In this study, the earth ellipsoid used is that of the SAO C-7 system and the gravity model used 
is the SAO- Ml modified by the 12th order terms of Gaposchkin and Yeis (References 2 and 7). The 
perturbations applied in NONAME are solar gravitation (M s /M e = 332951.25), lunar gravitation 
(M m /M e - .0123), and solar radiation pressure (4.5 x 10“ 6 Newtons/m 2 ). 


Discussion of Comparison 

The period of data that was intercompared covers from July 9th through July 26th and from 
July 31st through August 7th of 1966. The data were separated into thirteen 2-day arcs since any 
longer arcs are more affected by errors due to uncertainties in the earth’s gravity model, solar 
radiation pressure or other parameters. 

Both the optical and Doppler data were run separately on the NONAME ODP in the data reduc- 
tion mode. All optical observations were assigned a weight of 2 seconds of arc, and the right as- 
cension measurements were further down-weighted by the cosines of the corresponding declination 
measurements due to the geometry of the pair (Le. the higher the declination, the larger the un- 
certainties introduced into the right ascension). The Doppler data, converted to range rate, were 
assigned a weight of 10 cm/sec. Based on an editing criteria of 3^ there was approximately a 1 % 
rejection rate on both data types. 

The converged solution state vectors for both the optical and Doppler coincident arcs were 
then input into the NONAME ODP in the orbit generation mode. The orbits generated were com- 
pared to obtain the position differences every ten minutes. Table 3 presents the RMS of solution 
values for the data reduction runs and the summary data of the orbit intercomparisons. Appendix B 
shows position difference plots of five of these 2-day arc inter comparisons. 

See Appendix A for a complete breakdown of the optical and Doppler data used in these 2- day 

arc studies. 

As seen from Table 3, the RMS of total position differences over the 2-day arcs is on the order 
of 20 to 40 meters with maximum position differences on the order of 30 to 70 meters. It is 
important to note that only 5 stations in the Northern global hemisphere were used to obtain the 
Doppler orbits and all low elevation data was included, indicating good-quality refraction correc- 
tions. The large amplitude of some of the cross track position differences, might seem odd at first 
glance since the along track errors usually dominate; however, this could possibly be attributed to 
difficulty in determining the inclination of the orbit in the Doppler solution since, again, the five 
Doppler stations are poorly distributed around the globe. These consistent results show that the 
capability to handle Doppler data and determine orbits of the quality demanded in geodetic studies 
from such data exists in the NONAME System at Goddard. 
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Table 3 


Position Differences Between Orbits Generated From TRANET Doppler Data and Orbits 
Generated From Optical Data (2 day arcs). 


Arc of 
Compare 
(1966) 

*RMS of Fit 

RMS of Position Differences (meters) 

— — — — 

Maximum 

Difference 

(meters) 

Doppler 
(cm/ sec) 

Optical 
(secs, arc) 

Radial 

Cross 

Track 

Along 

Track 

Total 

July 9-10 

2.7 

1.8 

12 

7 

38 

41 

70 

July 11-12 

2.8 

1.9 

7 

16 

19 

26 

36 

July 13-14 

2.8 

1-8 

13 

9 

26 

30 

44 

July 15-16 

2.8 

1.8 

10 

23 

23 

33 

51 

July 17-18 

2.7 

1.9 

7 

10 

20 

23 

35 

July 19-20 

2.6 

1.9 

12 

10 

36 

39 

67 

July 21-22 

2.7 

2.1 

12 

10 

38 

41 

77 

July 23-24 

2.7 

1.8 

3 

14 

11 

19 

27 

July 25-26 

2.7 

1.7 

7 

9 

16 

20 

38 

July 31— August 1 

2.8 

1.9 

6 

5 

18 

20 

39 

August 2-3 

3.0 

2.1 

5 

9 

18 

21 

47 

August 4-5 

2.7 

2.0 

6 

25 

32 

41 

66 

August 6-7 

2.7 

1.9 

11 

21 

32 

40 

60 


* Approximate number of observations per arc: Doppler — 1880, Optical — 870. 


IV. EFFECT OF STATION POSITIONS ON THE DOPPLER ORBITAL SOLUTIONS 

As previously mentioned, the 5 Doppler tracking stations are in the Northern hemisphere of 
the globe; Figure 1 shows the locations of ail the optical and Doppler tracking stations from which 
data were used for this study. It is interesting that all optical stations tracked the satellite only as 
it passed from south to north over the station. Thus the South American optical stations were often 
tracking just minutes before the Doppler stations in Maryland and England started tracking. The 
reason for this consistent south to north tracking .by the optical stations lies in the facts that for 
the GEOS- 1 satellite, the right ascension of the ascending node traverses the celestial sphere at a 
rate of approximately 2-1/4 degrees per day, and only night tracking of the flash sequences is pos- 
sible. Thus, for this month long period, the right ascension of the node was in the earth's shadow, 
causing the consistent south to north tracking. The Doppler stations have no such restriction and 
were able to track whenever the satellite passed over their vicinity. 

It was thought that perhaps the poor global distribution of the available Doppler stations caused 
the Doppler orbit to be weakened in the southern hemisphere where no stations existed and 
thus where no (o-c) values could be determined to correct on the state vector. Thus it was decided 
to take 4 of the 2 -day Doppler and optical arcs and examine how accurately the Doppler determined 
orbits fit the optical data from each optical station around the entire globe. 
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D NAVY TRANET DOPPLER STATION © OPTICAL STATION ( SAO, SPEOPT, OR MOTS ) 



Figure 1 — Locations of all TRANET Doppler and optical stations used in intercomparisons 
and typical GEOS-I orbital paths July-August, 1966. 


The approach is to take the converged vector solution, for each of the arcs, as determined by 
the Doppler arc and input it as the state vector for the corresponding optical arc. The Doppler 
solution is passed through the optical data to note the optical station data fits to the Doppler de- 
termined orbit. The observation residual summary by optical station is examined to see which 
groups of stations representing areas of the globe have the largest increases in the RMS of solu- 
tion when compared to the final iteration of the original optical solution for that arc. For each 
station appearing in two or more of the four arcs, the RMS’s of fits were combined to yield one 
summary value for that station. As an example, if a station appeared in two of the arcs and had 
respective RMS values of RMS X and RMS 2 , where: 



then the combined RMS would equal: 

/(m - 1 )(KMS 1 ) 2 + (k - 1)(RMS 2 ) 2 

y m + k - 1 


I 36 


Table 4 shows all of the optical stations used ? the number of observations each has over the four 
2-day arcs of this station position study, what arcs the data appear. in, the combined EMS values 

Table 4 


Summary, by Station, of Optical Data Fits to the Optical and Doppler Orbital Solutions (July 9-16, 1966), 


Optical Station 

Latitude 

(degrees) 

Observations 

RMS of Fit (secs, arc) 

Total Number 

*From Arcs 

To Optical Orbit 

To Doppler Orbit 

1GFOKK 

48 

73 1 

i, n 

1.9 

2.0 



73 * 


1.5 

1.7 

1SUDBK 

46 

60 

II, IV 

1.9 

2,3 



60 


1.5 

1.9 

1DENVR 

40 

20 

I 

2.4 

2,7 



20 


1.3 

1.2 

1COLBA 

39 

93 

1,11 

1.3 

1.6 



94 


1.2 

1.7 

1BPOIN 

38 

45 

i.n 

1.3 

1.5 



44 


2.1 

2.1 

1 SPAIN 

36 

73 

I, II, XII, IV 

2.7 

3.3 



77 


1.7 

2.1 

1MOJAV 

35 

35 

n 

1.8 

2.4 



35 


1.5 

1.9 

1ROSMA 

35 

33 

i, n 

1.3 

1.3 



32 


2.3 

2.2 

10RGAN 

32 

224 

i, n, m, iv 

1.6 

1.6 



224 


1.6 

1.8 

MPsbebbi^ 

32 

49 

n, iv 

2.6 

2,6 



48 


3.2 

2,5 

1SHRAZ 

30 

28 

n, m 

1.8 

2.1 



25 


2.8 

2.4 

1NATOL 

29 

25 

in 

1.4 

1.4 



25 


1.7 

2,3 

1JUPTR 

27 

305 

i, n, m, iv 

1.7 

1.8 



304 


1.8 

2,3 

1EDINB 

26 

69 

i.n 

0.7 

1.2 



70 


1.1 

1.6 

1MAUIO 

21 

56 

I, IV 

2.2 

2,5 



56 


1.4 

1.4 

1JAMAC 

18 

35 

n 

1.3 

1.4 



34 


1.9 

2,6 

1CURAC 

12 

27 

i 

1.7 

1.6 



27 


2.3 

2,6 

1QUIPA 

-16 

69 

i, n, in, iv 

2.0 

2.0 



69 


2.7 

3.6 ! 

1TANAN 

-19 

14 

n 

1.8 

1.6 



14 


0.9 

1.2 1 


*1 - July 9-10, II - July 11-12, III - July 13-14, and IV - July 15-16. 

tdeclmation 

bright ascension 
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Table 4 (Continued) 


j Optical Station 

Latitude 

(degrees) 

Observations 

RMS of Fit (secs arc) 

Total Number 

*From Arcs 

To Optical Orbit 

To Doppler Orbit 

1JOBUR 

-26 

61 t 

I, IX, in, IV 

1.7 

3.4 

j 


.52* 


2.0 

2.9 

10LFAN 

-26 

132 

i, n, m, iv 

1.7 

3.4 



130 


2.1 

3.0 

lOOMER 

-31 

42 

i» n 

1.8 

3.0 

' 

! 

41 


2.1 

3.6 

AUSBAK 

-31 

69 

i, n, in, iv 

1.6 

2.5 



69 


2.2 

2.8 

i XVI LOO 

-32 

31 

II, IV 

3.0 

3.5 

!. 


33 


3.3 

3.9 


*1 * July 9-10, II - July 11-12, III - July 13-14, and IV - July 15-16. 

tdeclinaeion 
bright ascension 


4.0 


a 

Z 

O 

o 



3,6 

3.2 
2.8 
2.4 
2.0 
1 a 6 

1.2 
0.8 
0.4 


0.0 


z c 
o „ 

P C 
< 2 
Z C 
z c 

y * 

O 2 



m DOPPLER BASED RESIDUALS EXCEED 
11! OPTICAL RESIDUALS 

B OPTICAL RESIDUALS EXCEED DOPPLER 
BASED RESIDUALS 



OPTICAL STATIONS 


Figure 2 ■“■“-Summary of optical residuals based upon 
TRANET Doppler orbits for the period July 9-16, 1966. 


from the original optical solutions, and the com- 
bined RMS values when the Doppler state vector 
was passed through the optical data. The stations 
are grouped according to hemispheric location. 
By comparing the optical based RMS values 
against the Doppler based RMS values for a given 
station, some feeling for the relative strengths of 
the two solutions can be obtained. 

Figure 2 presents the RMS differences in 
graphical form. For a few stations, the Doppler 
orbit fits the optical data slightly better than the 
optical orbit. This could be due to errors in the 
station position, instrument errors, or other 
factors. However, the general trend for northern 
hemisphere optical stations is for the optical 
orbit to fit the optical data a little better than 
the Doppler orbit. In the southern hemisphere, 
the Doppler orbit has obvious difficulty fitting 
the optical data. Thus, it seems that the good 
orbital comparisons that were achieved would 
have been even better if data were available 
from Doppler stations in the southern hemis- 
phere to tie down that half of the GEOS- 1 orbit. 

Figure 2 should be interpreted with the in- 
formation presented in Table 4 in mind. As an 
example, the histogram for the station in 
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Madagascar (1TANAN) seems to disagree with those of the two nearby South African stations in 
Olifantsfontein and Johannesburg. However, data from Tananarive consists of 28 observations in 
only one arc while there are 113 observations from Johannesburg and 262 from Olifantsfontein and 
these data are distributed through all 4 arcs. Thus, this indicates that the Tananarive results 
should be discounted. 


V. FURTHER STUDY OF THE INTERCOMPARISONS 


Evaluation of Atmospheric Refraction Corrections on Low Elevation Doppler Observations 

It was thought that the optical and Doppler solutions might be brought into even better agree- 
ment by dropping the low elevation measurements from the Doppler data. The logic behind this is 
that perhaps the refraction effects, which are greatest at low elevations, are not sufficiently 
modelled. Additional NONAME solutions were generated with all Doppler observations of eleva- 
tion less than 20° dropped from the data; two 2-day arcs were tested. Table 5 indicates a very 
slight drop in the RMS of solution and shows the results of the position comparisons against the 
corresponding optical arcs. These "elevation cut-off" Doppler versus optical comparisons show an 


Table 5 

Effects of Removing Low Elevation Observations From Doppler Data on the Doppler 
Solutions and on the Doppler Versus Optical Orbital Comparisons.^ 


Navy TRANET Doppler Solutions 

2 -Day 
Arc 
1966 

Station 

Total 

Number of 
Observations 
Used 

Total 

Number of 
Observations >20° 
Elevation Used 

Original 
RMS of 
Fit 

(cm/ sec) 

RMS of 
Fit 
(>20°) 
(cm/sec) 

July 17-18 

ANCHOR 

429 

272 

3.1 

3.0 


WAfflWA 

227 

111 

2.5 

2.6 


LACRES 

387 

187 

2.9 

2.8 


LASHM2 

447 

298 

a.s 

2.2 


APLMND 

423 

266 

2.4 

2.3 


ALL 

1913 

1134 

2.7 

2.6 

July 19-20 

ANCHOR 

522 

303 

2.8 

2.6 


WAHIWA 

234 

149 

2.7 

2.5 


LACRES 

307 

190 

2.7 

2.3 


LASHM2 

410 

277 

2.1 

2.0 


APLMND 

483 

271 

2.6 

2.5 


ALL 

1956 

1190 

2.6 

2.4 
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Table 5 (Continued) 


Doppler Versus Optical Orbital Comparisons 


2 - Day Arc 
of Compare 
(1966) 

RMS of Fit 

RMS of Position Differences (meters) 

Maximum 

Difference 

(meters) 

Doppler 

(cm/sec) 

Optical 
(secs, arc) 

Radial 

Cross 

Track 

Along 

Track 

Total 

July 17-18 

2.6* 


7 

11 

21 

24 

38 



1.9 







2.7 


7 

10 

20 

23 

35 

July 19-20 

2.4* 


13 

10 

37 

41 

68 



1.9 







2.6 


12 

10 

36 

40 

67 


•Observations with > 20° elevation only. 

tNjote: the inclusion of low elevation Doppler data does not degrade the solution. 


increase in the total RMS of position differences of about 1 meter over the original Doppler versus 
optical comparisons presented in Table 3. Thus it is felt that the inclusion of low elevation Doppler 
data does not degrade the solution. In fact, the low elevation points add to the geometryof the Doppler 
solutions since they increase the tracking range in the GEOS -I satellite by approximately 60 miles/ 
degree of elevation on both ends of any specific pass. 


Effect of Arc Length on Solution 


Two day arc lengths may be less accurate than shorter arc lengths due to the effects of data 
distribution, gravity model errors, perturbation model errors, and unmodelled parameters. How- 
ever, two day arc lengths were chosen for this 
study because it was thought that none of 
these errors would have an adverse effect 
over a 2-day period. This section presents 
the results of a check made to determine 
the quality of 2-day versus shorter length 
orbital solutions. 


ARC 

LENGTH 


1st 

DAY 


1 1 / 2 DAYS 
11/2 DAYS 
1 DAY 
1 DAY 
1 / 2 ' DAY 
1/2 DAY 
1/2 DAY 
1/2 DAY 


NOON 


COVERAGE 

2nd 

DAY 


NOON 


* EIGHT ARCS WERE DETERMINED FOR EACH OF TWO 
2 DAY PERIODS, JULY 17-18, JULY 19-20. 


Figure 3— Reduction of two day arcs into 
shorter arc lengths.* 


In order to assess the accuracy of orbital 
arc solutions of less than 2-day length, two arcs 
(July 17-13, July 19-20) are broken into smaller 
arcs of length 1-1/2, 1, and 1/2 day length as 
shown in Figure 3. Doppler and optical orbits 
are determined for these 16 shorter arcs and 
then inter compared as in the initial 2- day orbi- 
tal solutions. 
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The results of 9 of these 16 inter comparisons along with the results of the two 2-day inter - 
comparisons are presented in Table 6. Seven results are not presented since very poor data dis- 
tribution over the arc length prevented the determination of sufficiently accurate orbits. The re- 
sults indicate that the fit of the data to the orbit does not appreciably improve as the arc length of 
solution decreases; thus, the orbital solutions of length less than 2-days are not of better quality 
than the 2 -day solutions. 


Table 6 

Position Differences Between TRANET Doppler Orbits and Optical 
Orbits for Arc Lengths Less Than 2 Days. 


Arc of Compare 

RMS of Fit 

RMS of Position Differences (meters) 


Epoch (1966) 


Doppler 

Optical 

Radial 

Cross 

Along 

Total 

Maximum 

Difference 

Day 

Hour 

Bj | 

(cm/sec) 

(secs, arc) 

Track 

Track 


July 17 

0 

2 

2.7 

1.9 

7 

mm 

20 

23 

35 

July 17 

0 

1-1/2 

2*7 

1.9 

5 

9 

18 

21 

34 

July 17 

12 

1-1/2 

2.7 

1.9 

8 

2 

27 

29 

54 

July 17 

0 

1 

2.5 

1.8 

8 

11 

23 

27 

53 

July 18 

0 

1 

2.8 

1.9 

8 

8 

22 

25 

49 

July 17 

0 

1/2 

2.5 

1.8 

1 

16 

7 

17 

26 

July 18 

0 

1/2 

2.6 

1.9 

4 

2 

13 

13 

26 

July 19 

0 

2 

2.6 

1.9 

12 

10 

36 

39 

67 

July 19 

0 

1-1/2 

2.7 

1.8 

10 

11 

29 

32 

53 

July 19 

0 

1/2 

2.5 

1.6 

11 

6 

38 

40 

65 

July 20 

! 

0 

1/2 

2.2 

1.7 

i 

12 

11 

33 

37 

62 


The position comparisons of the Doppler and optical orbits are generally in the same range. 

The two comparisons having very low RMS of position differences over the half day solutions on 
July 1 7th and July 18th are the results of very dense data in the period. Each of the thirteen 2 -day 
arcs composed of optical data have an average of 870 observations in it, but both of the above men- 
tioned half day arcs have over 400 optical observations in them. They also have over 400 Doppler 
observations. These comparisons axe further evidence of the quality of the Doppler orbits. 

The plots of position differences for these 9 shorter arcs appear in Appendix C. 

VL BIAS STUDY OF DOPPLER DATA 
Timing Errors 

One of the ancillary programs in the NONAME system processes the residuals calculated in the 
last iteration of a NONAME ODP data reduction run to determine zero-set and timing errors (Refer- 
ence 1). ' • 
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Table 7 


The residuals are investigated with the fol- 
lowing regression model: 


Summary of Timing Biases (At)* Found 
in the TRANET Doppler Data 


AS = AB + At 6 

where 

AR = the residual for a specific observation. 

AB = the zero-set error in the observing 
instrument. 

At = the timing error in the observing 
instrument. 

6 = the rate of change of the observation. 

This program, GEORGE, was used to de- 
termine the timing errors for the July 17-18 
and July 19-20 Doppler arcs. Table 7 sum- 
marizes the results which appear in full in Appendix D. Although zero-set errors were computed 
by GEORGE, they are neglected because such errors were absorbed into the adjustment of the 
range rate measurement bias (the base frequency adjustment). 

It is obvious from the table that no significant timing errors were found. 


2 -Day 
Arc 
(1968) 

. 

Station 

Number of 
Passes 

Timing 
Error (At) 

± 1 Standard 
Deviation 
(msecs.) 

July 17-18 

APLMND 

10 

0.6 ± 4.2 


LASHM2 

10 

0.0 ± 3.2 


ANCHOR 

10 

-0.8 ± 3.2 


LACHES 

9 

-0.4 ± 5.5 


WAHTWA 

5 

2.7 ± 5.5 

July 19-20 

APLMND 

11 

0.4 ± 5.8 


LASHM2 

9 

-2.1 ± 3.0 


ANCHOR 

12 

-1.7 ± 2.7 


LACHES 

7 

2.9 ± 1.5 


WAHIWA 

5 

3.5 ± 6.3 


*AR - AB + At 6 where AR " Residual (o-c), AB = Zero Set 
Bias* At ~ Timing Error, and O - Rate of Change of 

Observation. 


Results of the Range-Rate Bias Adjustment on TRANET Doppler Data 

As mentioned previously 3 the NONAME Orbit and Geodetic Parameter Estimation System is 
used to adjust on a range-rate or base frequency Mas for each pass of data in all Doppler orbital 

solutions (refer to Section II), 


Table 8 


Summary of Biases Found in the TRANET Doppler 
Data Over the Period July 9-28 
July 31 -August 7, 1968. 


Station 

Number of 
Passes 

Range Rate 
Bias ± 1 Standard 
Deviation 
(cm/sec) 

APLMND 

140 

10.2 ± 3.4 

LASHM2 

138 

8.2 ± 2.6 

ANCHOR 

138 

8.3 dfc 3.1 

LACHES 

94 

8.2 ± 3.2 

WAHIWA 

60 

9.3 ± 2.4 

ALL 

. J 

570 

8.8 ± 3.1 


The bias values that are determined in the 
NONAME System are consistent throughout all 
of the Doppler orbital solutions. Table 8 pre- 
sents a summary of all range-rate biases de- 
termined over the entire period of the study. 
Appendix D presents the determined biases in 
more detail. The biases for each station gen- 
erally range from 8 to 10 cm/sec with stand- 
ard deviations of 2 to 3 cm/sec. 

Analysis of simultaneous GEOS- II Doppler 
and laser data at Wallops Island by Herbert 
and Parker (Reference 8) also indicated the 
presence of a positive bias in the Doppler data. 
Using 10 passes of laser data to determine 
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reference orbits, Berbert and Parker noted TRANE T range rate biases averaging +16 cm/sec. 
Discussions between the investigators and NWL personnel uncovered a procedure in the preproc- 
essing program at NWL which was responsible for the large positive biases. After NWL provided 
corrections to the base frequency, the average bias dropped to +4 cm/sec. 


Effect of Arc Length on the Range-Rate Bias Adjustment 

The question of independence of the range-rate bias adjustment on arc length is investigated 
Two of the 2-day arcs, July 17-18 and July 19-20 are broken into smaller arcs of 1-1/2, 1, and 1/2 
day lengths (see Figure 3); 8 arcs are thus formed in each 2-day period. If a particular pass lasts 
from morning into the afternoon (G.M.T.), then only two of these shorter arcs can be used to get a 
range-rate bias for that pass, the reason being that the pass must fall completely within an arc in 
order to get a bias adjustment for the entire pass. For passes within the first and last quarters of 
the 2-day period, 3 of these shorter arcs can be used to determine a range-rate bias, and for passes 
within the center quarters of the 2-day period, 4 shorter arc solutions for range -rage bias are possible. 

Table 9 


Dependency of the TRANET Doppler Range Rate Bias Adjustment on Arc Length. 


Pass for Which 
Bias is Computed 

Arcs Used 

Adjusted 
R Bias 
(cm/sec) 

Mean Bias 
± 1 Standard 
Deviation 
(cm/sec) 

Station 

Day 

Epoch (1966) 

Length 
in Days 



Day 

Hour 

APLMND 

July 17 

0 

2 

5.6 


July 17 

July 17 

0 

1-1/2 

5.6 

C A J. A C 

9:48 

10:17 

July 17 

0 

1 

5.9 

O.U ± U.D 



July 17 

0 

1/2 

6.7 


ANCHOR 

July 17 

0 

2 

7.5 


July 17 

July 17 

0 

1-1/2 

7.7 


15:49 

16:18 

July 17 

12 

1-1/2 

6.5 

7.1 ± 0.5 



July 17 

0 

1 

6.8 




July 17 

12 

1/2 

6.9 


LASHM2 

July 17 

0 

2 

8.3 


July 18 

July 17 

0 

1-1/2 

7.8 


3:48 

4:17 

July 17 

12 

1-1/2 

8.0 

7.8 ± 0.4 



July 18 

0 

1 

7.7 




July 18 

0 

1/2 

7.3 


WAfflWA 

July 19 

0 

2 

12.4 

1 

July 19 

July 19 

0 

1-1/2 

12.7 

i 

16:06 

16:35 

July 19 

12 

1-1/2 

11.9 

11.8 ± 0.8 



July 19 

0 

1 

11.3 




July 19 

12 

1/2 

10.7 

I 

LACRES 

July 19 

0 

2 

11.9 


July 20 

July 19 

12 

1-1/2 

12.5 


12:04 

12:33 

i 

0 

1 

12.7 

11,9 db 0.9 




12 

1/2 

10.6 
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Table 9 presents the results for 5 representative passes to show the general consistency of the 
range -rate or base frequency adjustment. Appendix E presents the results for all 88 passes cover- 
ing the 4 days of July 17-20, 1966. Most of the passes show a mean bias value with a standard devia- 
tion of about 2-10% of the mean. This result indicates a strong level of independence between arc 
length and the Mas adjustment on range-rate, and it provides further evidence that the TRANET 
Doppler data are being processed properly in the NONAME System. 


VIL CONCLUSIONS 

The thirteen 2-day optical orbits and the thirteen 2-day Doppler orbits have RMS of fits in the 
ranges of 1.9 seconds of arc and 2.7 cm/sec respectively, indicating for both data types that the 
fits of the data to the orbits are almost down to the noise level of the data. The RMS of position 
differences between corresponding optical and Doppler orbits range from 20-40 meters, and max- 
imum position differences are from 27-77 meters for the tMrteen 2-day arcs. Since the noise 
level on both data types is from 10-15 meters, it is felt that these results are very consistent and 
that they demonstrate the ability to generate Doppler orbits of quality comparable to the optical 
orbits that are being used in geodetic studies at Goddard. 

Furthermore, the study of the effects of station positions on the Doppler orbital solutions in- 
dicates that even better agreement between the Doppler and optical orbits would have been attained 
if there had been data available from Doppler stations located in the southern hemisphere to tie 
down those halves of the Doppler orbits. 

It is felt that the NONAME ODP models tropospheric refraction well enough to enable the use 
of Doppler observations far below 20° in elevation. The use of observations in the range of lG°-20° 
elevation increases the geometry of any pass of the GEOS-1 satellite by approximately 60 miles per 
degree of elevation on either end of the pass; thus, this ability is quite valuable in determining 

orbits of geodetic quality. 

Biases in the base frequency values sent along with each pass of Doppler data were consistently 
on the order of 8-10 cm/sec ±2-3 cm/sec for each of the 5 stations used in this study. Also, the 
range-rate biases adjustment in the NONAME System was found to be highly independent of arc 
length of solution. No significant timing biases were found in the Doppler data. 

Finally, this report indicates that Navy TRANET Doppler data can currently be used to sup- 
plement optical data in geodetic studies using the NONAME Orbit and Geodetic Parameter Estima- 
tion System at Goddard or in a system having similar capacity to handle Doppler data. 

Work is currently in progress performing additional studies of this nature using data from 

other satellites (GEOS-H, BEB, BEC). 
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Appendix A 

Optical and Doppler Data Used in the Analysis 


The tables herein present the thirteen 2-day optical orbital solutions (Table Al) and the thirteen 
2-day TRANE T Doppler orbital solutions (Table A2) in detail. The number of observations from 
each station that were used to determine the orbit and the RMS of fit of the orbit to the data from 
each station is shown. 


Table Al 


2-Day Optical Orbital Solutions. 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

■ 

RMS of Fit 
(secs, arc) 

July 9-10 

AUSBAK 

14* 


1*5 



14 1 


2.2 


1BPOIN 

20 


1.5 



19 

1 

2.3 


ICO LB A 

28 


i.i 



28 


1.3 


1CURAC 

27 

1 

1.7 



27 

1 

2.3 


1DENVR 

21 


2.4 



21 


1.3 


1EDINB 

20 

1 

0.7 



21 


1.3 


1GFORK 

25 


2.1 



25 


2.0 I 


1GSFCP 

6 


- 



1 

5 

- 


1JOBUR 

21 


2.0 



18 

3 

2.5 


1JUPTR 

35 


1.8 



35 


1.3 


1MAUIO 

28 


2.0 



28 


1.5 


lOLFAN 

14 


2.0 



14 


!.6 | 


*First line is for declination observations. 
tSecond line is for right ascension observations. 
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Table A1 (Continued) 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(secs, arc) 

July 9-10 

lOOMEK 

21 


2.0 

(continued) 


21 


2.1 


10RGAN 

56 


1.8 



56 


1.6 


1QUIPA 

14 


1.2 



14 


3.7 


1ROSMA 

21 


1.3 



20 

1 

2.5 


ISP AIN 

14 


1.2 



14 


3.7 


1TANAN 

7 


- 



7 


- 


ALL 

775 

13 

1.8 

July 11-12 

AUSBAK 

21 


1.4 



21 


2.4 


1BERMD 

35 


2.9 



35 


2.7 


1BPGIN 

25 


1.2 



25 


2.0 


1COLBA 

65 

1 

1.4 



66 


1.2 


1EDINB 

49 


0.7 



49 


1.1 


1GFORK 

48 


1.8 


i 

48 


1.2 


1 JAM AC 

35 


1.3 



34 

1 

1.9 


1JOBUR 

14 


1.6 



10 

4 

1.4 


1JUPTR 

55 


2.1 



55 


2.4 


1MAUIO 

3 


_ 



3 


- 


1MOJAV 

35 


1.8 



35 


1.5 


XNATOL 

7 


- 



6 

1 

- 


lOLFAN 

21 


1.7 



20 

1 

2.8 
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Table A1 (Continued) 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(secs, arc) 

July 11-12 

lOOMER 

21 


1.7 

(continued) 


20 

1 

2.2 


1 ORGAN 

35 


1.5 



35 


1.6 


1PURIO 

7 


- 



7 


- 


1QUIPA 

14 


1.9 



14 


3.1 


1ROSMA 

12 


1.2 



12 


1.8 


1SHRAZ 

14 


2.2 



11 

3 

3.0 


ISP AIN 

24 

4 

3.6 



28 


2.1 


1SUDBR 

46 


2.0 



46 


1.7 


1TANAN 

14 


1.8 



14 


0.9 


1VILDO 

17 

4 

3.8 



19 

2 

3.6 


ALL 

1230 

22 

1.9 

July 13-14 

AUSBAK 

20 


1.8 



20 


1.3 


1JOBUR 

12 

1 

1.7 



10 

3 

1.8 


1JUPTR 

105 


1.7 



104 

1 

1.7 


1NATOL 

25 


1.4 



25 


1.7 


lOLFAN 

48 

1 

1.8 



47 

2 

2.0 


lOOMER 

7 


- 



7 


- 


lORGAN 

59 

2 

1.9 



59 

2 

1.7 


1QUIPA 

20 


2.7 



20 


2.0 


1SHRAZ 

14 


1.5 



14 


2.8 | 
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Table A1 (Continued) 


2 -Day Arc 


Station 


Number of 
Observations Used 


Number of 
Observations 
Rejected 


RMS of Fit 
(secs, arc) 


July 13-14 
(continued) 


July 15-16 


1 SPAIN 
1TANAN 
1VILDO 
ALL 

21 

21 

7 

7 

6 

0 

678 

AUSBAK 

14 


14 

1BERMD 

14 


13 

1CURAC 

6 


6 

1JOBUR 

14 


14 

1JUPTR 

no 


no 

1MAUIO 

28 


28 

10LFAN 

49 


49 

lORGAN 

74 


S3 




Table A1 (Continued) 


2 -Day Arc 

Station 



Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(secs, arc) 

July 17-18 

AUSBAK 

49 


1.7 



46 

3 

2.3 


1BERMD 

42 


2*9 



41 

1 

3,0 


1BPOIN 

12 


1,4 



12 


2,7 


1COLBA 

28 


1.3 



28 


1.5 


1EDINB 

35 


1.1 



35 


1.8 


1GFORK 

18 


2,5 



17 

1 

2.1 


1JAMAC 

13 


1.2 



13 


1.7 


1JOBUR 

13 


2,1 



13 


1.5 


1JUPTR 

112 


1.9 



112 


1.9 


1MAUIO 

29 


1.4 



29 


1.6 


1MOJAV 

14 


2,0 



14 


1.1 


lOLFAN 

35 


1.8 



35 


1.9 


10RGAN 

71 

2 

1.6 



71 

2 

1.4 


1PURIO 

13 


1.0 



12 

1 

3,3 


1QUIPA 

21 


1.7 



21 


2,8 


ISP AIN 

47 

2 

1.9 



47 

2 

2,0 


1TANAN 

14 


2,3 



14 


2,6 


1VXLDO 

14 


1.4 



14 


2,0 


ALL 

1154 

14 

1.9 
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Table A1 (Continued) 


2 -Day Arc 

Station 

Number of 

! Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(secs, arc) 

July 19-20 

AUSBAK 

6 





6 


- 


1BERMD 

7 


— 



7 


- 


1BPOIN 

33 


1.7 



33 


1.8 


1COLBA 

43 


1.0 



43 


1.4 


1EDINB 

56 


1.0 



56 


1.3 


1FTMYR 

•7 





7 


- 


1GFORK 

62 


1.5 



62 


1.3 


1GSFCP 

7 


_ 



7 


- 


1JAMAC 

14 


0.6 



14 


2.6 


1JOBUR 

32 

3 

1.9 



33 

2 

1.8 


1JUPTR 

62 

1 

2.3 



62 

1 

2.3 


1MAUIO 

12 


2.6 



12 


2.6 


1MOJAV 

35 


1.3 



35 


1.9 


lOLFAN 

21 


1.6 



21 


2.3 


1 ORGAN 

8 


_ 



8 


- 


1QUIPA 

27 

1 

3.0 



25 

3 

2.8 


1 SPAIN 

75 

1 

2.7 



75 

1 

1.9 


1SUDBR 

42 


1 2.1 



41 

1 

2.0 


1VILDO 

7 





1 

6 

- 


ALL 

1104 

20 

1.9 
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Table AX (Continued) 


2-Dav Are 


July 21-22 



AUSBAK 


Number of 
Observations Used 


Number of 
Observations 
Rejected 


RMS of Fit 
(secs, arc) 


1COLBA 


1EDINB 


1FTMYR 


1JOBUR 


1JUPTR 


1MAUIO 


lOLFAN 


10RGAN 


1QUIPA 


1ROSMA 


1 SPAIN 


1SUDBR 


1VILDO 


AUSBAK 


1BERMD 


1BPOIN 


1DENVR 


1EDINB 
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Table A1 (Continued) 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(secs, arc) 

July 23-24 

1FTMYR 

7 


- 

(continued) 


7 


“ 


1GFORK 

26 


2.2 



26 


1.5 


1GSFCP 

7 


- 



7 


- 


1JAMAC 

14 


1.6 



14 


2.2 


1JOBUR 

21 


1.1 



21 


1.5 


1MAUIO 

8 


- 

j 


8 


- 

i 

1MOJAV 

14 


1.5 



14 


2.1 


10LFAN 

21 


1.6 



21 


2.0 


lORGAN 

132 


1.7 



132 


1.6 


1PURIO 

14 


1.9 



14 


1.6 


1QUIPA 

28 


1.6 



28 


2.8 


1ROSMA 

38 


1.0 



38 


1.4 


1 SPAIN 

14 


1.8 



14 


1.7 


1SUDBR 

21 


2.3 



21 


1.8 


1VILDO 

21 


1.7 

1 


21 


2.2 

| 

1 

ALL 

989 

25 

1.8 

July 25-26 

AUSBAK 

35 


1.8 



33 

2 

2.3 


1BPOIN 

7 


- 



7 


- 


1COLBA 

21 


0.8 



21 


1.2 


. 1DENVR 

19 


1.4 



19 


1.3 
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Table A1 (Continued) 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(secs, arc) 

August 2-3 

JOHNST 

0 

2 

- 



1 

1 



1BERMD 

7 


- 



4 

3 

- 


1BPOIN 

21 


1.5 



21 


1.9 


1COLBA 

63 


2.0 



63 


1.2 


1FTMYR 

28 


1.7 



28 


1.5 


1GFORK 

42 


1.9 



42 


1.5 


1GSFCP 

19 


2.1 



17 

2 

3.8 


1JUPTR 

56 


1.6 



56 


1.4 


1MAUIO 

27 


2.0 



27 


2.6 


1MOJAV 

14 


2.4 



14 


2.0 


ISP AIN 

63 


2.3 



63 


2.5 


1SUDBR 

53 


2.1 



51 

2 

2.8 


ALL 

780 

10 

2.1 

August 4-5 

AUSBAK 

35 


2.3 



35 


1.7 


1BERMD 

7 


- 



0 

7 

- 


1COLBA 

41 


2.3 



41 


2.9 


1EDINB 

21 


1.6 



21 


1.9 


1GFORK 

32 


1.7 



32 


3.4 


1JOBUR 

21 


1.5 



21 


1.6 


1JUPTR 

75 


2.2 



74 

1 

1.9 
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2 -Day Arc 


Station 


Table A1 (Continued) 

Number of 
Observations Used 


Number of 
Observations 
Rejected 


RMS of Fit 
(secs, arc) 


August 4-5 
(continued) 


August 6-7 


1MAU10 

53 


53 

1MOJAV 

27 


27 

lOLFAN 

21 


21 

lORGAN 

56 


56 

1QUIPA 

14 


14 

1 SPAIN 

90 


90 

1VILDO 

7 


7 

ALL 

992 

AUSBAK 

14 


14 

JOHNST 

0 


1 

1BERMD 

14 


11 

1BPOIN 

28 


28 

ICO LB A 

35 


37 


1DENVR 


1EDINB 


1GFORK 


1JAMAC 


1JOBUR 


1MOJAV 


lOLFAN 



Table A1 (Continued) 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(secs, arc) 

August 6-7 

lORGAN 

7 


- 

(continued) 


7 


- 


1PURIO 

4 


- 



4 


- 


1QUIPA 

14 


2.2 



14 


2.1 


1 SPAIN 

63 


1.8 



62 

1 

1.8 


1SUDBR 

17 


1.6 



17 


2.0 


1VILDO 

14 


2.1 



14 


1.9 


ALL 

571 

19 

1.9 


Table A2 


2 -Day TRANET Doppler Orbital Solutions. 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(cm/sec) 

July 9-10 

ANCHOR 

363 

4 

2.9 


WAfflWA 

214 

6 

2.8 


LACHES 

352 

2 

2.5 


APLMND 

500 

3 

2.8 


LASHM2 

548 

3 

2.3 


ALL 

1977 

18 

2.7 

July 11-12 

ANCHOR 

449 

11 

3.2 


WAHIWA 

98 

0 

2.4 


LACRES 

420 

2 

2.7 


APLMND 

409 

4 

2.9 


LASHM2 

515 

5 

2.5 


ALL 

1891 

22 

2.8 

July 13-14 

ANCHOR 

522 

7 

3.1 


WAHIWA 

154 

1 

2.6 


LACRES 

326 

5 

2.6 


APLMND 

321 

5 

2.9 


LASHM2 

535 

5 

2.6 


ALL 

1858 

23 

2.8 
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Table A2 (Continued) 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RMS of Fit 
(cm/sec) 

July 15-16 

ANCHOR 

494 

3 

2.9 


WAHIWA 

182 

3 

2.9 


LACRES 

274 

3 

3.0 


APLMND 

487 

7 

2.9 


LASHM2 

528 

2 

2.4 


ALL 

1965 

18 

2.8 

July 17-18 

ANCHOR 

429 

25 

3.1 


WAHIWA 

227 

3 

2.5 


LACRES 

387 

5 

2.9 


APLMND 

423 

0 

2.4 


LASHM2 

447 

3 

2.5 


ALL 

1913 

36 

2.7 

July 19-20 

ANCHOR 

522 

13 

2.8 


WAHIWA 

234 

3 

2.7 


LACRES 

307 

2 

2.7 


APLMND 

483 

c 

• t4 ' 

2.6 


LASHM2 

410 

2 

2.1 


ALL 

1956 

25 

2.6 

July 21-22 

ANCHOR 

440 

11 

2.8 


WAHIWA 

192 

0 

3.0 


LACRES 

429 

3 

2.6 


APLMND 

520 

2 

2.8 


LASHM2 

307 

41 

2.0 


ALL 

1888 

57 

2.7 

July 23-24 

ANCHOR 

440 

11 

3.2 


WAHIWA 

287 

2 

2.7 


LACRES 

328 

2 

2.6 


APLMND 

467 

3 

2.6 


UfSHM2 

275 

2 

2.3 


ALL 

1797 

20 

2.7 

July 25-26 

ANCHOR 

375 

15 

3.2 


WAHIWA 

284 

1 

2.5 


LACRES 

208 

7 

2.9 


APLMND 

494 

31 

2.6 


LASHM2 

501 

3 

2.1 


ALL 

1862 

57 

2.7 

July 31 - 

ANCHOR 

438 

13 

3.3 

August 1 

WAHIWA 

316 

1 

2.3 


LACRES 

160 

2 

3.1 


APLMND 

484 

4 

2.8 


LASHM2 

471 

4 

2.1 


ALL 

1869 

24 

2.8 
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Table A2 (Continued) 


2 -Day Arc 

Station 

Number of 
Observations Used 

Number of 
Observations 
Rejected 

RtyS of Fit 
(cm/sec) 

August 2-3 

ANCHOR 

517 

6 

3.3 


WAfflWA 

268 

1 

2.5 


LACRES 

158 

2 

2.6 


APLMND 

519 

3 

3.1 


LASHM2 

391 

36 

2.7 


ALL 

1853 

48 

3.0 

August 4-5 

ANCHOR 

491 

20 

3.1 


LACRES 

407 

5 

2.8 


APLMND 

392 

0 

2.9 


LASHM2 

532 

4 

2.1 


ALL 

1822 

29 

2.7 

August 6-7 

ANCHOR 

442 

13 

3.1 


WAHIWA 

180 

0 

2.4 


LACRES 

174 

40 

2.8 


APLMND 

511 

4 

2.8 


LASHM2 

541 

3 

2.2 


ALL 

1848 

60 

2.7 
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Appendix B 


Position Differences Between Doppler and Optically Determined Orbits 


Figures are presented for five of the thirteen 2-day orbital intercomparisons showing plots of 
the satellite position differences between the optically determined and Doppler determined orbits 
over the span of the arc. 



CROSS TRACK DIFFERENCES 



Figure Bl— Position differences between TRANET Doppler orbit and optica! orbit 
for the 2-day arc, July 15-16, 1966. 
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METERS METERS 










HOURS FROM EPOCH 


Figure B2 (continued)— Position differences between 7RANET Doppler orbit and optical orbit 

for the 2-day arc, July 17-18, 1966. 











METERS METERS 



12,0 13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 21.0 22.0 23.0 24,0 

HOURS FROM EPOCH 


Figure B4— Position differences between TRANE I Doppler orbit and optica! orbit 
for the 2-day arc, July 23-24, 1966. 
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RADIAL DIFFERENCES 

CROSS TRACK DIFFERENCES 

* ALONG TRACK DIFFERENCES 



14.0 15*0 16*0 17.0 18.0 19.0 20.0 21.0 22.0 

HOURS FROM EPOCH 


B5— Position differences between TRANET Doppler orbit and optical orbit 
for the 2-day arc, August 4-5, 1966. 





METERS METER 



ALONG TRACK DIFFERENCES 



Figure B5 (continued)— Position differences between TRANET Doppler orbit and optical orbit 

for the 2 -day arc, August 4-5, 1966. 




Appendix C 


Position Differences Between Doppler and Optically Determined 
Orbits for Selected Arc Lengths 

Figures are presented for 9 selected arc lengths showing plots of the satellite position dif- 
ferences between the optically determined and Doppler determined orbits over the span of the 
particular arc. 
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METERS METERS 













METE 



13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 21.0 22.0 

HOURS FROM EPOCH 


23.0 24.0 


Figure C3— Position differences between TRANET Doppler orbit and optical orbit 
for the 1-day arc, July 17, 0 through 24 hrs, 1966. 




MET 



CROSS TRACK DIFFERENCES 

— ALONG TRACK DIFFERENCES 






METERS 



RADIAL DIFFERENCES 

CROSS TRACK DIFFERENCES 

ALONG TRACK DIFFERENCES 


Figure C5— Position differences between TRANET Doppler orbit and optica! orbit 
for the 1/2 -day arc, July 17 , 0 through 12 hrs, 1966. 



RADIAL DIFFERENCES 

CROSS TRACK DIFFERENCES 

ALONG TRACK DIFFERENCES 


Figure C6— Position differences between TRANET Doppler orbit and optica! orbit 
for the 1/2 -day arc, July 18, 0 through 12 hrs, 1966. 
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HOURS FROM EPOCH 


Figure C7— Position differences between 7RANET Doppler orbit and optical orbit 
for the 1-1/2 -day arc, July 19, 0 hrs— July 20, 12 hrs, 1966. 






METERS METERS 



RADIAL DIFFERENCES 

CROSS TRACK DIFFERENCES 

ALONG TRACK DIFFERENCES 


Figure C8— Position differences between TRANET Doppler orbit and optical orbit 
for the 1 /2-day arc, July 19, 0 through 12 hrs, 1966. 



RADIAL DIFFERENCES 

CROSS TRACK DIFFERENCES 

ALONG TRACK DIFFERENCES 


Figure C9— Position differences between TRANET Doppler orbit and optical orbit 
for the 1/2 -day arc, July 20, 0 through 12 hrs, 1966. 






Appendix D 


Timing Biases and Range-Rate (Base Frequency) Biases 
in the Doppler Data 

Table D1 presents timing biases,* calculated for each pass of TRANET Doppler data in the pe- 
riod July 17-20, 1966. The timing biases were determined by fitting the residuals from each pass 
to the formula: 


AR = AB + AtR , 


where 

ar = residual (o-c) 
ab = zero set bias 
At = timing error 

AR = rate of change of the observation. 

Table D2 presents a summary of the (base frequency) range rate biases determined for each 
pass of TRANET Doppler data over the entire period of the study, July 9-26, July 31— August 7, 

1966. In each 2-day Doppler orbital solution, the range rate biases for every pass were dynamically 
determined along with the state vector. 


* These timing biases could possibly be attributed to hardware or to an orbital error or a combination of both. 



Table D1 


Timing Biases Found in the TRANET Doppler Data, July 17-20, 1966. 


Station 

Day 

Pass 

Timing Error 
(msecs.) 

Mean Timing 
Error ± 1 Standard 
Deviation 
(msecs.) 

(1966) 

Start 

End 

APLMND 

July 17 

3:28 

3:47 

-0.6 




5:36 

5:54 

-1.5 




7:40 

8:03 

1.8 




9:49 

10:14 

-7.3 




11:55 

12:23 

-4.1 




14:02 

14:26 

4.1 

0.6 ± 4.2 


July 18 

3:30 

3:52 

1.2 




5:38 

5:59 

1.4 




7:45 

8:08 

4.4 




14:06 

14:28 

6.6 


LASHM2 

July 17 

3:42 

4:09 

-2.2 




5:48 

6:16 

-1.8 




7:54 

8:20 

-1.8 




21:33 

21:51 

4.4 




23:34 

23:58 

5.5 



July 18 

1:40 

2:05 

1.9 

0.0 ± 3.2 



3:49 

4:14 

-0.5 




5:53 

6:20 

0.6 




7:58 

8:24 

-5.0 




' 21:38 

21:56 

-1.2 


ANCHOR 

July 17 

7:34 

7:52 

-1.7 




9:35 

9:58 

3.8 




11:40 

12:04 

-1.7 




13:43 

14:10 

4.8 




15:49 

16:15 

-1.0 




17:53 

18:17 

-5.4 

-0.8 ± 3.2 


July 18 

9:39 

10:03 

-3.3 




11:45 

12:09 

-0.4 




13:48 

14:15 

0.6 




15:53 

16:19 

-3.6 


LACRES 

July 17 

9:43 

10:01 

7.4 




11:51 

12:15 

-3.1 




13:57 

14:25 

5.6 



July 18 

5:29 

5:50 

1.6 




7:37 

7:57 

1.9 

-0.4 ± 5.5 



9:48 

10:07 

-0.1 




11:56 

12:21 

-7.3 




14:01 

14:29 

-0.6 




16:08 

16:32 

-9.3 


WAHIWA 

July 17 

15:59 

16:20 

9.3 




18:03 

18:31 

2.1 



July 18 

7:24 

7:43 

-3.9 

2.7 ± 5.5 



16:03 

16:26 

7.0 




18:07 

18:35 

-1.1 
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Table D1 (Continued) 


Station 

Day 

(1966) 

Pass 

Timing Error 
(msecs.) 

Mean Timing 
Error ± 1 Standard 
Deviation 
(msecs.) 

Start 

End 

APLMND 

July 19 

3:34 

3:57 

2.2 




5:41 

6:03 

0.5 




7:50 

8:14 

1.6 




9:58 

10:25 

~6.7 




12:04 

12:32 

4.9 




14:11 

14:32 

10.6 

0.4 ± 5.8 


July 20 

5:46 

6:08 

1.1 




7:55 

8:18 

1.7 




10:04 

10:30 

-9.3 




12:08 

12:35 

4.0 




14:16 

14:35 

-6.5 


LASHM2 

July 19 

3:55 

4:18 

0.1 




5:58 

6:25 

-3.6 




8:03 

8:27 

-4.8 




21:40 

22:00 

2.0 



July 20 

1:49 

2:15 

-1.1 

-2.1 ± 3.0 



3:55 

4:23 

-2.1 




6:02 

6:29 

-4.4 




8:07 

8:31 

-6.5 




21:45 

22:05 

1.9 


LACRES 

July 19 

5:33 

5:55 

1.6 




12:00 

12:26 

1.2 




14:08 

14 :34 

3.5 



July 20 

5:38 

5:59 

1.9 

2.9 ± 1.5 



7:47 

8:06 

4.6 




12:05 

12:31 

4.9 




14:10 

14:39 

1 

2.7 


WAfflWA 

July 19 

7:29 

7:47 

-3.1 




16:07 

16:31 

11.3 




18:12 

18:40 

2.4 

3.5 ± 6.3 


July 20 

16:11 

16:36 

8.5 




18:12 

18:44 

-1.8 


ANCHOR 

July 19 

1 7:41 

8:01 

-2.6 




9:43 

10:07 

-3.9 




11:47 

12:12 

3.4 




13:52 

14:19 

-0.5 




15:57 

16:24 

-2,7 




18:02 

18:24 

-4.4 



July 20 

7:45 

8:05 

-0.1 

-1.7 ± 2.7 



9:48 

10:11 

-2.6 




11:52 

12:18 

1.8 




13:56 

14:22 

0.5 




16:01 

16:27 

-5.1 




18:07 

18:28 

-4.0 
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Table D2 


Summary of Range Rate Biases Found in the TRANET Doppler Data. 


2 -Day Arc 
(1966) 

Station 

Number of 
Passes 

Mean Range Rate 
Bias ± 1 Standard 
Deviation 
(cm/sec) 

2 -Day Arc 
(1966) 

Station 

Number of 
Passes 

Mean Range Rate 
Bias ± 1 Standard 
Deviation 
(cm/sec) 

July 9-10 

ANCHOR 

8 

8.0 ± 2.3 

July 23-24 

ANCHOR 

11 

7.9 ± 3.4 


WAHIWA 

5 

7,7 ± 2.2 


WAHIWA 

7 

8.3 ± 3.1 


LACHES 

8 

6.0 ± 2.9 


LACRES 

8 

7.3 =fc 3.3 


APLMND 

11 

8.7 ±3.1 


APLMND 

11 

9.0 ± 3.3 


LASHM2 

12 

7.7 ± 3.2 


LASHM2 

7 

7.0 dr 2.1 

July 11-12 

ANCHOR 

11 

8.2 ± 3.4 






WAHIWA 

2 

11.1 ± 1.8 

July 25-26 

ANCHOR 

9 

8.7 ± 3.3 


LACHES 

10 

8.9 ± 2.1 


WAHIWA 

7 

9.3 =fc 2.6 


APLMND 

10 

10.6 ± 3.1 


LACRES 

5 

10.5 dr 1.7 


LASHM2 

11 

9.3 ± 2.3 


APLMND 

12 

11.5 ± 3.2 






LASHM2 

u 

8.8 ± 2.5 

July 13-14 

ANCHOR 

12 

7.8 ± 3.2 






WAHIWA 

4 

9.7 ± 0.4 






LACRES 

8 

7.5 rfc 2.8 

July 31 — 

ANCHOR 

11 

9.1 dr 2.4 


APLMND 

8 

10.1 ± 5.3 

August 1 

WAHIWA 

7 

9.7 ± 2.7 


LASHM2 

12 

8.5 ± 2.9 


LACRES 

4 

8.9 ± 3,1 






APLMND 

11 

10.8 ± 3.5 

July 15-16 

ANCHOR 

11 

8.3 ± 3.5 


LASHM2 

11 

9.7 ± 3.5 


WAHIWA 

4 

9.2 ± 2.7 






LACRES 

6 

7.7 ± 2.8 






APLMND 

11 

10.4 =fc 3.0 

August 2-3 

ANCHOR 

12 

8.7 =t 3.3 


LASHM2 

12 

8.0 ± 2.9 


WAHIWA 

6 

9.2 dr 3.7 






LACRES 

4 

9.0 d: 2.4 


ANCHOR 

10 

9.3 ± 2.2 


APLMND 

12 

11.3 =fc 3.5 

July 17-18 

WAHIWA 

5 

9.9 ± 1.5 


LASHM2 

10 

8.4 rfc 2.7 


LACRES 

9 

8.6 ± 2.7 






APLMND 

10 

10.2 ± 2.9 






LASHM2 

10 

9.2 ± 2.6 

August 4-5 

ANCHOR 

12 

7.9 dr 5.1 






WAHIWA 

0 



ANCHOR 

12 

8.5 ± 2.3 


LACRES 

10 

7.5 d= 4.5 

July 19-20 

WAHIWA 

5 

10.3 ± 2.0 


APLMND 

9 

9.2 ± 3.1 


LACRES 

7 

11.1 ± 1.2 


LASHM2 

12 

6.9 ± 2.5 


APLMND 

11 

10.4 ± 4.0 






LASHM2 

9 

8.3 ± 1.5 


- 1 







August 6-7 

ANCHOR 

10 

8.0 rfc 2.9 


ANCHOR 

9 

7.7 ± 3.0 


WAHIWA 

4 

8.5 ± 2.6 

July 21-22 

WAHIWA 

4 

9.3 ± 1.0 


LACRES 

5 

6.6 dr 4.5 


LACRES 

10 

8.5 ± 3.6 


APLMND 

12 

10.5 dr 3.4 


APLMND 

12 

9.9 ± 3.0 


LASHM2 

13 

7.7 rfc 2.5 


[ LASHM2 

8 

6.7 ± 0.5 
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Data as Determined in Varying Arc Length Solutions 


The two tables of this appendix present the results of the study to determine the independence 
of the range rate bias adjustment on the TRANET Doppler data with respect to orbital arc length. 
Table El presents the 44 Doppler passes in the 2-day period of July 17-18, 1966 which were used 
in this report and the range rate biases that were dynamically determined for each pass in differ- 
ing arc length solutions (2, 1-1/2, 1, and 1/2 day arcs). Table E3 presents the same information 
for the 44 Doppler passes used in the 2- day period of July 19-20, 1966. 


Table El 


TRANET Doppler Range Rate Bias Adjustment as Determined in Various 
Arc Length Solutions Over July 17-18, 1966. 


Adjusted 
R Bias 
(cm/sec) 

Mean R 

Bias ± 1 Standard 

Deviation 

(cm/sec) 

11.8 

I 

12.1 

12.3 

12.3 ± 0.5 

12.9 


10.8 


11.2 

11.0 

11.4 ± 0.7 

12.4 


8.3 


8.3 

8.2 ± 0.1 

8.1 


12.6 


12.6 


11.5 

12.3 ± 0.5 

12.1 


12.5 

| 


Pass for Which R Bias 
is Computed 


Station 


Day 


Start 


End 


Arcs Used 


Epoch 
(July, 1966) 


Day 


Hour 


Length 
in Days 


ANCHOR 


July 17 

7:34 8:03 


July 17 

9:35 10:03 


July 17 

11:40 12:08 


July 17 

13:43 14:12 


17 

17 

17 

17 

17 

17 

17 

17 

17 

17 

17 

17 


2 

1 - 1/2 

1 

1/2 

2 

1 - 1/2 

1 

1/2 

2 

1 - 1/2 

1 

2 

1 - 1/2 

1 - 1/2 

1 


17 

17 

17 


12 

0 

12 


1/2 












Table El (Continued) 


Pass for Which R Bi 
is Computed 

as 

Arcs Used 

Adjusted 
R Bias 
(cm/sec) 

Mean R 

Bias ± 1 Standard 
Deviation 
(cm/sec) 


Epoch 
(July, 1966) 

Length 
in Days 


Day 

Station 

Start 

End 

Day 

Hour 

ANCHOR 

Jub 

17 

17 

0 

2 

7.5 


(continued) 

15:49 


16:18 

17 

0 

1-1/2 

7.7 

7.1 ± 0.5 




17 

12 

1-1/2 

6.5 





17 

0 

1 

6.8 






17 

12 

1/2 

6.9 



July 

17 

17 

0 

2 

5.5 



17:53 


18:22 

17 

0 

1-1/2 

5.5 






17 

12 

1-1/2 

4.5 

4.8 ± 0.6 





17 

0 

1 

4.3 






17 

12 

1/2 

4.3 



July 

18 

17 

0 

2 

8.5 



9:38 


10:07 

17 

0 

1-1/2 

8.0 






17 

12 

1-1/2 

8.7 

8.8 ± 0.7 





18 

0 

1 

8.7 






18 

0 

1/2 

10.0 



July 

18 

17 

0 

2 

9.1 



11:44 


12:13 

17 

12 

1-1/2 

9.2 

9.2 ± 0.1 





18 

0 

1 

9.4 



July 

18 

17 

0 

2 

10.7 



13:47 


14:16 

17 

12 

1-1/2 

11.1 

11.1 ± 0.3 





18 

0 

1 

11.3 





18 

12 

1/2 

11.2 



July 

18 

17 

0 

2 

7.7 



15:52 


16:21 

17 

12 

1-1/2 

8.2 

8.2 ± 0.3 





18 

0 

1 

8.5 





18 

12 

1/2 

8.3 


WAH1WA 

July 

17 

17 

0 

2 

11.4 



15:58 


16:27 

17 

0 

1-1/2 

11.3 






17 

12 

1-1/2 

10.3 

10.8 ± 0.6 





17 

0 

1 

10.9 





* 

17 

12 

1/2 

10.0 



July 

17 

17 

0 

2 

9.9 



18:02 


18:31 

17 

0 

1-1/2 

10.0 






17 

12 

1-1/2 

8.9 

9.2 ± 0.8 





17 

0 

1 

8.8 






17 

12 

1/2 

8.2 
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Table El (Continued) 


Pass for Which R Bias 
is Computed 

Arcs Used 

Adjusted 
R Bias 
(cm/sec) 

Mean R 

Bias ± 1 Standard 
Deviation 
(cm/sec) 

Epoch 
(July, 1966) 

Length 
in Days 

Station 

Day 

Start 

End 

Day 

Hour 

WAfflWA 

July 18 

17 

0 

2 

8.1 


(continued) 

7:24 

7:53 

17 

0 

1-1/2 

7.9 





17 

12 

1-1/2 

8.0 

8.4 ± 0.7 




18 

0 

1 

8.3 





18 

0 

1/2 

9.6 



July 18 

17 

0 

2 

11.3 



16:02 

16:31 

17 

12 

1-1/2 

11.5 





18 

0 

1 

11.7 

11.3 ± 0.3 




18 

12 

1/2 

10.9 



July 18 

17 

0 

2 

8.9 



18:07 

18:36 

17 

12 

1-1/2 

9.2 

Q T i A O 




18 

0 

1 

9.5 

U.l ± U.O 




18 

12 

1/2 

8.9 


LACRES 

July 17 

17 

0 

2 

7.3 



9:42 

1 10:11 

17 

0 

1-1/2 

7.0 

7 O t A Q 




17 

0 

1 

7.3 

/ «o ± U.O 




17 

0 

1/2 

7.7 



July 17 

17 

0 

2 

5.9 



11:51 

12:20 

17 

0 

1-1/2 

5.8 

5*9 ± 0.1 




17 

0 

1 

5.9 



July 17 

17 

0 

2 

11.9 



13:57 

14:25 

17 

0 

1-1/2 

11.9 





17 

12 

1-1/2 

10.5 

11*4 ± 0.6 




17 

0 

1 

11.5 





17 

12 

1/2 

11.2 



July 18 

17 

0 

2 

12.1 



5:28 

5:57 

17 

0 

1-1/2 

11.6 





17 

12 

1-1/2 

12.0 

12.0 ± 0.3 




18 

0 

1 

12.2 

| 




18 

0 

1/2 

12.3 



July 18 

17 

0 

2 

11.1 



7:36 

8:05 

17 

0 

1-1/2 

10.4 





17 

12 

1-1/2 

11.4 

11.1 ± 0.4 




18 

0 

1 

11.4 





18 

0 

1/2 

11.2 



July -18 

17 

0 

2 

8.0 



9:47 

10:16 

17 

0 

1-1/2 

7.2 





17 

12 

1-1/2 

8.2 

7.9 ± 0.4 




18 

0 

1 

8.2 





18 

0 

1/2 

8.1 j 
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Pass for Which R Bias 
is Computed 


Adjusted 
R Bias 
(cm/sec) 

Mean R 

Bias ± 1 Standard 
Deviation 
(cm/sec) 

Epoch 
(July, 1966) 

Length 
in Days 

1 Station 

Day 

Start 

End 

Day 

Hour 

1 LACRES 

July 18 

17 

0 

2 

5.9 


(continued) 

11:55 

12:24 

17 

12 

1-1/2 

6.0 

6.0 ± 0.1 

1 



18 

0 

1 

6.1 


i 

July 18 

17 

0 

2 

9.7 



14:01 

14:30 

17 

12 

1-1/2 

9.9 

9.8 ± 0.3 

I 



18 


1 

10.1 


1 

! 



18 


1/2 

9.6 


| 

! 

July 18 

17 


2 

5.2 



16:08 

16:37 

17 


1-1/2 

5.5 

5.5 ± 0.2 

1 ' 



18 

0 

1 

5.8 





18 

12 

1/2 

5.5 


APLMND 

July 17 

17 

0 

2 

9.6 



j 3:28 

3:57 

17 

0 

1-1/2 

10.0 

10.4 ± 0.7 




17 

0 

1 

11.2 



. 


17 

0 

1/2 

10.7 



July 17 

17 

0 

2 

7.0 

5:36 6:05 

17 

0 

1-1/2 

7.0 


17 

0 

1 

8.0 


17 

0 

1/2 


July 17 

17 

0 

2 


7:40 8:09 

17 

0 

1-1/2 



17 

0 

1 

10.6 


17 

0 

1/2 

10.8 

July 17 

17 

0 

2 

5.6 

9:48 10:17 

17 

0 

lrl /2 

5.6 


17 

0 

1 

5.9 

i 

17 

0 

1/2 

6.7 

July 17 

17 

0 

2 

7.6 

11:54 12:23 

17 

0 

1-1/2 

7.7 


17 

0 

1 

7.5 

July 17 

17 

0 

2 

11.8 

14:01 14:30 

17 

0 

1-1/2 

12.1 


17 

12 

1-1/2 

10.5 


17 

0 

1 

11.3 

' 

17 

12 

1/2 

11.5 


7.4 ± 0.5 


10.4 ± 0.3 


6.0 ± 0.5 


7.6 db 0.1 


11.4 ± 0.6 
















Table El (Continued) 


Pass for Which R Bias 
is Computed 

Arcs Used 

Adjusted 
R Bias 
(cm/sec) 

. 

Mean R 

Bias ± 1 Standard 
Deviation 
(cm /sec) 

Epoch 
(July, 1966) 

Length 
in Days 

Station 

Day 

Start 

. _ 

End 

Day 

Hour 

APLMNB 

July 18 

17 

0 

2 

12.5 

I 

(continued) 

3:30 

3:59 

17 

0 

1-1/2 

12.1 

1 




17 

12 

1-1/2 

12.2 

12.2 ± 0.2 




18 

0 

1 

12.2 

1 




18 

0 

1/2 

12.0 



July 18 

17 

0 

2 

9.6 

j 


5:38 

6:07 

17 

0 

1-1/2 

8.9 





17 

12 

1-1/2 

9.6 

9.3 ± 0.3 




18 

0 

1 

9.5 

| 




18 

0 

1/2 

9.0 

j 


July 18 

17 

0 

2 

13.8 



7:45 

8:14 

17 

0 

1-1/2 

13.2 





17 

12 

1-1/2 

13.9 

13.7 ± 0.3 




18 

0 

1 

13.9 





18 

0 

1/2 

13.7 



July 18 

17 

0 

2 

. 13.9 



14:06 

14:34 

17 

12 

1-1/2 

14.1 

14.1 ± 0.2 




18 

0 

1 

14.3 





18 

12 

1/2 

14.1 


LASHM2 

July 17 

17 

0 

2 

8.0 



3:41 

4:10 

17 

0 

1-1/2 

8.1 

8.6 ± 0.7 




17 

0 

1 

9.3 





17 

0 

1/2 

8.9 



July 17 

17 

0 

2 

8.0 



5:47 

6:16 

17 

0 

1-1/2 

8.2 

8.5 ± 0.6 

1 



17 

0 1 

1 

8.9 


i 

i 



17 

0 

1/2 

9.2 

1 

' 


July 17 

17 

0 

2 

8.8 

| 


7:54 

8:23 

17 

0 

1-1/2 

9.1 

9.4 ± 0.7 




17 

0 

1 

9.4 





17 

0 

1/2 

1 

10.4 



July 17 

17 

0 

2 

13.9 



21:32 

22:01 

17 

0 

1-1/2 

13.8 





17 

12 

1-1/2 

12.9 

12.5 ± 1.5 




17 

0 

1 

10.9 

1 




17 

12 

1/2 

10.9 j 

5 

j 


July 17 

17 

0 

2 

‘ 13.0 



23:34 

24:03 

17 

0 

1-1/2 

13.0 

12.9 ± 0.2 




17 

12 

1-1/2 

12.6 
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Table El (Continued) 


Pass for Which R Bias 
is Computed 

Arcs Used 

Adjusted 
R Bias 
(cm/sec) 

Mean R 

Bias ± 1 Standard 
Deviation 
(cm/sec) 

Epoch 
(July, 1966) 

Length 
in Days 

Static n 

Day 

Start 

End 

Day 

Hour 

LASHM2 

July 18 

17 

0 

2 

9.2 


(continued) 

1:39 

2:08 

17 

0 

1-1/2 

8.7 





17 

12 

1-1/2 

8.7 

8.5 ± 0.7 




18 

0 

1 

8.3 





18 

0 

1/2 

7.4 



July 18 

17 

0 

2 

8.3 



3:48 

4:17 

17 

0 

1-1/2 

7.8 





17 

12 

1-1/2 

8.0 

7.8 db 0.4 




18 

0 

1 

7.7 





18 

0 

1/2 

7.3 











July 18 

17 

0 

2 

10.4 



5:52 

6:21 

17 

0 

1-1/2 

10.0 





17 

12 

1-1/2 

10.1 

10.2 ± 0.2 




18 

0 

1 

10.0 





18 

0 

1/2 

10.5 



July 18 

17 

0 

2 

7.8 



7:58 

8:27 

17 

0 

1-1/2 

7.6 





17 

12 

1-1/2 

7.7 

7.9 ± 0.5 




18 

0 

1 

7.6 





18 

0 

1/2 

8.8 



July 18 

17 

0 

2 

5.1 



21:37 

22:06 

17 

12 

1-1/2 

6.4 

iR 4- 1 1 




18 

0 

1 

7.5 

0.0 ± 1.1 




18 

12 

1/2 

7.3 
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Table E2 


TRANET Doppler Range Rate Bias Adjustment as Determined in Various Arc Length 

Solutions Over July 19-20, 1966. 


Pass for Which R Bias 
is Computed 

Arcs Used 


n 

Mean R 

Bias ± 1 Standard i 
Deviation 
(cm/sec) 

Vr 

loch 
, 1966) 

Length 

Adjusted 
R Bias 
(cm/sec) 

Station 

Day 

(July 

Start 

End 

Day 

Hour 

in .uays 

ANCHOR 

July 19 

19 

0 

2 

9.7 



7:41 

8:09 

19 

0 

1-1/2 

9.7 





19 

0 

1 

9.9 

10.0 ± 0.5 




19 

0 

1/2 

10.8 



July 19 

19 

0 

2 

7.7 



9:42 

10:11 

19 

0 

1-1/2 

7.7 





19 

0 

1 

7.6 

8.1 ± 0.9 




19 

0 

1/2 

9.5 



July 19 

19 

0 

2 

11.9 



11:46 

12:15 

19 

0 

1-1/2 

11.9 

11.8 ± 0.1 




19 

0 

1 

11.8 



July 19 

19 

0 

2 

9.5 



13:51 

14:20 

19 

0 

1-1/2 

9.6 





19 

12 

1-1/2 

8.6 

9.0 ± 0.5 




19 

0 

1 

8.9 





19 

12 

1/2 

8.4 



July 19 

19 

0 

2 

7.1 



15:56 

16:25 

19 

0 

1-1/2 

7.3 





19 

12 

1-1/2 

6.3 

6.7 ± 0.5 




19 

0 

1 

6.6 





19 

12 

1/2 

6.2 



July 19 

' 

19 

o 

2 

5.0 



18:02 

18:31 

19 

0 

1-1/2 

5.3 





19 

12 

1-1/2 

4.2 

4.7 ± 0.4 




19 

0 

1 

4.6 





19 

12 

1/2 

4.4 



July 20 

19 

0 

2 

11.1 



7:45 

8:14 

19 

0 

1-1/2 

12.0 





19 

12 

1-1/2 

10.6 

11.4 ± 0.7 




20 

0 

1 

11.0 





20 

0 

1/2 

12.4 



July 20 

19 

0 

2 

6.4 



9:47 

10:16 

19 

0 

1-1/2 

7.4 





19 

12 

1-1/2 

6.2 

6.9 ± 0.9 




20 

0 

1 

6.1 





20 

0 

1/2 

8.3 

I 
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Table E2 (Continued) 


Pass for Which R Bias 
is Computed 

Arcs Used 


Mean R 

Bias ± 1 Standard 
Deviation 
(cm/sec) 

Et 

Doch 
, 1966) 

Length 

Adjusted 
R Bias 
(cm/sec) 

Station 

Day 

(July 

Start 

End 

Day 

Hour 

in juay s 

LACRES 

July 19 

19 

0 

2 

12.7 



5:33 

6:02 

19 

0 

1-1/2 

12.6 





19 

0 

1 

13.1 

12.9 ± 0.2 




19 

0 

1/2 

13.0 



July 19 

19 

0 

2 

9.6 



12:00 

12:29 

19 

0 

1-1/2 

9.7 





19 

12 

1-1/2 

9.0 

8.9 ± 0.8 




19 

0 

1 

8.5 

| 




19 

12 

1/2 

7.8 



July 19 

19 

0 

2 

10.0 



14:07 

14:36 

19 

0 

1-1/2 

10.2 





19 

12 

1-1/2 

9.0 

9.4 ± 0.7 




19 

0 

1 

9.2 





19 

12 

1/2 

8.4 



July 20 

19 

0 

2 

12.0 



5:37 

6:06 

19 

0 

1-1/2 

12.8 





19 

12 

1-1/2 

11.4 

12.0 ± 0.5 




20 

0 

1 

11.9 





20 

0 

1/2 

12.0 



July 20 

19 

0 

2 

10.1 



7:46 

8:15 

19 

0 

1-1/2 

10.7 





19 

12 

1-1/2 

10.4 

10.7 ± 0.5 




20 

0 

1 

10.7 





20 

0 

1/2 

11.5 



July 20 

19 

0 

2 

11.9 



12:04 

12:33 

19 

12 

1-1/2 

12.5 





20 

0 

1 

12.7 

11.9 ± 0.9 




20 

12 

1/2 

10.6 



July 20 

19 

0 

2 

11.2 



14:10 

14:39 

19 

12 

1-1/2 

11.5 





20 

0 

1 

11.4 

11.1 ± 0.5 




20 

12 

1/2 

10.3 


APLMND 

July 19 

19 

0 

2 

14.3 



3:34 

4:03 

19 

0 

1-1/2 

13.9 





19 

0 

1 

14.5 

14.0 ± 0.5 




19 

0 

1/2 

13.4 

| 

| 


July 19 

19 

0 

2 

11.0 



5:41 

6:09 

19 

0 

1-1/2 

10.7 

I 




19 

0 

1 

10.6 

10.7 ± 0.3 

1 




19 

0 

1/2 

10.3 
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Table E2 (Continued) 




Pass for Which R Bias 
is Computed 


Station 


APPMNI) July 20 

(continued) 14:1(5 14:44 


Day 

Start End 


July 19 

3:55 4:24 


July 19 

5:57 6:26 


July 19 

8:02 ' 8:21 


July 19 

21:40 22:08 


July 20 

1:48 ‘ 2:17 


July 20 

3:55 4.24 


July 20 

0:02 6:30 


Table E2 (Continued) 
Arcs Used 


Epoch 
(July, 1906) 


Length 


Adjusted 
R Bias 
(cm/sec) 


Mean R 

Bias :i: 1 Standard 
Deviation 


Day 

Hour 




19 

0 

2 

6.5 

— 

19 

12 

1-3/2 

6.5 


20 

0 

1 

6.6 

6.4 ± 0.2 

20 

12 

1/2 

6.2 


19 

0 

2 

6.2 

| 

19 

0 

1-1/2 

6.0 


19 

0 

1 

5.4 

5.7 1 0.5 

19 

0 

1/2 

5.1 


19 

0 

2 

7.2 


19 

0 

1-1/2 

7.1 


19 

0 

1 

6.6 

7.1 i- 0.4 

19 

0 

1/2 

7.6 



8.8 ± 0,9 


19 

0 

2 

10.6 


19 

0 

1-1/2 

11.1 


19 

12 

1-1/2 

9.3 

10.4 -k 0.7 

19 

0 

1 

10.5 


19 

12 

1/2 

10.5 


19 

0 

2 

8.1 


19 

0 

1-1/2 

8.6 


19 

20 

12 

0 

1-1/2 

1 

7.9 

9.1 

8.5 ± 0.5 

20 

0 

1/2 

8.6 j 



9.2 ± 0,6 


8.4 ± 0.9 
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GLOBAL C-BAND RADAR NETWORK 
CALIBRATION UTILIZING GEOS-II 
SATELLITE 


by 

Ronald L. Brooks 

Wolf Research and Development Corporation 


H. Raymond Stanley 

National Aeronautics and Space Administration 


ABSTRACT 

A calibration project for the C-Band Radar Network is 
currently in progress at NASA/Wallops Station. The 
project objectives, unique data processing problems, 
and accomplishments are summarized. Recommendations 
are presented for future radar calibration projects. 


1. GEOS-II C-BAND PROJECT 
1 . 1 Introduction 

Prior to the launch of the GEOS-II satellite, it was realized 
that data from the worldwide C-Band radars were potentially of great 
benefit to geodesy. The C-Band radar sites not only have good geo- 
graphic distribution (Figure 1), but are capable of all-weather, 
horizon- to -'horizon tracking at rates as high as 40 observations per 
second . Further , the majority of these radars provide unambiguous 
ranging to 32,000 nautical miles. However, the accuracy of these 
radars had not been established due to the lack of a suitable cali- 
bration target at satellite distances . 

Since February of 1968 , both land-based and shipborne C-Band 
radars have successfully tracked a C-Band transponder on board the 
GEOS-II satellite . The primary goals of this GEOS-II Project are : 
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a) To provide geodetic scale for the interim NASA unified 
network . 

b) To better determine the absolute accuracy of instrumentation 
radar systems, develop refined methods of calibrating these 
systems, and improve the techniques employed in processing 
the associated data. 

c) To better determine the geodetic location of the C-Band 
radar sites and their intersite distances. 

d) To compare and correlate results obtained from other GEOS-B 
systems with those obtained by the C-Band system. 

e) To make generally available the results of both the C-Band 
system calibration and geodetic investigations. 

These objectives will assist specifically in meeting the GEOS Program 
objectives stated in the GEOS-A Mission Plan. 

1 . 2 Accomplishments to Date 

A team comprised of radar engineers and geodetic software special- 
ists from NASA/Wallops Station, Wolf Research and Development*, and 
RCA/Moores town** , has designed controlled radar studies involving data 
acquisition, reduction and analysis . The accomplishments to date 
include : 

a) Ascertained the range accuracy of properly operated and 

ground calibrated C-Band radars . Wallops Station's results 

indicate that the AN/FPS-16 and AN/FPQ-6 radars are capable 

of providing range data at satellite distances with an 

accuracy of 1 - 4 meters , with an RMS noise of approximately 

1,2 

one meter. 

b) Recovered a composite set of resonant geopotential coeffic- 
ients for GEOS- II , from C-Band radar data alone , which 
improved orbital accuracies ? 

* Contracts NAS6-1467 and NAS6-1628 with NASA/Wallops 

Subcontractor to Wolf Research and Development Corporation 
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c) Recovered geodetic positions of Apollo tracking ships 
utilizing data from ship-borne and land-based radars!! 

d) Utilized pulse doppler system of Wallops AN/FPQ-6 to 
generate range-rate-derived range measurements. The 
resulting ranges illustrate RMS noise values as low 
as 0.0006 meters from short arc orbital fits? 

e) Improved procedures for radar calibrations which are 
applicable to all radar systems. 

f) Recovered significant range, angle, and timing biases 
for several of the C-Band radars. 

g) Determined improved geodetic coordinates and intersite 
distances for participating radars? 

h) Filtered radar data to reduce storage requirements and 
computer time , but retains Statistical content . 

2. DATA PROBLEMS ENCOUNTERED AND RESOLVED 

It was of course anticipated that each radar site could have 
time or measurement biases , or that the station coordinates may be 
in error . In the course of this study , several other problem areas 
have been encountered and resolved. The most significant of these 

have been 

a) pulse-width mismatch 

b) range calibration target refraction 

c) radar operator/ data analyst interface 

2 . 1 Pulse-Width Mismatch 

Since the radar can only be calibrated us ing its own transmitter ' s 
output pulse width, a transponder- generated return pulse will introduce 
some pulse-width dependent range bias error . The range bias , B , 
caused by the pulse-width mismatch is approximately 7 ’ 8 

Bp = 0 . 5 (P^ - P^) (150 meters/ysec) 
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where 


= pulse-width received from transponder, in usee 
= radar calibration pulse-width, in ysec 

The effect of the pulse-width mismatch was particularly brought 
to light when it was observed that a 15-meter range bias existed be- 
tween the Bermuda AN/FPS-16 and AN/FPQ -6 radars with the AN/FPQ -6 
measuring longer. Upon our request , Bermuda personnel measured the 
radars' pulse widths and found that their nominal 0.5 usee pulse- 

9 

widths actually were 0.38 usee (AN/FPQ- 6 ) and 0.55 usee (AN/FPS-16). 

The pulsewidth mismatch, its cause and its effect, are discussed in 
considerably more detail in reference 7. 

2 . 2 Range Calibration Target Refraction 

Radar range measurements are made on a well -surveyed calibration 
target, both prior to and after each satellite track to obtain zero- 
set corrections for the track data. The multi-station reductions 
currently being performed have promoted the realization that signi- 
ficant range errors (range measurements are short) will result if 
refraction effects are not taken into account during the range target 
calibration. This error is nominally as large as 13 meters (Table 1) 
for the Carnarvon AN/FPQ -6 radar, due to the fact that their range cali- 
bration target is approximately 30 miles distant. For other radars, this 
error is normally 1-2 meters . 

2 . 3 Radar Operator/Data Analyst Interface 

As a result of the data used in the Network evaluation coming 
from many different radars and agencies, a communication problem 
developed. Information pertinent to the data reduction and analysis 
such as detailed calibration procedures at individual sites was not 
readily available. 
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Preliminary results concerning the individual radars have 
been sent to the participating radar sites and the radar system 
operators are now more cognizant of the objectives of the study 
being performed at NASA/Wallops Station. More information con- 
cerning on-site data corrections, radar timing systems, and site 
operating procedures is being provided. 

3. EXPERIMENT DESCRIPTIONS 
3 . 1 Achievable Range Accuracy 

Generally speaking, the C-Band radars are unmatched as geodetic 
instrumentation in terms of dependability, consistency, or precision 
The calibration project personnel at Wallops Island have been con- 
cerned with assessing or achieving accuracy . To achieve this goal, 
the following studies have been performed: 

a) Radar parameters have been varied during tracks of GEOS- II 
After short-arc orbital reductions, the residuals show 
effects on the range measurements of changing parameters 
such as pulsewidth , bandwidth and prf (Figure 2) . 

b) Collocated radar range residuals have been compared to 
assess relative accuracies. Five meters of the observed 
Wallops FPS-16 minus FPQ-6 range differences (shown in 
Figure 3) are attributable to the FPS-16 ' s use of a water 
tank for calibration , rather than a point source . 

c) Collocated laser and radar range residuals have been com- 
pared as an independent accuracy check (Table 2) . The 
average Wallops Island Collocation Experiment (WICE) 
Laser/FPQ- 6 range difference was only 0 . 5 meters . 

d) Unweighted WICE radar and laser range measurement were com 
pared with an optically -determined orbit . The results are 
shown in Table 3 . 
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DIFFERENCES BETWEEN AN/FPQ-6 RANGE § LASER RANGE 



206 


10 


WALLOPS ISLAND COLLOCATION 




OPTICAL SOLUTION 
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3 . 2 Recovery of Resonance Coefficients 


th 

A composite set of 13 order resonance coefficients have been 
recovered for GEOS-II using only radar data. Table 4 lists these 
recovered coefficients, along with a residual fit comparison for 
orbital solutions before and after coefficient recovery. 

3 . 3 Ship Radar Calibration and Position Recovery 

Dr. Martin is presenting a separate paper at this conference 
on this subject. 

3 . 4 Coherent Signal Processing (CSP) Ranges 

The CSP range-rates from the Wallops AN/FPQ-6 radar have been 
integrated to CSP ranges. Not only do these CSP ranges provide very 
low noise ranges, but even short bursts of CSP range data provide 
relatively accurate orbits. 

In a recent experiment , three different determinations of an 
orbit were made , each using only data from the same 20 second time 
span, including PCA. The first orbit was determined using the 
normal range , azimuth , and elevation measurements from the skin- track 
data . The second orbit determination was made using CSP ranges in 
place of the normal ranges used on orbit 1 . The third orbit was com- 
puted using filtered CSP ranges . The RMS of fit for each solution 
is shown in Figure 4 . 

Unweighted collocated laser data was used to assess the orbital 
accuracies of each solution. The laser residuals plotted in Figure 4 
clearly show a marked improvement in orbital determination with the 

use of CSP ranges . 

3 . 5 Improved Radar Site Geodetic Coordinates 

Improved station coordinates on a common datum have been determined 
for many of the C-Band sites. Mr . Leitao is presenting a separate 
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RECOVERED VALUES OF LUMPED 13TH 
ORDER RESONANCE COEFFICIENTS 




paper on this subject at this Conference. 


3 . 6 Data Filtering 

Due to the data gathering ability fo the C-Band radars, filtering 
programs have been developed to reduce the data storage requirements 
and computer time for orbital solutions ; however , the statistical 
content of the data is retained. 

Figure 5 compares filtered and unfiltered range residuals. A 
time-tag truncation problem in the NGSP data format is revealed . 


4. CONCLUSIONS AND RECOMMENDATIONS 

4 . 1 Software Development 

The radar hardware provides data with excellent inherent pre- 
cision. The proper software development and data analysis is necessary , 
however, to assess and achieve accuracy. 

4 . 2 Future C-Band Satellite Transponders 

It would not have been possible to recover the radar calibration 
information without a satellite such as GEOS- II . Due to the install- 
ation of new radars and the modification of existing radars , future 
satellites with C-Band transponders are required . It would be very 
desirable if these transponders could be manufactured such that they 
were pulse-width sensitive ; that is , they would respond with the same 
pulsewidth as the interrogation pulsewidth from the radar. It would 
be of benefit also , if future satellite transponders were coherent 
so that valid range-rate measurements would not have to be restricted 
to relatively low- signal skin tracking. 

4 . 3 Radar Site Operations 

The C-Band radars have all been built to government specifications ; 
these specifications include multiple switches for the radar operator 
to change to allow for a wide dynamic range of tracking conditions . 
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However , the calibrations performed at Wallops Island illustrate 
clearly that the majority of these switches affect the measurements. 

It is a pre-requisite to all radar calibration that all satellite 
tracking be performed in a consistent manner . 

4 . 4 Radar Data 

Steady progress has been made during this radar calibration 
project. We have become increasingly aware of the factors which 
affect the radar measurements , such as pulsewidth mismatch , con- 
figuration of the range calibration target , and calibration techni - 
ques. In view of this increased knowledge , we offer to provide fine - 
grain data corrections , where necessary, for that C-Band radar data 
which has been earlier distributed to others in the geode tic community . 
We also invite other geodetic investigators to utilize our data . 
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THE JUPITER CAMERA INTERCOMPARISON TEST 


by 

Frank G. Rawlins on 
David Wo Harris 
John H. Berbert 
John D. Oosterhout 
Goddard Space Flight Center 

From Novembers 1965 to May, 1966, a test was conducted at Jupiter, 
Florida to compare the various cameras used in the tracking of GEOS-I. 

The cameras used in this test were: 

Camera Focal Length 


1) 

SAO Baker -Nunn 

500 mm 

2) 

SAO K50 

914 rm 

3) 

PC- 1000 

1000 mm 

4) 

MOTS -40 

1000 mm 

5) 

PTH-100 

1000 mm 

6) 

MOTS -24 

610 mm 

7) 

BC4 

305 mm 


Three of these cameras , the MOTS -40, the PC-1000, and the PTH-100 have 
identical lens systems . 

The SAO cameras were operated by Smithsonian Astrophysical Observatory 
(SAO) personnel, the PC-1000 was operated by Air Force personnel and the 
remaining four cameras were operated under contract from GSFC by Bendix 
Field Engineering personnel. 

All of these cameras were located within thirty meters of each other 
to minimize ground parallax. The GEOS-I satellite was tracked by all of the 
cameras for a period of several months in order to accumulate a set of 
simultaneous observations for the intercomparison s tudy . 

The data taken was reduced by the various agencies using the techniques 
described in the "Preprocessing of Optical Satellite Observations" by 
Frank Donald Hotter. This is the Ohio State University Department of Geodetic 
Science Report No. 82 as revised in May, 1968. 

The data from the four cameras operated for GSFC was reduced by New 
Mexico State University under a contract from (SFC. The technique described 
in the above report for the MOTS -40 reduction was used on all of these 
cameras with appropriate modification for the two ballistic cameras . 
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Operational problems were encountered with the four cameras operated 
for GSFC which increased the difficulty of the reduction and led to a delay 
in reporting of the results. These problems will be detailed more fully 
in a report to be issued on the Jupiter experiment in the near future. 

A special problem occurred with the BC4 camera* Due to operator error 
a possible time error may have been introduced into the star calibration 
times of the plates taken by this camera. This error is in the range of 
0 to 300 milliseconds and varies from plate to plate. This leads to a varying 
bias in Right Ascension and is probably not recoverable. 

The data from this experiment did not provide as high a degree of 
simultaneous observations as was desired. Currently there are 65 passes 
in the Data Center from Jupiter which have two or more cameras observing 
the same pass. Of these 65 passes about 50 percent consist of two or three 
simultaneous camera sets. 

The results presented in this paper' are preliminary as all of the 65 
passes have not been analyzed. A more complete paper will be published 
in the near future giving a more extensive analysis of the Jupiter experiment. 
Forty-two (42) passes have been analyzed in part and the results are reported 
here. The data in this set for the two SAO cameras is the current data in 
the Data Center. The corrections for parallactic refraction and diurnal 
aberration were not applied to the SAO data in this preliminary analysis. 4 
These corrections are small. Paralactic refraction is on the order of 
0.150 seconds of arc and diurnal aberration is a maximum of 0.290 seconds 
of arc. These corrections will be applied to the SAO data in the final 
report. The data in this set for the PC-1000 was obtained from the Data • 
Center in December s 1969. Since that time the Aeronautical Chart and 
Information Center (ACIC) has resubmitted their Jupiter data to the Data 
Center. In this recent submission the total number of passes was markedly 
reduced. It was decided to use the data from the earlier set for this 
initial investigation in order to obtain a reasonable number of passes. 

A preliminary comparison was done with the two submissions (on four of the 
passes in both) and the changes were found to be small. The Right Ascension 
did not change. The change in Declination was about 0.003 seconds of arc. 

The data included for the four cameras operated for GSFC consists of 
passes that are in the Data Center currently. However^ for this preliminary 
investigation they were re- reduced using a reduction program which solves 
for the elements of interior orientation separately on each plate. In 
this reduction the correction for decentering distortion is not included. 

This was done as a special test to investigate some of the operational 
problems that occurred with these cameras * This will be discussed in the 
final report. 
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The initial approach adopted to compare the various cameras was to fit 
the observations from each pass to a common trajectory. All the cameras 
were given a weight of one in this curve fitting except the BC4 which was 
given a weight of zero because of the possible timing problem. The Noname 
orbit generator program was used to do the curve fitting and produce the 
res iduals . 


The table included presents the results of this effort. The residuals 
given are in Right Ascension and Declination. The statistic quoted in the 
table is the RMS of the residuals for each coordinate of the flash images 
on the plate. In most cases there were seven images present. The Right 
Ascension values are in seconds of arc and had been modified by the usual 
cosine of Declination factor. The Declination values are also in seconds 
of arc. At the bottom of the second page of the table the averages of 
each column in the table are given. 

As can be seen from the table, three of the cameras produced a larger 
number of observations for the 42 passes then the others . The breakdown 
is as follows : 


1) SAO Baker-Nunn 

2) MOTS -40 

3) PC-1000 

4) SAO K50 

5) MOTS -24 

6) PTH- 100 

7) BC4 


32 passes 
26 passes 
20 passes 
10 passes 
10 passes 
9 passes 
8 passes 


6 


The data contained in the table indicates the presence of three basic 
camera groupings. This can be seen from the average values given at the 
end of the table. There is* a four camera group consisting of the K50, the 
Baker-Nunn, the MOTS -40, and the PC-1000 whose average values are slightly 
smaller than those of the other groups . This group's average values are 
on the order of 1.5 seconds of arc. The two camera group , consisting of 
the MOTS -24 and the PTH- 100, has average values on the order of 2.0 seconds 
of arc. 


The BC4 camera ' s average values are higher than those of the other 
two groups being on the order of 3.3 seconds of arc. This is due in 
part to the possible timing problem. What part, if any, is due to the 
characteristics of the camera is being investigated. 

The analysis of the results of the Jupiter experiment is still in 
process . As was stated earlier a more detailed report will be published 
in the near future. 
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Geos “I and Geo® -IX MOTS Camera Operations , 
Data Analysis and Experiments 
David W. Harris 
Goddard Space Flight Center 


Operations 

MOTS camera operations on Geos -II were concluded December 22, 1969. 

By that time, the 40- inch focal length, 8- inch aperture cameras had compiled 
a substantial amount of data on both Geos -I and Geos-II. The percentage of 
successes to attempted operations and assignments for each MOTS station is 
provided in Table 1. 


Table 1 
Geos Results 


Geos -I Geos-II 


STATION 

MOTS 

OPERATIONS/ 

SUCCESSES/ 

OPERATIONS/ 

SUCCESSES/ 

NUMBER 

STATION 

ASSIGNMENTS 

OPERATIONS 

ASSIGNMENTS 

OPERATIONS 

1021 

B POINT 

537. 

45% 



1022 

FTMYSS 

727. 

65% 

51% 

79% 

1024 ■ 

W00M1A 

72% 

56% 



1038 

0R0RAL 

— 

— 

50% 

76% 

1025 

QUITOE 

6% 

15% 

3% 

15% 

1026 

LIMAFU 

1% 

0% 

0% 

0% 

1028 

SNTAGO 

47% 

43% 

53% 

70% 

1030 

MOJAVE 

70% 

77% 

71% 

88% 

1031 

JOBURG 

74% 

78% 

66% 

82% 

1043 

MADGAR 

33% 

69% 

27% 

67% 

1032 

NEWFLD 

19% 

23% 

6% 

33% 

1033 

COLEGE 

50% 

417. 



1036 

ALASKA 



57% 

44% 

1035 

WHKFLD 

22% 

43% 

15% 

61% 

1042/1037 

M5SMAN 

49% 

54% 

57% 

71% 

7036 

EDINBG 

60% 

76% 

547. 

86% 

7037 

COIMM 

49% 

747. 

46% 

79% 

7039 

BERMDA 

48% 

49% 

35% 

78% 

7040 

PURIC0 

677. 

59% 

52% 

79% 

7043 

GSFCPO 

35% 

50% 

MM 


7077 

GSFCNF 



22% 

75% 

7045 

DENVER 

56% 

527. 

55% 

88% 

7044 

CKVLLE 

18% 

1% 

_ 

«H» 

7072 

JUPTER 

59% 

38% 

MM. 



7075 

SUDBRY 

43% 

75% 

23% 

697. 

7076 

JAMACA 

48% 

75% 

42% 

75% 

1034/7034 

UND AKA 

53% ■ 

67% 

35% 

707. 

7078 

WA1M0T 


MM 

317. 

65% 

7079 

CAW0N 

- 

— 

43% 

62% 


TOTALS 

487. . 

59% 

40% 

721 


';V .'-jTV;' 
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For Geos -I, 9,500 plates were exposed and 5,576 were found to contain 
the satellite flashes. Results from more than 3,300 of these plates, 
representing more than 21,000 flashes are now on tape in the Data Center. 
Results for Geos -IX show that 7,891 plates were exposed, 5,714 of them 
containing flashes. By June 15, 1970, reductions of more than 2,500 
Geos -II plates were submitted to the Data Center. Approximately 1,000 
more plates are to be reduced. It is anticipated that all the MOTS data 
for these two spacecraft will be supplied to the Data Center prior to 
January, 1971. 

Validation 

Prior to submission to the Data Center, the reduced optical observations 
are processed through a validation program. The right ascensions and 
declinations from the reductions are compared to a reference orbit obtained 
from an independent tracking system. The gates applied, which determine 
whether or not an observation has passed the validation routine, is 
purposely set wide enough to remove only obvious outliers . A comprehensive 
report on the validation procedure is presented in the Goddard Space Flight 
Center (GSFC) document , NASA-GSFC GEOS -I MOTS Optical Validation Report . 
X-514-69-83 . 

In January.. 1969, while analyzing the Geos-I validation results , it 
became apparent that a subtle error existed in the MOTS reduction routine 
as employed by the reduction group at New Mexico State University (NMSU) . 

The error, occurring only when the camera was pointed at or near the 
station's local meridian (i.e. an azimuth at or near 0 or 180 degrees ) , 
varied in magnitude from only hundredths of an arc second to nearly 50 
seconds of arc depending upon the azimuth-elevation angles of the camera 
axis. The data containing the larger errors were flagged by the validation 
routine and were never submitted to the Data Center. • However , more than 
300 plates containing smaller errors of less than 10 arc seconds were 
provided to the Data Center before the problem was detected. Re-reductions , 
using corrected formulations were conducted and the new results submitted 
to the Data Center in June 1969. MOTS observation tapes provided by the 
center after October 9, 1969, contain the correct information and should 
be used in lieu of the earlier data. A list of the MOTS plates affected 
by this problem is provided in the above mentioned GSFC document X-514-69-83. 

As is so often the case, the error was inadvertently introduced while 
attempting to remove a nearly insignificant source of bias . It was realized 
that during a Geos operation, the sidereally driven MOTS camera images the 
reference stars at very slightly different positions due to differential 
atmospheric refraction during the 24 seconds of plate exposure. The method 
introduced to handle this condition was to compute the camera orientation 
angles , azimuth, elevation and tilt for the center flash of the sequence 
(equivalent to the measured center of the stellar images) and by differential 
form obtain the orientation angles to the earlier and later flashes. In 
this technique, when the camera had certain angles (near the local meridian), 
the trigonometric functions used in the computations rapidly approached 
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infinity and accuracy was lost. The program has been corrected by exchanging 
the differential form of orientation angle computation for an exact computa- 
tion for each desired flash time. 

MOTS Camera Sensitivity 

During the active tracking period of Geos- I, a camera sensitivity and 
tracking accuracy experiment was conducted at Jupiter, Florida. Six cameras 
were placed within 30 meters of the Baker- Nunn camera located there. By 
simultaneously observing the satellite flashes, the reduced observations 
may be compared and the relative accuracies of the cameras determined. A 
comprehensive report on this phase of the experiment is expected to be 
published by the GSFC shortly. 

A photometer was also situated at the Jupiter station during the 
experiment. By observing the stars in the region of the satellite passage 
with the photometer, a value for the 'atmospheric attenuation may be 
computed® Measurement of the diameters of the recorded flash images then 
allow for a determination of the sensitivity of the camera system. The 
results were published in the GSFC document, GEOS Tracking Camera Sensitivity 
Report , X-514-69-533. “ 

Briefly, the MOTS -40 camera results from Jupiter indicate that the 
camera- emulsion combination used appears to be more sensitive providing 
larger diameter images, than estimated prior to the launch of Geos -I. 

However, because of operational difficulties and changes made at the Jupiter 
station during the test period, it is possible that the larger images are 
due to focus problems and therefore the results may not be truly representative 
of all the MOTS cameras. Until such time when another test can be conducted 
under more uniform operational conditions, the a priori values used to 
determine the amount of illumination needed to obtain a reducible flash image 
should continue to be used. 

MOTS Camera Drive Evaluation 


The accuracy of a sidereally driven camera system has been an unanswered 
question. In connection with the evaluation of the accuracy of the MOTS camera 
data, an analysis of the camera drive system was undertaken . A report covering 
this subject has been prepared under the title Analyses of the MOTS Camera 
Drive , X-514-69-482. The report describes the tests conducted and presents 
the results for a number of Individual MOTS cameras. In general, the 
conclusions reached for the "average" MOTS camera reveals that; 

1. The camera exhibits a periodic drive oscillation on the order of 
nearly 90 seconds of time. 

2. For a short Geos type exposure, depending upon the position of the 
drive gear in its oscillation, the error in the flash positions will vary 
from nearly 0 arc seconds to a maximum of +: 2 arc seconds in right ascension. 

3. The error in declination caused by drive oscillations will be less 
than +1 second of arc. 

Before any operations are conducted on Geos-G, it is anticipated that 
a technique for smoothing out the drive oscillations will be incorporated 
into the camera system. By this means, errors in the reduced data due to 
the sidereal drive will be reduced to the order of +1 arc second in both 
coordinates • 
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